Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 9,839 G→A P178P (CCG→CCA)  mog → molybdopterin adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009139,8390GA96.4% 163.3 / ‑2.6 57P178P (CCG→CCA) mogmolybdopterin adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (19/34);  total (20/37)
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.35e-01

GCCTGGCTTTGGTGAACAGATGCGCCAGATCAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTGGGCGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGAGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGAT  >  NC_000913/9593‑10077
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              
gCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGaa                                                                                                                                                                                                                                             <  1:516798/250‑1 (MQ=255)
    ggCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGt                                                                                                                                                                                                                                        <  1:179175/251‑1 (MQ=255)
        ttGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTg                                                                                                                                                                                                                                      <  1:540519/249‑1 (MQ=255)
                     gcggcAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGGCAGGTTGGTGTGATGCGCAAACAGGCGCTGATCCGTAGCTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCGGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTTCCGTACTGCATTCAGTTGCTGTAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGa                                                                                                                                                                                                                       <  1:626855/247‑1 (MQ=255)
                           gataaGCCTGCATTTTGCACCAACTGCGATCCTTTCGCGTCAGGTGGGTGTGATTCGCAAACAGGCGCTGATCCTTACCGTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTTCATTCAGACCGAagaggg                                                                                                                                                                                                                 <  2:326766/247‑3 (MQ=255)
                              cAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTGGGCGTGATTCGCAAACAGGCGCCGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGAGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCa                                                                                                                                                                                                               >  1:387338/1‑250 (MQ=255)
                                     cATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGAcgc                                                                                                                                                                                                         <  1:167814/249‑1 (MQ=255)
                                     cATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGGGGTAGCTGCATTCAGACCGAAGAGTGCAAGAcgc                                                                                                                                                                                                         >  2:167814/1‑249 (MQ=255)
                                      gttttGTACCAACTGCGATCCTTTCGAGTCAGGTGGGTGTGATGCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTAGTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATGTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGAc                                                                                                                                                                                                      <  2:17926/250‑1 (MQ=255)
                                      aTTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGAc                                                                                                                                                                                                      <  1:47024/251‑1 (MQ=255)
                                           gTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGtt                                                                                                                                                                                                   >  1:267934/1‑249 (MQ=255)
                                            tACCAACTGCGATCCTTTCGCGTCAGGTTGGTATGATTCGCAAACGGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGCAGGACGCTGCGGGTAACGTTGAGGGACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGCAGGGCCATAC‑GTAGAACCGGCACAAGAAGTGGTTGCTGCATTAAGACCTAAGAGTGCAAGCCGCGACggtac                                                                                                                                                                                                 >  2:351658/1‑246 (MQ=255)
                                             aCCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGc                                                                                                                                                                                                <  1:399955/250‑1 (MQ=255)
                                               cAACTGCGATCCTTTCGCGGCAGGTTGGTGTGATTCGCAAACAGGCGCTGAGCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGaa                                                                                                                                                                                             <  1:98341/251‑1 (MQ=255)
                                                 aCTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATa                                                                                                                                                                                           <  1:187018/251‑1 (MQ=255)
                                                 aCTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATa                                                                                                                                                                                           <  1:242817/251‑1 (MQ=255)
                                                    gCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGTTGATCCTTAACTTACCCGGTCAGTCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATaa                                                                                                                                                                                          <  1:553062/249‑1 (MQ=255)
                                                         ccTTTCGCGTCAGGTTGGTGTGATACGAAAAAAGGCGCTGAGCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATc                                                                                                                                                                                   <  2:367728/251‑1 (MQ=255)
                                                            ttCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATcccc                                                                                                                                                                                >  2:482374/1‑251 (MQ=255)
                                                             tCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGACGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATcccc                                                                                                                                                                                <  2:231068/250‑1 (MQ=255)
                                                                 gTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAg                                                                                                                                                                           <  1:289267/251‑1 (MQ=255)
                                                                    gggtgggtgggagtCGCAAACAGGCGCTGATCCTTAACTGTCCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATGCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCgg                                                                                                                                                                        <  2:139556/238‑1 (MQ=255)
                                                                          gagTGATGCGCACACAGTCGCTTATCCTTAACGTACCCGGGCAGTCTAATTCTATTAAAGAGACGCTCGAAGGTGTGAAGGACCCTGCGTGTAACGTTGTGGTACACAGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCCTCC‑GTTGAAACGGCACCCGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAATTCGCGAGGTTTGCGAATAAAAAAATCCCCCCGCGCGGGGGGAt                                                                                                                                                                  <  2:253074/249‑1 (MQ=255)
                                                                             tGATTCGCAAACAGGCGCTGATCCTTACCTTACCCGGTCCGCCGCAGTATATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtctc                                                                                                                                                               >  1:569337/1‑251 (MQ=255)
                                                                                             gcTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGACGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTggg                                                                                                                                                >  1:446762/1‑250 (MQ=255)
                                                                                                      tAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCa                                                                                                                                        <  1:595333/249‑1 (MQ=255)
                                                                                                        actgaCCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGTATTCACCAATCg                                                                                                                                    <  1:476921/247‑1 (MQ=255)
                                                                                                            aCCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCa                                                                                                                                 <  1:45602/250‑1 (MQ=255)
                                                                                                               catctaGCCGAATGCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCAGTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTGAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCGCCAATCGGCAGAAc                                                                                                                             <  2:637675/246‑1 (MQ=255)
                                                                                                                ggTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACg                                                                                                                            <  2:569337/251‑1 (MQ=255)
                                                                                                                ggTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGGGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGCCGGTAGCGAATAAAAAAATCCACCCGAGCGGGGGGATCCCACAACAATTACTGGGATTGCCCAATCGTCAGaccg                                                                                                                            >  2:103787/1‑248 (MQ=255)
                                                                                                                    aGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAAc                                                                                                                             >  1:637675/1‑246 (MQ=255)
                                                                                                                        gAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGAcc                                                                                                                    <  2:267934/251‑1 (MQ=255)
                                                                                                                          aggcTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAAGTAGTGGGATTCACCAATCgg                                                                                                                                   <  2:565018/231‑1 (MQ=255)
                                                                                                                          aGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCgg                                                                                                                                   >  1:565018/1‑234 (MQ=255)
                                                                                                                                   aaGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGt                                                                                                         <  1:442312/251‑1 (MQ=255)
                                                                                                                                     gagaCGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTc                                                                                                       <  1:16652/251‑1 (MQ=255)
                                                                                                                                             ggAAGGTGTGAAGGACGCTGAGGGTAACGGTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGTCGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTTCtt                                                                                                <  2:387338/250‑1 (MQ=255)
                                                                                                                                                   tgtgAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGACGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCAt                                                                                         >  1:98634/1‑251 (MQ=255)
                                                                                                                                                                           tgtgGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAAAAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTACGTACTTCACCCATCGCCAGATAGATTGCACTGGCAcc                                                                 >  2:67067/1‑251 (MQ=255)
                                                                                                                                                                             tgGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCAccac                                                               >  1:108100/1‑249 (MQ=255)
                                                                                                                                                                                 acacGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGa                                                            <  1:482374/250‑1 (MQ=255)
                                                                                                                                                                                    cGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGATCGTTCCGTACTTCACCCATCGCCAGATAGATTGCACTGGCAccaccgatc                                                          >  2:637237/1‑242 (MQ=255)
                                                                                                                                                                                            ttGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCTGGGGGATCTCAAAACAATTATTGGGATTCACCAATCGGCAGAACGGTGCTACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATc                                                 <  1:263375/250‑1 (MQ=255)
                                                                                                                                                                                              gCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCa                                               <  1:228040/250‑1 (MQ=255)
                                                                                                                                                                                                  gCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCgg                                          <  1:268313/251‑1 (MQ=255)
                                                                                                                                                                                                           ccTGCATTCAGGTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGga                                 <  1:359739/250‑1 (MQ=255)
                                                                                                                                                                                                           acTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtct                                                                                                                                                                <  2:546511/124‑1 (MQ=255)
                                                                                                                                                                                                           acTGCATTCAGTTGCTGGAAGGGCAATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtct                                                                                                                                                                >  1:546511/1‑124 (MQ=255)
                                                                                                                                                                                                            cTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCCTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCCGAACGGTGCGACCAAACTGCTCGTTCAGCACTTCACCCATCGCCAGATAGATTGCACTGGCACACAAGATCGGACCAATCCCGCCGGCAAAGTGga                                 >  2:359739/1‑250 (MQ=255)
                                                                                                                                                                                                            cTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGgat                                >  1:271460/1‑251 (MQ=255)
                                                                                                                                                                                                                          tGGAAGGGCCATACCGTT‑AAACGGCAACAGAAGTGGGTTCTGCAATCAGACCGCAGAGTGCAAGACGCGACGTTAGCTAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGCATGCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTAAGTACTTCACCCACCGCCAGATAGATTGCACTGTCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTa                   <  2:446762/250‑1 (MQ=255)
                                                                                                                                                                                                                            gAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTAc                  >  1:48425/1‑249 (MQ=255)
                                                                                                                                                                                                                               gggCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACTGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCa             <  1:614304/251‑1 (MQ=255)
                                                                                                                                                                                                                                    aTAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGtt        >  1:653212/1‑251 (MQ=255)
                                                                                                                                                                                                                                       c‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGAGCGACCAAACTGCTCGTTCAGCACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCgacactgtaacc     >  1:669897/1‑240 (MQ=255)
                                                                                                                                                                                                                                            gAAACGGCACTCGAAGGGGTTGAGGCAATCAGACCGATGAGTGCAAGACGCGTCGTTATCGAATATAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCAACCAGCGCCAGATAGATTGCACTGGCACCACAGAGCAGATCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCTATGTTAGCGAt  <  2:240856/250‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              
GCCTGGCTTTGGTGAACAGATGCGCCAGATCAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTGGGCGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGAGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGAT  >  NC_000913/9593‑10077

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