Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 9,851 A→T A182A (GCA→GCT)  mog → molybdopterin adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009139,8510AT95.0% 178.9 / ‑0.6 61A182A (GCA→GCT) mogmolybdopterin adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base T (22/35);  total (23/38)
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.84e-01

GCGCCAGATCAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTGGGCGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGAGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAG  >  NC_000913/9614‑10095
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                    
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:242817/251‑1 (MQ=255)
                            aCTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATa                                                                                                                                                                                                             <  1:187018/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)
                                                                                                                                                                                                                               accagAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAg                                                                                                                       >  1:606423/5‑143 (MQ=255)
                                                                                                                                                                                                                               accagAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAg                                                                                                                       <  2:606423/139‑1 (MQ=255)
                                                                                                                                                                                                                               accagAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAACGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCcagc             >  1:581663/5‑249 (MQ=255)
                                                                                                                                                                                                                                  agAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGc        >  1:110624/2‑251 (MQ=255)
                                                                                                                                                                                                                                      gTGGTTGCTGCAGTCAGACCGAAGAGGGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGTGGGATCTCAAAACAATTATGGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGTGCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACATATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAAc    <  2:271460/251‑1 (MQ=255)
                                                                                                                                                                                                                                       gggTTGCAGCATTCAGAACGAAGAGTGCAAGACGCGACGTTTGCGGATAAAAAAATCCCCCCGAGCGGGTGGATCTCAAACCAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCGAGCAGGGCAAACa   <  2:340916/250‑1 (MQ=255)
                                                                                                                                                                                                                                        ggTTGCTGCATTCAGACCGATGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAATTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAg  <  1:351658/251‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                    
GCGCCAGATCAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTGGGCGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGAGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAC‑GTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAG  >  NC_000913/9614‑10095

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: ‑