Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 6,933 G→A L343L (CTA→TTA)  yaaJ ← putative transporter YaaJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009136,9330GA97.5% 114.9 / ‑6.4 41L343L (CTA→TTA) yaaJputative transporter YaaJ
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (31/8);  total (32/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

AAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAAC‑GATGGAGCTGAAGGCAAACAGAATAACCACAAGGGTAACAAACTCA‑GCACCCCAGGAACCCA‑TT‑AGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTC‑CAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGAACCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTT  >  NC_000913/6689‑7176
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
aaaTTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGACGAGATTGTTTTCGGCGTAAATGTAGTTGGCAAC‑GATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCA‑GCACCCCAGGAACCCA‑TT‑AACACCCg                                                                                                                                                                                                                                                >  1:472681/1‑251 (MQ=255)
            gCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGCAGACTTGGCACGGTGCCGCCGGTGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAAC‑GATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCA‑GCACCCCAGGAACCCA‑TT‑AAAACCCGCATCGCCTTc                                                                                                                                                                                                                                      >  2:532575/1‑249 (MQ=255)
             ccATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAAC‑GATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCA‑GCACCCCAGGAACCCA‑TT‑AACACCCGCATCGCCTTCTgg                                                                                                                                                                                                                                   >  1:556621/1‑251 (MQ=255)
             ccATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAAC‑GATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCA‑GCACCCCAGGAACCCA‑TT‑AACACCCGCATCGCCTTCTgg                                                                                                                                                                                                                                   >  1:410889/1‑251 (MQ=255)
             ccATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAAC‑GATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCA‑GCACCCCAGGAACCCA‑TT‑AACACCCGCATCGCCTTCTgg                                                                                                                                                                                                                                   <  1:410701/251‑1 (MQ=255)
             ccATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAAC‑GATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCA‑GACCCCCAGGAACCCA‑TT‑AACACCCGCATCGCCTTCTgg                                                                                                                                                                                                                                   >  2:56726/1‑251 (MQ=255)
              cATGCAGGCCATTATGAGATCTGCCAGGTGCCACATCAGCGGAATACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAAC‑GATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCA‑GCACCCCAGGAACCCA‑TT‑AACACCCGCATCGCCTTCTGGa                                                                                                                                                                                                                                  <  2:62207/251‑1 (MQ=255)
                       cATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAAC‑GATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCA‑GCACCCCAGGAACCCA‑TT‑AACACCCGCATCGCCTTCTGGAt                                                                                                                                                                                                                                 >  1:230314/1‑243 (MQ=255)
                       cATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAAC‑GATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCA‑GCACCCCAGGAACCCA‑TT‑AACACCCGCATCGCCTTCTGGAt                                                                                                                                                                                                                                 <  2:230314/243‑1 (MQ=255)
                        aTTATGATATCTGCCAGATGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAAC‑GATGGAGCTGTAGGCAAACAGAATAGCCACGAGGGTAACCAACTCA‑GCACCCCAGGAACCCA‑TT‑AACACCCGCATCGCCTTCTGGATAAGCTGAAt                                                                                                                                                                                                                        >  2:560600/1‑251 (MQ=255)
                          tATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAAC‑GATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCA‑GCACCCCAGGAACCCA‑TT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATAc                                                                                                                                                                                                                      >  1:183429/1‑251 (MQ=255)
                                             cacaTCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAAC‑GATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCA‑GCACCCCAGGAACCCA‑TT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTC‑CAGCGGCATGTAGGt                                                                                                                                                                                                  >  1:38219/1‑251 (MQ=255)
                                             cacaTCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAAC‑GATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCA‑GCACCCCAGGAACCCA‑TT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTC‑CAGCGGCATGAAGGt                                                                                                                                                                                                  >  1:328806/1‑251 (MQ=255)
                                                aTCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAAC‑GATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCA‑GCACCCCAGGAACCCA‑TT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTC‑CAGCGGCATGTAGGTtgt                                                                                                                                                                                               >  1:632483/1‑251 (MQ=255)
                                                   aGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAAC‑GATGGAGCTGAAGACAAACAGAATAACCACGCGGGTAACAAACTCA‑GCACCCCAGGAACCCA‑TT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTC‑CAGCGGcctggagtgtgtgcc                                                                                                                                                                                            >  2:289806/1‑237 (MQ=255)
                                                                  taGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGACAACACCAGATAGAATTAGGGTTGGTGAGGCGTGAAAAGAAGAGAATGTTATAGGCGTAAATGTAGTGGGCAAC‑GCAGGAGCTGAAGACAAACAGAAGAACCACGAGGGAAATAAACACA‑GCACCCCAGGAACACA‑GT‑AGCACCACCATCGCCTTCTGGACAAGCTGAATCCCTAC‑CAGCGGCAGGTAGATTGGGCCGTTACCCGCCAgga                                                                                                                                                                              >  2:265600/2‑248 (MQ=255)
                                                                      tGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAAC‑GATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCA‑GCACCCCAGGAACCCA‑TT‑AACACCCGCATCGCCTTCTGGATAAGCTGACTACCTTC‑CAGCGGCATGTGGGTTGTGCCGGTACCCGCCAGTAatacc                                                                                                                                                                         >  2:284433/1‑249 (MQ=255)
                                                                                     ttGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAAC‑GATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCA‑GCACCCCAGGAACCCA‑TT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTC‑CAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGcc                                                                                                                                                          >  1:55406/1‑251 (MQ=255)
                                                                                       gCGAAGTTGCAGTCCCGCAAACCCCAGATTGTTTTGGGGTTGTCCGGGCAAAAAAAAAAGAGTTTGTTTGCGTCTTAAATGGAGTTGGCAAC‑GATGGAGCGGAGGCCAAAAAGAATACCCACGAGGGTAAAAAATTAA‑GCCCCCCAG‑ACCCCA‑TTTAACACCCGCTTCGCTTTTTGGAGAAGCTAATTCCCTTC‑CAGCGGCATTTGGGTTGTGCGGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGt                                                                                                                                                        <  1:316688/251‑1 (MQ=255)
                                                                                         gAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAAC‑GATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCA‑GCACCCCAGGAACCCA‑TT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTC‑CAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTAc                                                                                                                                                      >  1:419498/1‑251 (MQ=255)
                                                                                                cAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAAC‑GATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCA‑GCACCCCAGGAACCCA‑TT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTC‑CAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGAc                                                                                                                                               >  1:578512/1‑251 (MQ=255)
                                                                                                cAGATCCGCAAACACCAGATCGCTTTAGGGATGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAAC‑GATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCA‑GCACCCCAGGAAACCA‑TT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTC‑CAGCGGCATGTAGGTTCTGCCGTTACCCGCCAGGAATACCCGCATGgccacttccgtccagatgcc                                                                                                                                               >  2:369636/1‑233 (MQ=255)
                                                                                                  aacccgaaacctccAGAGCGCTTTGGGTTTTTTCGGGCGTAAAAAGAAGATTTTGTTTGCGGCAAAAATGTGGGTGG‑ACCAGAGGGAGCGGAAGCCAAACATAACACCCGGGAGGGAAACAAACCCA‑CCCCCCCAG‑AACCCAATA‑AACCCCCGCTTCGCCTTCTGGAAAAGCTAAAT‑CCTTCACAGCGGAATGTGGGTTGGCCCTTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCa                                                                                                                                             <  1:166230/239‑1 (MQ=255)
                                                                                                           acacCAGATCGCTTTAGGGTTGTTCCGGCGTAAAAAGAAGAGATTGCTTTCGGCGTAAATGTAGTTGGCAAC‑GATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCA‑GCACCCCAGGAACCCA‑TT‑AACAcc                                                                                                                                                                                                                                                  >  2:277168/1‑142 (MQ=255)
                                                                                                           acacCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAAC‑GATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCA‑GCACCCCAGGAACCCA‑TT‑AACAcc                                                                                                                                                                                                                                                  <  1:277168/142‑1 (MQ=255)
                                                                                                                aGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAAC‑GATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCA‑GCACCCCAGGAACCCA‑TT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTC‑CAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCa                                                                                                                                             >  1:51919/1‑237 (MQ=255)
                                                                                                                    cgttttAGGGTTGTGAAGGCAAAAAAAAAAGAGATTTTTTTCGGCGAAAATGTATTTGGAAAC‑AATGGACCTGAGGACAAACAAAAAAACCACGAGGTTAACAACCTCA‑CCCCCCCAGGACCCAA‑TA‑ACCCCCCGCATCGCCTTCGGAATAAGCGAAATCCCTCC‑CAGGGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAAt                                                                                                                             <  1:474523/246‑1 (MQ=255)
                                                                                                                        ttAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAAC‑GATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCA‑GCACCCCAGGAACCCA‑TT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTC‑CAGCGGCATGTAGGTTGTGcc                                                                                                                                                                                            >  1:494045/1‑182 (MQ=255)
                                                                                                                        gtAGGGTTGTTCATGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAAC‑GATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCA‑GCACCCCAGGATCCCA‑TT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTC‑CAGCGGCATGTAGGTTGTGcc                                                                                                                                                                                            <  2:494045/181‑1 (MQ=255)
                                                                                                                                        gTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAAC‑GATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCA‑GCACCCCAGGAACCCA‑TT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTC‑CAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCt                                                                                                       >  1:426086/1‑251 (MQ=255)
                                                                                                                                             aagaagAGATAGTTTTCGGAGTAAATGTGGTTGGCAAC‑GATGGCGCTGAAGACAAACCGAATAACCACGAGGGTAACAAACTCA‑GCACCCCAGGAACCCA‑TT‑AACACCCGCAACGCCTTCTGGATAAGCTGAATACCTTC‑CAGCGGCATGCAGGCTGTGCAGTTACCCGCCAGTAAAACCAGCAGGGCGCTTTCGGTACAGATGCGCAGGGATTCGGTAAACATGCCACACATATGGCCCATCCCTTgccc                                                                                                  >  2:423325/1‑249 (MQ=255)
                                                                                                                                                        ttttcgggcgaaaaTGTAGTTGGAAAC‑GATGGAGCTAAAGGAAAACAGAATACCCTCGAGGGTAAAAAAC‑CACGCCCCCCAGAACCCCA‑TT‑AACCCCCGCATCGCCTTCTGGATAAGCTAAATCCCTTC‑CAGCGGTATGTAGGTTGGGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGc                                                                                        <  1:182574/239‑1 (MQ=255)
                                                                                                                                                                 gtAAATGTAGTTGGCAAC‑GATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCA‑GCACCCCAGGAACCCA‑TT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTC‑CAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAgg                                                                              >  1:108881/2‑251 (MQ=255)
                                                                                                                                                                                     aTGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCA‑GCACCCCAGGAACCCA‑TT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTC‑CAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAgg                                                                                   >  1:565976/1‑227 (MQ=255)
                                                                                                                                                                                      tGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCA‑GCACCCCAGGAACCCA‑TT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTC‑CAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGAGCTGCCGGGTGCGGAGGCCAGGACGCCCCTGCCGCTGCcgcg                                                          >  2:573852/1‑251 (MQ=255)
                                                                                                                                                                                              aaGGCAAACAGAATAACCACGAGGGTAACAAACTCA‑GCACCCCAGAAACCCA‑TT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTC‑CAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAAACCCTTGAGCTGCCGGATGCGGAGGCCAGGAAGCCGCTGCCGCTGCCGCGTTTGGCgg                                                  >  1:517001/1‑250 (MQ=255)
                                                                                                                                                                                                    aaCAGAATAACCACGAGGGTAACAAACTCA‑GCACCCCAGGAACCCA‑TT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTC‑CAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCCGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGaccc                                            >  1:62464/1‑248 (MQ=255)
                                                                                                                                                                                                                          cAAACTCA‑GCACCCCAGGAACCCA‑TT‑AACACCCGCATCGCCTTCTGGCTAAGCTGAATACCTTC‑CAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACCGATGACCCGGGCGTCGATAAAAAAGCCAATCATCTGGACAATCCCTTACGCTGCCGGATGCGGAGGCCCGGAAGCCGCTGCCGCTGCCGCGTTTGGCGTCGCCCCCATTACCGCCCCATTGGAAAAc                      >  2:57250/1‑251 (MQ=255)
                                                                                                                                                                                                                                cA‑GCACCCCAGGCACCCA‑TT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTC‑CAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTAGACCCCATTCCCGCCTCATTGGAAAACATACTg                >  1:14699/1‑251 (MQ=255)
                                                                                                                                                                                                                                               cccA‑TT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTC‑CAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGtt  >  1:624936/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 cA‑TT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTC‑CAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGATGAAAACCGtt  >  1:642856/1‑249 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
AAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAAC‑GATGGAGCTGAAGGCAAACAGAATAACCACAAGGGTAACAAACTCA‑GCACCCCAGGAACCCA‑TT‑AGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTC‑CAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGAACCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTT  >  NC_000913/6689‑7176

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑