Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 876,933 T→A 6.6% Q368L (CAG→CTG)  rimO ← ribosomal protein S12 methylthiotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913876,9330TA 6.6% 214.7 / 5.6 61Q368L (CAG→CTG) rimOribosomal protein S12 methylthiotransferase
Reads supporting (aligned to +/- strand):  ref base T (30/27);  new base A (2/2);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.45e-01

GTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGATCCAGACGCGCTTCTTTCAGGAAGTCGAGTA  >  NC_000913/876785‑877080
                                                                                                                                                    |                                                                                                                                                   
gTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATctgct                                                                                                                                                     >  3:364281/1‑149 (MQ=255)
 ttGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATctgctg                                                                                                                                                    <  4:392023/149‑1 (MQ=255)
   ggTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGATCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATctgctgct                                                                                                                                                  >  1:593257/1‑149 (MQ=255)
     tttCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCtgctgctgc                                                                                                                                                >  4:580315/1‑149 (MQ=255)
         aCCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTgctgcagct                                                                                                                                            >  2:232848/1‑149 (MQ=255)
            aTTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCAgctgcagcagca                                                                                                                                         >  4:674718/1‑149 (MQ=255)
            aTTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCAgctgcagcagca                                                                                                                                         <  3:157702/149‑1 (MQ=255)
            aTTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCAgctgcagcagca                                                                                                                                         >  4:157702/1‑149 (MQ=255)
            aTTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCAgctgcagcagca                                                                                                                                         <  3:674718/149‑1 (MQ=255)
                 gaTGAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCAt                                                                                                                                        >  4:588751/1‑145 (MQ=255)
                 gaTAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCAt                                                                                                                                        <  3:588751/145‑1 (MQ=255)
                 gaTAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGaaa                                                                                                                                    >  2:145578/1‑149 (MQ=255)
                        cgcgCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTcc                                                                                                                             <  2:9242/149‑1 (MQ=255)
                          cgcCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAg                                                                                                                           <  1:232848/149‑1 (MQ=255)
                           gcCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGc                                                                                                                          <  4:566828/149‑1 (MQ=255)
                              gTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGtt                                                                                                                               <  1:671658/141‑1 (MQ=255)
                              gTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGtt                                                                                                                               >  2:671658/1‑141 (MQ=255)
                                 gATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGttctt                                                                                                                    <  2:107614/149‑1 (MQ=255)
                                        ggTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAAc                                                                                                             <  1:252828/149‑1 (MQ=255)
                                        ggTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAAc                                                                                                             <  2:519387/149‑1 (MQ=255)
                                         gTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAAct                                                                                                            >  3:307372/1‑149 (MQ=255)
                                          tGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAActt                                                                                                           <  1:651630/149‑1 (MQ=255)
                                                tGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGaaa                                                                                                                                    <  4:141854/118‑1 (MQ=255)
                                                tGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGaaa                                                                                                                                    >  3:141854/1‑118 (MQ=255)
                                                   cATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGaa                                                                                                  >  3:131926/1‑149 (MQ=255)
                                                   cATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGaa                                                                                                  >  2:173097/1‑149 (MQ=255)
                                                         gcGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGt                                                                                            <  1:277292/149‑1 (MQ=255)
                                                           gACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACgg                                                                                                                                 >  2:374452/1‑110 (MQ=255)
                                                           gACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACgg                                                                                                                                 <  1:374452/110‑1 (MQ=255)
                                                           gACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTca                                                                                          >  1:524859/1‑149 (MQ=255)
                                                           gACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTca                                                                                          >  3:338883/1‑149 (MQ=255)
                                                            aCCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCgg                                                                                                    <  3:434207/138‑1 (MQ=255)
                                                            aCCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCgg                                                                                                    >  4:434207/1‑138 (MQ=255)
                                                              cAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTcaggc                                                                                       >  2:80566/1‑149 (MQ=255)
                                                                 tCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCGGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTcaggcagg                                                                                    <  3:6975/149‑1 (MQ=255)
                                                                 tCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTcaggcagg                                                                                    <  2:550812/149‑1 (MQ=255)
                                                                       cttcttCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTg                                                                              >  3:220626/1‑149 (MQ=255)
                                                                       cttcttCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTg                                                                              <  4:364281/149‑1 (MQ=255)
                                                                              gTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTg                                                                       <  3:113923/149‑1 (MQ=255)
                                                                                ccACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCa                                                                     >  4:247175/1‑149 (MQ=255)
                                                                                  aCTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTcaggcag                                                                                     <  1:128486/131‑1 (MQ=255)
                                                                                  aCTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTcaggcag                                                                                     >  2:128486/1‑131 (MQ=255)
                                                                                   cTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCt                                                                  >  2:493315/1‑149 (MQ=255)
                                                                                             aaTCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTcaggcagg                                                                                    <  2:235114/121‑1 (MQ=255)
                                                                                             aaTCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTcaggcagg                                                                                    >  1:235114/1‑121 (MQ=255)
                                                                                              aTCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCt                                                       >  3:276801/1‑149 (MQ=255)
                                                                                                 aCCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTa                                                    >  3:447367/1‑149 (MQ=255)
                                                                                                 aCCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTa                                                    >  2:613994/1‑149 (MQ=255)
                                                                                                  ccAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTAt                                                   >  3:375116/1‑149 (MQ=255)
                                                                                                  ccAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTAt                                                   <  4:375116/149‑1 (MQ=255)
                                                                                                   cAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTAtt                                                  >  1:517767/1‑149 (MQ=255)
                                                                                                   cAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTAtt                                                  >  1:621314/1‑149 (MQ=255)
                                                                                                        tttCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTc                                                                                                                              <  3:460894/68‑1 (MQ=255)
                                                                                                        tttCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTc                                                                                                                              >  4:460894/1‑68 (MQ=255)
                                                                                                                  cACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAAcgc                                   <  4:307372/149‑1 (MQ=255)
                                                                                                                     tttCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGAt                                >  1:495315/1‑149 (MQ=255)
                                                                                                                     tttCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGAt                                >  1:468368/1‑149 (MQ=255)
                                                                                                                       tctcTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGATcc                              >  3:255803/1‑149 (MQ=255)
                                                                                                                         tctTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTcaggcagg                                                                                    >  3:311727/1‑93 (MQ=255)
                                                                                                                         tctTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTcaggcagg                                                                                    <  4:311727/93‑1 (MQ=255)
                                                                                                                                     tCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCTGAACATGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGGAGCCAACGCGATCCAGACGCGCTTCttt                <  2:495315/149‑1 (MQ=255)
                                                                                                                                        gcggAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGATCCAGACGCGCTTCTTTCAg             <  2:267395/149‑1 (MQ=255)
                                                                                                                                          ggAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGATCCAGACGCGCTTCTTTCAGGa           <  4:338883/149‑1 (MQ=255)
                                                                                                                                            aaaTCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGATCCAGACGCGCTTCTTTCAGGAAg         >  2:390204/1‑149 (MQ=255)
                                                                                                                                              aTCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGATCCAGACGCGCTTCTTTCAGGAAGTc       <  1:173097/149‑1 (MQ=255)
                                                                                                                                                   ctgctgcAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGATCCAGACGCGCTTCTTTCAGGAAGTCGAGTa  <  4:447367/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGATCCAGACGCGCTTCTTTCAGGAAGTCGAGTA  >  NC_000913/876785‑877080

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