Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,745,0890AC30.6% 131.9 / 13.2 62V84G (GTG→GGG) trmDtRNA m(1)G37 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (38/5);  new base C (0/19);  total (38/24)
Fisher's exact test for biased strand distribution p-value = 9.93e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCG  >  NC_000913/2744868‑2745298
                                                                                                                                                                                                                                |                                                                                                                                                                                                                 
ccGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTgc                                                                                                                                                                                                            >  1:752017/1‑229 (MQ=255)
ccGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTgc                                                                                                                                                                                                            >  1:940450/1‑229 (MQ=255)
ccGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCAGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTgc                                                                                                                                                                                                            >  1:1780151/1‑229 (MQ=255)
 tgaaaTAAACCGGGAATCGAAGTCAATCAGCTTCATTGCTGTTACCTCACCACCTCTGAGAACGTAATCGCCTATTGACCTTTCTGCGTCAATTTCGGTTTGGATCACGCGCTCATCTATCCCTTCGTAGCGACCGCACCCCAGAATCATTTTTTGATTCTTTGCCAGTTCGCTGACGCCCCCTTGATCAAGCTTGCGTCCCTGGGGTGAAGGAAAAATCCCCTTTggg                                                                                                                                                                                                           <  2:1736837/226‑3 (MQ=255)
                         aaTCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATCCCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTtgc                                                                                                                                                                                    <  1:1939214/228‑1 (MQ=255)
                                 tcgttGCTGGTCACTCACCCACCGTGAGAACTCAATCGCGTATACACTCTTCTCCGTCATTTTCGGGTGGGACCACCCGCCAATCTACCCTTGCGTACCACCCGCCCCCAGGAATAATTTTTTTTGTTTTTGCCAGCTCGCGGACGCCCGTTGGATCAAGCTGGCGCCCTTGTGGTGACAAAAAAACCCCCTTGGCCCCTCCCCCCGCCGCGGCTTTTGCTGCATGaa                                                                                                                                                                            <  2:88629/225‑1 (MQ=255)
                                   aTTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATgg                                                                                                                                                                         >  2:14108/1‑229 (MQ=255)
                                     tGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCg                                                                                                                                                                       >  2:842626/1‑229 (MQ=255)
                                           tAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCttt                                                                                                                                                                                                              <  1:1445033/184‑1 (MQ=255)
                                           tAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCttt                                                                                                                                                                                                              >  2:1442558/1‑184 (MQ=255)
                                                 ccctccacTCAGAACGTAATCCCCGAGTTATCGTTCTTTGTCAATTTAGTCTTGGTTCACGCACACATCTATACCTTCTTAGCGACCGCACACCAGAAACAATTTTTGATTCGTTGCCAGTTCGCTGACCCCCGCCTGATCAAGCTTCCGTCCCTGTGGTGCCAGATAAATCATCTTTGCGCCTTCCCCCTCCTCGGCCTTTGCCACCTGAAGCTCCTCCCGCAAgc                                                                                                                                                             >  2:746616/5‑226 (MQ=255)
                                                           agaACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACcatc                                                                                                                                                  >  2:1817982/1‑228 (MQ=255)
                                                           agaACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACcatc                                                                                                                                                  >  2:1367172/1‑228 (MQ=255)
                                                                gTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCACGCAAGGGTTGCACCATCATTaa                                                                                                                                             >  2:1632717/1‑228 (MQ=255)
                                                                  aaTCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACa                                                                                                                                           >  1:1247692/1‑228 (MQ=255)
                                                                               cttgttcgcGTCAGTTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCAACCGCCCACCAAAATAATTTTTGGATGCGTTGCCGTTCCGCTGACCCCGGCTGGACCAAGCTGGCGCCCCTGGGGTGAAAAATAAATCACTTTGGGCCCTTCCCCCGCCGCGGCTTTTGCGGAATGAATGGCGTCCCGAAAGGGTTGCACCATAATTAACATCCCCGGTccgccg                                                                                                                             <  2:33405/221‑1 (MQ=255)
                                                                                        gTCCATTTCGGTTTGGATCACGCCCTCATCTATACCTTCGTAGCGACCTCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTTCGTCCCTGTGGTGACAGCTAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGCCGTCCCGCAAGGGCTGCACCAACATTGACATCCCCGGTCACCCGCCGTCAgg                                                                                                                     >  2:565997/1‑228 (MQ=255)
                                                                                         tCAATGGCGGGTTGGACCACGCGCCCGACTATTCCTTCGTTACGCCCCCACCCCAAAAACAAATTTTGATTCTGTTCCAATTCGCTTGCCCCCGGTTTGATAACGTTGGCGCCCCTGTGGAGACAAAAAAACCCCTTTTGCGCCTCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGa                                                                                                                    <  2:1777125/229‑1 (MQ=255)
                                                                                          cAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGTATGGCGTCCCGCAAGGGCTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACg                                                                                                                  >  2:231582/1‑229 (MQ=255)
                                                                                                       agtCACGCTCTACTCTATCTCTTCTGAGCGACACCACACAAGAATAATTTTTTGCGCCTTTTCCATTTCGTCGCCGCCCGCTTGATCAAGTTTGGGCTCGTGTGCTGACAGAAAAACCCCCTTGGCCCCTCCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTACCATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCacggta                                                                                                     <  2:1642614/227‑1 (MQ=255)
                                                                                                        ttcctgggttCAACTATCCCTTCGAGCCGACCCCACGCCAAAATCATTTTTTGTGTGGGTCCCAGTCCGCTGCCCCCCGTTTGTCCATGCTGCCTTCCGTGGGGTGCAAGATAAACCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCacggtac                                                                                                    <  2:499858/220‑1 (MQ=255)
                                                                                                                 tCATCTATCCCTTCGTAGCGCCCGCACCCAAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGACCAAGCTTGCGCCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCa                                                                                                                                                       <  2:1391982/169‑1 (MQ=255)
                                                                                                                 tCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCa                                                                                                                                                       >  1:1394378/1‑169 (MQ=255)
                                                                                                                  cATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCCTCCTTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGCGCCGGTc                                                                                          >  1:1963500/1‑229 (MQ=255)
                                                                                                                             tcgacgcgcCCGCAGCCCAGAATAATTTTTGATTTCTTTCCCAGTTCCCTGCCGCCCGTTTGACCAAGCTGGCGTCCTGTGGGTGACAAATAACTCCCTTTTGCCCCTCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGCCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATc                                                                                                               <  2:1244276/192‑1 (MQ=255)
                                                                                                                             tCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACccgcc                                                                                                                                                                                               >  2:340775/1‑117 (MQ=255)
                                                                                                                             tCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACccgcc                                                                                                                                                                                               <  1:341391/117‑1 (MQ=255)
                                                                                                                             tCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCcgccgc                                                                                                                                                                                             >  2:1602795/1‑119 (MQ=255)
                                                                                                                             tCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCcgccgc                                                                                                                                                                                             <  1:1605527/119‑1 (MQ=255)
                                                                                                                             tCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATc                                                                                                               >  1:1246444/1‑197 (MQ=255)
                                                                                                                             tCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATCCGTGAAg                                                                                >  2:689865/1‑228 (MQ=255)
                                                                                                                             tCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGCAAGGACGATCGTCCACGGTCCGGTGCCGGTCATCCGTGAAg                                                                                >  2:1849034/1‑228 (MQ=255)
                                                                                                                             tCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATACCCGGTCCGCCGCCGTCAGGACGATCGTCCACGGTACGGTGCCGGTcgcccttgaag                                                                                >  2:40169/1‑218 (MQ=255)
                                                                                                                              cGTGGCGCCCGCACACAAGAATCAATTTTTGTTTCGTTCCCGTTCCGCTGCCCCCCGCTGGATCAAGCTTGCGCCCGTGTGGTGAAAAAAAAACCCCCTTTGCGCCTCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTccgc                                                                                                                               <  2:860356/180‑1 (MQ=255)
                                                                                                                              cGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTccgc                                                                                                                               >  1:861859/1‑180 (MQ=255)
                                                                                                                                tAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGGTCGTCCACGGTACGGTGCCGGTCATGCGttcagtagc                                                                            >  2:1816662/1‑221 (MQ=255)
                                                                                                                                    gACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGAGCGTCCACGGGACGGTGCCGGTCCTGCGTGAAGTCGCGAgg                                                                        >  1:1742227/1‑229 (MQ=255)
                                                                                                                                      ccGCACACCAGAATAAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGAAAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCa                                                                                                                                                       >  1:503617/1‑148 (MQ=255)
                                                                                                                                      ccGCACACCAGAACAAATTTTGCATTCTTTCCCATTCCGCTGACGCCCGCTTGTTCAAGCTTGTGCCCCGGTGTTCCCAGATAATCCCCCTTCGCCCTTCCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCa                                                                                                                                                       <  2:502733/148‑1 (MQ=255)
                                                                                                                                       cGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTACCATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACt                                                                     >  1:1079581/1‑229 (MQ=255)
                                                                                                                                       cGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTACCATCCCCGGTCCGCCGCCGTAAGGACGATCGCCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACt                                                                     >  1:320614/1‑229 (MQ=255)
                                                                                                                                                        ttttGATTCGTTGCCAGTTCGCTGCCCCCCTCTTGATCAAGCTGGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCa                                                                                                                                                                <  2:1511249/121‑1 (MQ=255)
                                                                                                                                                        ttttGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCa                                                                                                                                                                >  1:1513866/1‑121 (MQ=255)
                                                                                                                                                                              cTGCCGCCCGCTTGATCAAGCTTGCGTCCTTGTGGTGACAGATAAATCCCCTTGGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCa                                                                                                                                                       <  2:1852638/109‑1 (MQ=255)
                                                                                                                                                                              cTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCa                                                                                                                                                       >  1:1855795/1‑109 (MQ=255)
                                                                                                                                                                                         cTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGCCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTGACTGCCCGGCCAGTTAc                     >  2:1684298/1‑229 (MQ=255)
                                                                                                                                                                                                  gtttGCGCCCCTGGGGGGACAGAAAAACCCCTTTGGCGCCTCCCCCCGCCGCGGCTTTTGCTGCATGAAGGGCGTCCCGAAAGGGTTGCCCCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATc                                                                                                               <  2:1388326/128‑1 (MQ=255)
                                                                                                                                                                                                  gCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATc                                                                                                               >  1:1390716/1‑130 (MQ=255)
                                                                                                                                                                                                       cGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAgg                                               >  1:758436/1‑189 (MQ=255)
                                                                                                                                                                                                            cTGTGGTGACAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATgg                                                                                                                                                                         <  1:1427998/62‑1 (MQ=255)
                                                                                                                                                                                                            cTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATgg                                                                                                                                                                         >  2:1425550/1‑62 (MQ=255)
                                                                                                                                                                                                            cTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCAt                                                                                        >  1:622863/1‑143 (MQ=255)
                                                                                                                                                                                                            cTGGGGGGACAGATAAACCCCTTTTGCGCCTCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCAt                                                                                        <  2:621770/143‑1 (MQ=255)
                                                                                                                                                                                                             tgtgGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACcatc                                                                                                                                                  >  2:473373/1‑84 (MQ=255)
                                                                                                                                                                                                             tgtgGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACcatc                                                                                                                                                  <  1:474209/84‑1 (MQ=255)
                                                                                                                                                                                                              gtggtcACATATATATCCCGTTCGCCCCTTCCCCCGCCGCGGCTTTTGCGGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTgcg  <  2:1739262/227‑1 (MQ=255)
                                                                                                                                                                                                                     aCAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTccgccgcc                                                                                                                           <  2:921141/100‑1 (MQ=255)
                                                                                                                                                                                                                     aCAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTccgccgcc                                                                                                                           >  1:922749/1‑100 (MQ=255)
                                                                                                                                                                                                                     aCAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAg                                                                         >  1:1537274/1‑150 (MQ=255)
                                                                                                                                                                                                                     aCAGAAAAACCCCTTTTGCGCTTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAg                                                                         <  2:1534621/150‑1 (MQ=255)
                                                                                                                                                                                                                         aTAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACcatc                                                                                                                                                  <  2:559903/73‑1 (MQ=255)
                                                                                                                                                                                                                         aTAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACcatc                                                                                                                                                  >  1:560887/1‑73 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                 
CCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCG  >  NC_000913/2744868‑2745298

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: