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→GCTmogmolybdopterin 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

GCGCCAGATCAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTGGGCGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGAGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAG  >  NC_000913/9614‑10095
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                    
gcggcAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGGCAGGTTGGTGTGATGCGCAAACAGGCGCTGATCCGTAGCTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCGGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTTCCGTACTGCATTCAGTTGCTGTAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGa                                                                                                                                                                                                                                         <  1:626855/247‑1 (MQ=255)
      gataaGCCTGCATTTTGCACCAACTGCGATCCTTTCGCGTCAGGTGGGTGTGATTCGCAAACAGGCGCTGATCCTTACCGTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTTCATTCAGACCGAagaggg                                                                                                                                                                                                                                   <  2:326766/247‑3 (MQ=255)
         cAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTGGGCGTGATTCGCAAACAGGCGCCGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGAGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCa                                                                                                                                                                                                                                 >  1:387338/1‑250 (MQ=255)
                cATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGAcgc                                                                                                                                                                                                                           <  1:167814/249‑1 (MQ=255)
                cATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGGGGTAGCTGCATTCAGACCGAAGAGTGCAAGAcgc                                                                                                                                                                                                                           >  2:167814/1‑249 (MQ=255)
                 gttttGTACCAACTGCGATCCTTTCGAGTCAGGTGGGTGTGATGCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTAGTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATGTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGAc                                                                                                                                                                                                                        <  2:17926/250‑1 (MQ=255)
                 aTTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGAc                                                                                                                                                                                                                        <  1:47024/251‑1 (MQ=255)
                      gTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGtt                                                                                                                                                                                                                     >  1:267934/1‑249 (MQ=255)
                       tACCAACTGCGATCCTTTCGCGTCAGGTTGGTATGATTCGCAAACGGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGCAGGACGCTGCGGGTAACGTTGAGGGACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGCAGGGCCATACGTAGAACCGGCACAAGAAGTGGTTGCTGCATTAAGACCTAAGAGTGCAAGCCGCGACggtac                                                                                                                                                                                                                   >  2:351658/1‑246 (MQ=255)
                        aCCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGc                                                                                                                                                                                                                  <  1:399955/250‑1 (MQ=255)
                          cAACTGCGATCCTTTCGCGGCAGGTTGGTGTGATTCGCAAACAGGCGCTGAGCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGaa                                                                                                                                                                                                               <  1:98341/251‑1 (MQ=255)
                            aCTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATa                                                                                                                                                                                                             <  1:242817/251‑1 (MQ=255)
                            aCTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATa                                                                                                                                                                                                             <  1:187018/251‑1 (MQ=255)
                               gCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGTTGATCCTTAACTTACCCGGTCAGTCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATaa                                                                                                                                                                                                            <  1:553062/249‑1 (MQ=255)
                                    ccTTTCGCGTCAGGTTGGTGTGATACGAAAAAAGGCGCTGAGCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATc                                                                                                                                                                                                     <  2:367728/251‑1 (MQ=255)
                                       ttCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATcccc                                                                                                                                                                                                  >  2:482374/1‑251 (MQ=255)
                                        tCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGACGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATcccc                                                                                                                                                                                                  <  2:231068/250‑1 (MQ=255)
                                            gTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAg                                                                                                                                                                                             <  1:289267/251‑1 (MQ=255)
                                               gggtgggtgggagtCGCAAACAGGCGCTGATCCTTAACTGTCCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATGCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCgg                                                                                                                                                                                          <  2:139556/238‑1 (MQ=255)
                                                     gagTGATGCGCACACAGTCGCTTATCCTTAACGTACCCGGGCAGTCTAATTCTATTAAAGAGACGCTCGAAGGTGTGAAGGACCCTGCGTGTAACGTTGTGGTACACAGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCCTCCGTTGAAACGGCACCCGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAATTCGCGAGGTTTGCGAATAAAAAAATCCCCCCGCGCGGGGGGAt                                                                                                                                                                                    <  2:253074/249‑1 (MQ=255)
                                                        tGATTCGCAAACAGGCGCTGATCCTTACCTTACCCGGTCCGCCGCAGTATATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtctc                                                                                                                                                                                 >  1:569337/1‑251 (MQ=255)
                                                                        gcTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGACGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTggg                                                                                                                                                                  >  1:446762/1‑250 (MQ=255)
                                                                                 tAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCa                                                                                                                                                          <  1:595333/249‑1 (MQ=255)
                                                                                   actgaCCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGTATTCACCAATCg                                                                                                                                                      <  1:476921/247‑1 (MQ=255)
                                                                                       aCCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCa                                                                                                                                                   <  1:45602/250‑1 (MQ=255)
                                                                                          catctaGCCGAATGCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCAGTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTGAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCGCCAATCGGCAGAAc                                                                                                                                               <  2:637675/246‑1 (MQ=255)
                                                                                           ggTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACg                                                                                                                                              <  2:569337/251‑1 (MQ=255)
                                                                                           ggTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGGGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGCCGGTAGCGAATAAAAAAATCCACCCGAGCGGGGGGATCCCACAACAATTACTGGGATTGCCCAATCGTCAGaccg                                                                                                                                              >  2:103787/1‑248 (MQ=255)
                                                                                               aGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAAc                                                                                                                                               >  1:637675/1‑246 (MQ=255)
                                                                                                   gAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGAcc                                                                                                                                      <  2:267934/251‑1 (MQ=255)
                                                                                                     aggcTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAAGTAGTGGGATTCACCAATCgg                                                                                                                                                     <  2:565018/231‑1 (MQ=255)
                                                                                                     aGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCgg                                                                                                                                                     >  1:565018/1‑234 (MQ=255)
                                                                                                              aaGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGt                                                                                                                           <  1:442312/251‑1 (MQ=255)
                                                                                                                gagaCGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTc                                                                                                                         <  1:16652/251‑1 (MQ=255)
                                                                                                                        ggAAGGTGTGAAGGACGCTGAGGGTAACGGTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGTCGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTTCtt                                                                                                                  <  2:387338/250‑1 (MQ=255)
                                                                                                                              tgtgAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGACGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCAt                                                                                                           >  1:98634/1‑251 (MQ=255)
                                                                                                                                                      tgtgGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAAAAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTACGTACTTCACCCATCGCCAGATAGATTGCACTGGCAcc                                                                                   >  2:67067/1‑251 (MQ=255)
                                                                                                                                                        tgGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCAccac                                                                                 >  1:108100/1‑249 (MQ=255)
                                                                                                                                                            acacGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGa                                                                              <  1:482374/250‑1 (MQ=255)
                                                                                                                                                               cGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGATCGTTCCGTACTTCACCCATCGCCAGATAGATTGCACTGGCAccaccgatc                                                                            >  2:637237/1‑242 (MQ=255)
                                                                                                                                                                       ttGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCTGGGGGATCTCAAAACAATTATTGGGATTCACCAATCGGCAGAACGGTGCTACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATc                                                                   <  1:263375/250‑1 (MQ=255)
                                                                                                                                                                         gCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCa                                                                 <  1:228040/250‑1 (MQ=255)
                                                                                                                                                                             gCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCgg                                                            <  1:268313/251‑1 (MQ=255)
                                                                                                                                                                                      ccTGCATTCAGGTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGga                                                   <  1:359739/250‑1 (MQ=255)
                                                                                                                                                                                      acTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtct                                                                                                                                                                                  <  2:546511/124‑1 (MQ=255)
                                                                                                                                                                                      acTGCATTCAGTTGCTGGAAGGGCAATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtct                                                                                                                                                                                  >  1:546511/1‑124 (MQ=255)
                                                                                                                                                                                       cTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCCTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCCGAACGGTGCGACCAAACTGCTCGTTCAGCACTTCACCCATCGCCAGATAGATTGCACTGGCACACAAGATCGGACCAATCCCGCCGGCAAAGTGga                                                   >  2:359739/1‑250 (MQ=255)
                                                                                                                                                                                       cTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGgat                                                  >  1:271460/1‑251 (MQ=255)
                                                                                                                                                                                                     tGGAAGGGCCATACCGTTAAACGGCAACAGAAGTGGGTTCTGCAATCAGACCGCAGAGTGCAAGACGCGACGTTAGCTAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGCATGCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTAAGTACTTCACCCACCGCCAGATAGATTGCACTGTCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTa                                     <  2:446762/250‑1 (MQ=255)
                                                                                                                                                                                                       gAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTAc                                    >  1:48425/1‑249 (MQ=255)
                                                                                                                                                                                                          gggCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACTGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCa                               <  1:614304/251‑1 (MQ=255)
                                                                                                                                                                                                               aTACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGtt                          >  1:653212/1‑251 (MQ=255)
                                                                                                                                                                                                                  cGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGAGCGACCAAACTGCTCGTTCAGCACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCgacactgtaacc                       >  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)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                    
GCGCCAGATCAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTGGGCGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGAGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAG  >  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: