Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,913,956 T→C 5.1% K349E (AAG→GAG)  rlmD ← 23S rRNA m(5)U1939 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,913,9560TC 5.1% 266.8 / 3.7 78K349E (AAG→GAG) rlmD23S rRNA m(5)U1939 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (27/47);  new base C (2/2);  total (29/49)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

GCTATCCCGAGCCAGCGTTGCAGGGTTACAGGATACATAAACTATACGAATAGGTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCTTTTTCCACCAGCGCCGGAACACCTTCTACACCGACCACACTGGCAGCTTGTGTCGCCAATGGCAGT  >  NC_000913/2913807‑2914098
                                                                                                                                                     |                                                                                                                                              
gCTATCCCGAGCCAGCGTTGCAGGGTTACAGGATACATAAACTATACGAATAGGTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCtt                                                                                                                                                >  1:648988/1‑150 (MQ=255)
   aTCCCGAGCCAGCGTTGCAGGGTTACAGGATACATAACCTAGCCGAAGAGGTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGATTTTCGCCCACGGCTGCTTtgt                                                                                                                                             <  1:840536/150‑1 (MQ=255)
       cGAGCCAGCGTTGCAGGGTTACAGGATACATAAACTATACGAATAGGTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACa                                                                                                                                         >  1:916053/1‑150 (MQ=255)
       cGAGACAGCGTTGCAGTTTTTCTGGATACATAAACTATACGATTAGGTTCCAGTTTTATTATTTGCTTCATAACACCTGCGTCACGTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCTTTTTTCGCGGACGGCTGCTTTGTGACa                                                                                                                                         <  2:637063/150‑1 (MQ=255)
           ccAGCGTTGCAGGGTTACAGGATACATAAACTATACGAATAGGTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACAtctt                                                                                                                                     >  1:1377225/1‑150 (MQ=255)
           ccAGCGTTGCAGGGTTACAGGATACATAAACTATACGAATAGGTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACAtctt                                                                                                                                     >  1:565658/1‑150 (MQ=255)
           ccAGCGTTGCAGGGTTACAGGATACATAAACTATACGAATAGGTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTCTGTGACAtctt                                                                                                                                     >  1:565818/1‑150 (MQ=255)
           ccAGCGTTGCAGGGTTACAGGATACATAAACTATACGAATAGGTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTCTGTGACAtctt                                                                                                                                     >  1:565702/1‑150 (MQ=255)
                gTTGCAGGGTTACAGGATACATAAACTATACGAATAGGTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTGTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCa                                                                                                                                >  1:745880/1‑150 (MQ=255)
                    cAGGGTTACAGGATACATAAACTATACGAATAGGTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCGCACGGCTGCTTTGTGTCATCTTCTTCAAGAt                                                                                                                            <  1:1074516/150‑1 (MQ=255)
                     aGGGTTACAGGATACATAAACTATACGAATAGGTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGAtt                                                                                                                           >  2:237106/1‑150 (MQ=255)
                      gggTTACAGGATACATAAACTATACGAATAGGTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACGTCTTCTTCAAGAttt                                                                                                                          >  1:1224429/1‑150 (MQ=255)
                      gggTTACAGGATACATAAACTATACGAATAGGTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGAttt                                                                                                                          <  1:761930/150‑1 (MQ=255)
                       ggTTACAGGATACATAAACTATACGAATAGGTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGAtttt                                                                                                                         <  1:659966/150‑1 (MQ=255)
                        gTTACAGGATACATAAACTATACGAATAGGTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTc                                                                                                                        <  2:1176750/150‑1 (MQ=255)
                         ttACAGGATACATAAACTATACGAATAGGTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCg                                                                                                                       <  1:1093604/150‑1 (MQ=255)
                          tACAGGATACATAAACTATACGAATAGGTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGt                                                                                                                      <  2:690751/150‑1 (MQ=255)
                          tACAGGATACATAAACTATACGAATAGGTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGt                                                                                                                      <  2:1238342/150‑1 (MQ=255)
                                    aTAAACTATACGAATAGGTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGt                                                                                                            <  1:3696/150‑1 (MQ=255)
                                    aTAAACTATACGAATAGGTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGt                                                                                                            <  1:825237/150‑1 (MQ=255)
                                     tAAACTATACGAATAGGTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGAGACATCTTCTACAAGATTTTCGTGATAGAACGTc                                                                                                           >  1:14207/1‑150 (MQ=255)
                                      aaaCTCTACGAATAGGTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTca                                                                                                          <  2:1365157/150‑1 (MQ=255)
                                      aaaCTATACGAATAGGTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTca                                                                                                          <  2:353918/150‑1 (MQ=255)
                                      aaaCTATACGAATAGGTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTca                                                                                                          <  1:348640/150‑1 (MQ=255)
                                       aaCTATACGAATAGGTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTcac                                                                                                         >  2:1227779/1‑150 (MQ=255)
                                        aCTATACGACTAGGTTCCTGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGTTTTTCGTGATAAAACGTcaca                                                                                                        <  2:473007/150‑1 (MQ=255)
                                          tataCGAATAGGTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACAtt                                                                                                      <  1:293312/150‑1 (MQ=255)
                                          tataCGAATAGGTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACAtt                                                                                                      <  1:71305/150‑1 (MQ=255)
                                            taCGAATAGGTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCt                                                                                                    >  1:454004/1‑150 (MQ=255)
                                             aCGAATAGGTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTg                                                                                                   <  2:1039384/150‑1 (MQ=255)
                                                aaTAGGTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTaa                                                                                                >  2:713638/1‑150 (MQ=255)
                                                  tAGGTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGc                                                                                              <  2:35021/150‑1 (MQ=255)
                                                    ggTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCa                                                                                            <  1:961123/150‑1 (MQ=255)
                                                      ttCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCAtt                                                                                          >  1:51021/1‑150 (MQ=255)
                                                          aGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGa                                                                                      >  1:562683/1‑150 (MQ=255)
                                                          aGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGa                                                                                      >  1:299964/1‑150 (MQ=255)
                                                          aGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGa                                                                                      <  1:1198171/150‑1 (MQ=255)
                                                            ttttATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGAcg                                                                                    >  2:1269871/1‑150 (MQ=255)
                                                               tataATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACgcgc                                                                                 <  2:459515/150‑1 (MQ=255)
                                                                       gctgcATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTctgct                                                                         <  2:799180/150‑1 (MQ=255)
                                                                       gctgcATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTctgct                                                                         >  2:618206/1‑150 (MQ=255)
                                                                              aaCACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCttt                                                                  >  1:327418/1‑150 (MQ=255)
                                                                               acacCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCtttt                                                                 <  2:623223/150‑1 (MQ=255)
                                                                                  ccTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGTTTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCTTTTTcc                                                              <  1:1300807/150‑1 (MQ=255)
                                                                                     gCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCTTTTTccacc                                                           >  1:249931/1‑150 (MQ=255)
                                                                                            cTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCTTTTTCCACCAGCGCCg                                                    <  2:216519/150‑1 (MQ=255)
                                                                                                  gcgcCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCTTTTTCCACCAGCGCCGGAacac                                              <  2:640040/150‑1 (MQ=255)
                                                                                                  gcgcCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCTTTTTCCACCAGCGCCGGAacac                                              <  2:639880/150‑1 (MQ=255)
                                                                                                  gcgcCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCTTTTTCCACCAGCGCCGGAacac                                              >  2:245302/1‑150 (MQ=255)
                                                                                                     ccGGGTCCAGCTACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGTTTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCTTTTTCCAGCTGCGCCGGAACACCtt                                           <  2:832851/150‑1 (MQ=255)
                                                                                                      cGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCTTTTTCCACCAGCGCCGGAACACCTTc                                          <  2:160142/150‑1 (MQ=255)
                                                                                                      cGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCTTTTTCCACCAGCGCCGGAACACCTTc                                          <  2:160116/150‑1 (MQ=255)
                                                                                                       gggTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCTTTTTCCACCAGCGCCGGAACACCTTCt                                         >  2:326897/1‑150 (MQ=255)
                                                                                                         gTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCTTTTTCCACCAGCGCCGGAACACCTTCTac                                       <  2:634417/150‑1 (MQ=255)
                                                                                                             aGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCTTTTTCCACCAGCGCCGGAACACCTTCTACACCg                                   >  1:513768/1‑150 (MQ=255)
                                                                                                             aGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCTTTTTCCACCAGCGCCGGAACACCTTCTACACCg                                   >  1:83127/1‑150 (MQ=255)
                                                                                                             aGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCTTTTTCCACCAGCGCCGGAACACCTTCTACACCg                                   >  1:83120/1‑150 (MQ=255)
                                                                                                               cAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACGTCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCTTTTTCCACCAGCGCCGGAACACCTTCTACACCGAc                                 <  2:1224429/150‑1 (MQ=255)
                                                                                                                aaCACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCTTTTTCCACCAGCGCCGGAACACCTTCTACACCGAcc                                <  2:975137/150‑1 (MQ=255)
                                                                                                                     tttATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAAAGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCTTTTTCCACCAGCGCCGGAACACCATCTACACCGACCACACt                           >  1:138065/1‑150 (MQ=255)
                                                                                                                     tttATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATATTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCTTTTTCCACCAGCGCCGGAACACCTTCTACACCGACCacaca                           >  1:138060/1‑149 (MQ=255)
                                                                                                                       tATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCTTTTTCCACCAGCGCCGGAACACCTTCTACACCGACCACACTgg                         >  2:1403452/1‑150 (MQ=255)
                                                                                                                           gAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCTTTTTCCACCAGCGCCGGAACACCTTCTACACCGACCACACTGGCAGc                     <  2:565658/150‑1 (MQ=255)
                                                                                                                           gAAGCCGTTTTTCGCCCACGGCTGCTCTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCTTTTTCCACCAGCGCCGGAACACCTTCTACACCGACCACACTGGCAGc                     <  2:565702/150‑1 (MQ=255)
                                                                                                                           gAAGCCGTTTTTCGCCCACGGCTGCTCTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCTTTTTCCACCAGCGCCGGAACACCTTCTACACCGACCACACTGGCAGc                     <  2:565818/150‑1 (MQ=255)
                                                                                                                            aaGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCTTTTTCCACCAGCGCCGGAACACCTTCTACACCGACCACACTGGCAGCt                    <  2:691661/150‑1 (MQ=255)
                                                                                                                            aaGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCTTTTTCCACCAGCGCCGGAACACCTTCTACACCGACCACACTGGCAGCt                    <  2:691549/150‑1 (MQ=255)
                                                                                                                             aGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCTTTTTCCACCAGCGCCGGAACACCTTCTACACCGACCACACTGGCAGCtt                   <  2:462032/150‑1 (MQ=255)
                                                                                                                             aGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCAGTCTGCTGGCCTTTTTCCACCAGCGCCGGAACACCTTCTACACCGACCACACTGGCAGcgg                   >  2:1359345/1‑148 (MQ=255)
                                                                                                                                 gTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTTATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCTTTTTCCACCAGCGCCGGAACACCTTCTACACCGACCACACTGGCAGCTTgtgt               <  2:270719/150‑1 (MQ=255)
                                                                                                                                 gTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCTTTTTCCACCAGCGCCGGAACACCTTCTACACCGACCACACTGGCAGCTTgtgt               <  2:1192744/150‑1 (MQ=255)
                                                                                                                                 gTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCTTTGTCCACCAGCGCCGGAACACCTTCTACACCGACCACCCTGGCAGCTTgtgt               <  2:745880/150‑1 (MQ=255)
                                                                                                                                   ttttCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCTTTTTCCACCAGCGCCGGAACACCTTCTACACCGACCACACTGGCAGCTTGTGTCg             <  1:396342/150‑1 (MQ=255)
                                                                                                                                       cGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCTTTTTCCACCAGCGCCGGAACACCTTCTACACCGACCACACTGGCAGCTTGTGTCGCCaa         <  2:1094666/150‑1 (MQ=255)
                                                                                                                                           cACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCTTTTTCCACCAGCGCCGGAACACCTTCTACACCGACCACACTGGCAGCTTGTGTCGCCAATGGc     <  1:166953/150‑1 (MQ=255)
                                                                                                                                           cACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCTTTTTCCACCAGCGCCGGAACACCTTCTACACCGACCACACTGGCAGCTTGTGTCGCCAATGGc     <  1:166941/150‑1 (MQ=255)
                                                                                                                                           cACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCCTTTTCCACCAGCGCCGGAACACCTTCTACACCGACCACACTGGCAGCTTGTGTCGCCAATGGc     >  1:782366/1‑150 (MQ=255)
                                                                                                                                            aCGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCTTTTTCCACCAGCGCCGGAACACCTTCTACACCGACCACACTGGCAGCTTGTGTCGCCAATGGCa    <  2:703852/150‑1 (MQ=255)
                                                                                                                                              ggCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCTTTTTCCACCAGCGCCGGAACACCTTCTACACCGACCACACTGGCAGCTTGTGTCGCCAATGGCAgt  <  2:115682/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
GCTATCCCGAGCCAGCGTTGCAGGGTTACAGGATACATAAACTATACGAATAGGTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACACTTTATCGAAGCCGTTTTTCGCCCACGGCTGCTTTGTGACATCTTCTTCAAGATTTTCGTGATAAAACGTCACATTCTGTAAGCCATTAAGACGCGCATTCTGCTGGCCTTTTTCCACCAGCGCCGGAACACCTTCTACACCGACCACACTGGCAGCTTGTGTCGCCAATGGCAGT  >  NC_000913/2913807‑2914098

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