Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,643,989 C→T 6.4% V99I (GTT→ATT) ‡ rlmN ← 23S rRNA m(2)A2503 methyltransferase/tRNA m(2)A37 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,643,9890CT 6.4% 447.0 / 11.6 141V99I (GTT→ATT) ‡rlmN23S rRNA m(2)A2503 methyltransferase/tRNA m(2)A37 methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (64/68);  new base T (7/2);  total (71/70)
Fisher's exact test for biased strand distribution p-value = 1.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 142 total aligned reads displayed.

AATTTCCGACACCCGCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATAC  >  NC_000913/2643840‑2644137
                                                                                                                                                     |                                                                                                                                                    
  tttCCGACACCCGCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACgg                                                                                                                                                    >  1:196134/1‑150 (MQ=255)
  tttCCGACACCCGCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACgg                                                                                                                                                    <  1:2329569/150‑1 (MQ=255)
     ccGACACCCGCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGa                                                                                                                                                 >  2:1244636/1‑150 (MQ=255)
       gACACCCGCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATc                                                                                                                                               <  2:852918/150‑1 (MQ=255)
         cacCCGCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGc                                                                                                                                             <  1:2379312/150‑1 (MQ=255)
          acCCGCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCGTTGGCGATcagc                                                                                                                                            <  1:492162/150‑4 (MQ=255)
           cccGCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCGTTGGCGATcagcg                                                                                                                                           <  2:1462025/150‑5 (MQ=255)
            ccGCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCa                                                                                                                                          <  1:1096254/150‑1 (MQ=255)
              gCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCAtt                                                                                                                                        <  2:666299/150‑1 (MQ=255)
               cAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGGCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCAttt                                                                                                                                       <  2:809830/150‑1 (MQ=255)
                   ttGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATg                                                                                                                                   <  1:1404438/150‑1 (MQ=255)
                   ttGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATg                                                                                                                                   >  1:2333107/1‑150 (MQ=255)
                    tGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATgg                                                                                                                                  <  2:939477/150‑1 (MQ=255)
                     gCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGt                                                                                                                                 >  1:1285000/1‑150 (MQ=255)
                     gCGGTTAAAGCCCTGATGGGCGTTGGAACAGAATTTACACTCCTGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGt                                                                                                                                 <  2:974965/150‑1 (MQ=255)
                      cGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGGCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTg                                                                                                                                <  1:1020760/150‑1 (MQ=255)
                       ggTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGc                                                                                                                               >  2:749519/1‑150 (MQ=255)
                                ccTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATg                                                                                                                      <  2:22381/150‑1 (MQ=255)
                                      gggCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCt                                                                                                                <  1:466105/150‑1 (MQ=255)
                                        gCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGt                                                                                                              <  2:2056395/150‑1 (MQ=255)
                                          ggTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGttc                                                                                                            <  1:839494/150‑1 (MQ=255)
                                          ggTGGAACAGAATTTACACTCCAGAGAACACACCACCTGCGAAGAGACGAAGAGCGTGGCAAGGTCGTATTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGtta                                                                                                            >  2:1078374/1‑149 (MQ=255)
                                           gTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGttct                                                                                                           >  1:1423084/1‑150 (MQ=255)
                                           gTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGttct                                                                                                           <  1:1158360/150‑1 (MQ=255)
                                           gTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGttct                                                                                                           >  2:216942/1‑150 (MQ=255)
                                              gAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCa                                                                                                        <  2:766108/150‑1 (MQ=255)
                                               aaCAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCaa                                                                                                       >  2:1285821/1‑150 (MQ=255)
                                                aCAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAAc                                                                                                      <  2:2257054/150‑1 (MQ=255)
                                                 cAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAAcc                                                                                                     >  1:25362/1‑150 (MQ=255)
                                                 cAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAAcc                                                                                                     <  1:230117/150‑1 (MQ=255)
                                                    aaTTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACt                                                                                                  >  1:1044321/1‑150 (MQ=255)
                                                     aTTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACtt                                                                                                 >  2:851756/1‑150 (MQ=255)
                                                      tttACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTc                                                                                                >  1:910730/1‑150 (MQ=255)
                                                       ttACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTcc                                                                                               <  2:2057699/150‑1 (MQ=255)
                                                        tACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCg                                                                                              >  1:2108840/1‑150 (MQ=255)
                                                        tACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCg                                                                                              >  2:129021/1‑150 (MQ=255)
                                                         acacTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCgg                                                                                             <  2:1875781/150‑1 (MQ=255)
                                                             tCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTgcg                                                                                         <  1:154351/150‑1 (MQ=255)
                                                               cAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGcgcg                                                                                       <  1:1388025/150‑1 (MQ=255)
                                                                aGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCgg                                                                                      <  2:1161524/150‑1 (MQ=255)
                                                                 gcgcACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGa                                                                                     <  1:2022152/150‑1 (MQ=255)
                                                                  cccacaCCCCACCTGAGAACAGACGCAGAGCGTGGAACGGACGTCTTCCGGGATATACACCGTTTCGACGCGCTCATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGAt                                                                                    >  2:1541508/3‑150 (MQ=255)
                                                                    cacaCCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGAttt                                                                                  >  1:2465200/1‑150 (MQ=255)
                                                                    cacaCCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGAttt                                                                                  >  1:1863540/1‑150 (MQ=255)
                                                                        ccccACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGGTCTTCAACCACTTCCGGTGCGCGGATTTCCGc                                                                              >  1:264737/1‑150 (MQ=255)
                                                                         cccACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGcc                                                                             >  2:1337011/1‑150 (MQ=255)
                                                                         cccACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGcc                                                                             <  2:254012/150‑1 (MQ=255)
                                                                         cccACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGcc                                                                             <  2:351610/150‑1 (MQ=255)
                                                                           cACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCAc                                                                           <  2:218950/150‑1 (MQ=255)
                                                                             ccTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACct                                                                         >  2:1006063/1‑150 (MQ=255)
                                                                                 cGAAGAGACGCAGATCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCttt                                                                     >  2:2305099/1‑150 (MQ=255)
                                                                                  gAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTc                                                                    <  2:1492776/150‑1 (MQ=255)
                                                                                  gAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTc                                                                    >  1:1671521/1‑150 (MQ=255)
                                                                                     gagaCGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCGCTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGt                                                                 >  2:1240799/1‑150 (MQ=255)
                                                                                     gagaCGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGt                                                                 >  2:2321769/1‑150 (MQ=255)
                                                                                       gaCGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTGATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGttt                                                               >  2:46639/1‑150 (MQ=255)
                                                                                       gaCGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGttt                                                               >  2:1417094/1‑150 (MQ=255)
                                                                                       gaCGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGttt                                                               >  2:1736237/1‑150 (MQ=255)
                                                                                         cGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGc                                                             >  1:177692/1‑150 (MQ=255)
                                                                                                      aCGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACttt                                                >  1:2186802/1‑150 (MQ=255)
                                                                                                         gtcgtcTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACTCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTttgtt                                             <  1:70631/150‑1 (MQ=255)
                                                                                                         gtcgtcTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTttgtt                                             >  2:1604152/1‑150 (MQ=255)
                                                                                                         gtcgtcTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTttgtt                                             <  1:79629/150‑1 (MQ=255)
                                                                                                         gtcgtcTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTttgtt                                             <  2:2465200/150‑1 (MQ=255)
                                                                                                         gtcgtcTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTttgtt                                             <  1:1037079/150‑1 (MQ=255)
                                                                                                          tcgtcTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTttgttg                                            <  2:196134/150‑1 (MQ=255)
                                                                                                          tcgtcTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTttgttg                                            <  1:1873931/150‑1 (MQ=255)
                                                                                                           cgtgTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGa                                           <  2:1259514/150‑1 (MQ=255)
                                                                                                           cgtcTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGa                                           >  1:1023514/1‑150 (MQ=255)
                                                                                                           cgtcTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGa                                           <  1:2128983/150‑1 (MQ=255)
                                                                                                              cTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGt                                        >  1:1513824/1‑150 (MQ=255)
                                                                                                               ttCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTc                                       <  2:1423084/150‑1 (MQ=255)
                                                                                                               ttCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTc                                       >  2:39756/1‑150 (MQ=255)
                                                                                                               ttCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTc                                       <  2:2342138/150‑1 (MQ=255)
                                                                                                                 ccGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCgg                                     <  2:2333551/150‑1 (MQ=255)
                                                                                                                  cGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGt                                    >  1:1307725/1‑150 (MQ=255)
                                                                                                                  cGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGt                                    >  1:1419810/1‑150 (MQ=255)
                                                                                                                   gggATATAAACAGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTAAGATGAACGCTGTTCTTCCACCACTTCCGGTGCGCGGATTTCCGCCCCCGCTTTCAGTTTGCCGCGCAACACTTGGTTGATGTCGGtc                                   >  2:469557/1‑150 (MQ=255)
                                                                                                                      atatACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGtcatc                                <  1:394024/150‑1 (MQ=255)
                                                                                                                      atatACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGtcatc                                <  1:848675/150‑1 (MQ=255)
                                                                                                                        ataCACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGtcatctc                              <  1:1746203/150‑1 (MQ=255)
                                                                                                                         taCACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGtcatctca                             <  2:177692/150‑1 (MQ=255)
                                                                                                                           cacCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGtcatctcatc                           >  1:1065217/1‑150 (MQ=255)
                                                                                                                               gTTTCGACGCGCTGATCGGGAACGGCGATCGCCCATTTAATGGTGCAGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGGGGATTTCGGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTGGGTCATCTCATCAAAg                       >  2:1467388/1‑150 (MQ=255)
                                                                                                                               gTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAg                       <  2:1303037/150‑1 (MQ=255)
                                                                                                                                   cGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGt                   >  2:2111217/1‑150 (MQ=255)
                                                                                                                                     aatgggcgATCGCCGTTGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCg                 >  1:646899/8‑150 (MQ=255)
                                                                                                                                       tgggcgATCGCCGTTGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCgca               >  1:1056145/6‑150 (MQ=255)
                                                                                                                                        cgcTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCGCTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCgcag              <  1:1240799/150‑1 (MQ=255)
                                                                                                                                        cgcTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCgcag              >  2:640062/1‑150 (MQ=255)
                                                                                                                                        cgcTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCgcag              >  1:1648353/1‑150 (MQ=255)
                                                                                                                                         ggcgATCGCCGTTGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCgcagc             >  2:1093108/4‑150 (MQ=255)
                                                                                                                                         ggcgATCGCCGTTGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCgcagc             >  1:2174844/4‑150 (MQ=255)
                                                                                                                                         ggcgATCGCCGTTGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCgcagc             >  2:377010/4‑150 (MQ=255)
                                                                                                                                          gcgATCGCCGTTGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCgcagca            >  1:754374/3‑150 (MQ=255)
                                                                                                                                            gATCGCCGTTGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAAt          >  1:19310/1‑150 (MQ=255)
                                                                                                                                            gATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAAt          >  2:1714510/1‑150 (MQ=255)
                                                                                                                                             aTCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATa         >  2:2418291/1‑150 (MQ=255)
                                                                                                                                                   aaCGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATa   <  1:639377/150‑1 (MQ=255)
                                                                                                                                                    aCGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATAc  >  1:1581214/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AATTTCCGACACCCGCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATAC  >  NC_000913/2643840‑2644137

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