Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,153,464 T→G 6.2% intergenic (+10/‑182) sibA → / → sibB small RNA SibA/small RNA SibB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,153,4640TG 6.2% 247.6 / 7.6 81intergenic (+10/‑182)sibA/sibBsmall RNA SibA/small RNA SibB
Reads supporting (aligned to +/- strand):  ref base T (34/42);  new base G (2/3);  total (36/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

GGGAGAGGTTTCCCCCTCCCCCTGGTGTTCTTAGTAAGCCTGGAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGATTATTGCGTAACGACTCGCAAAGAATATTTGCTGGCAGGATCGCAGACTACAAAGCCTGCGGATTG  >  NC_000913/2153319‑2153612
                                                                                                                                                 |                                                                                                                                                    
gggAGAGGTTTCCCCCTCCCCCTGGTGTTCTTAGTAAGCCTGGAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCg                                                                                                                                                  <  2:477979/150‑1 (MQ=255)
gggAGAGGTTTCCCCCTCCCCCTGGTGTTCTTAGTAAGCCTGGAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCg                                                                                                                                                  <  2:473209/150‑1 (MQ=255)
gggAGAGGTTTCCCCCTCCCCCTGGTGTTCTTAGTAAGCCTGGAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCg                                                                                                                                                  >  2:440804/1‑150 (MQ=255)
 ggAGAGGTTTCCCCCTCCCCCTGGTGTTCTTAGTAAGCCTGGAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGc                                                                                                                                                 >  2:791526/1‑150 (MQ=255)
 ggAGAGGTTTCCCCCTCCCCCTGGTGTTCTTAGTAAGCCTGGAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGc                                                                                                                                                 >  2:546105/1‑150 (MQ=255)
  gagaGGTTTCCCCCTCCCCCTGGTGTTCTTAGTAAGCCTGGAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCt                                                                                                                                                >  2:660506/1‑150 (MQ=255)
      ggTTTCCCCCTCCCCCTGGTGTTCTTAGTAAGCCTGGAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGaaa                                                                                                                                            >  1:1216075/1‑150 (MQ=255)
        tttCCCCCTCCCCCTGGTGTTCTTAGTAAGCCTGGAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGGGACGCTGaaacg                                                                                                                                          <  2:870775/150‑3 (MQ=255)
         ttCCCCCTCCCCCTGGTGTTCTTAGTAAGCCTGGAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAatca                                                                                                                                         <  2:1267039/150‑1 (MQ=255)
         ttCCCCCTCCCCCTGGTGTTCTTAGTAAGCCTGGAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGGGACGCTGaaacgg                                                                                                                                         <  1:1020422/150‑4 (MQ=255)
          tccccctcccccTGGTGTTCTTAGTAAGCCTGGAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTACTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAatcat                                                                                                                                        >  2:951666/1‑150 (MQ=255)
            cccctccccctGGTGTTCTTAGTAAGCCTGGAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCg                                                                                                                                      <  1:1347997/150‑1 (MQ=255)
             ccctccccctGGTGTTCTTAGTAAGCCTGGAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGGGACGCTGaaacgggaaa                                                                                                                                     <  2:604392/150‑8 (MQ=37)
              cctccccctGGTGTTCTTAGTAAGCCTGGAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGc                                                                                                                                    <  1:991357/150‑1 (MQ=255)
              cctccccctGGTGTTCTTAGTAAGCCTGGAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGc                                                                                                                                    <  1:168326/150‑1 (MQ=255)
                tccccctGGTGTTCTTAGTAAGCCTGGAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCtt                                                                                                                                  <  2:1100028/150‑1 (MQ=255)
                   ccctGGTGTTCTTAGTAAGCCTGGAAGCTAATCCCTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTAcc                                                                                                                               <  1:1000269/150‑1 (MQ=255)
                   acatGGTTTTCTTAGTAAGCCTTGAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTAcc                                                                                                                               <  1:1100955/147‑1 (MQ=255)
                      tGGTGTTCTTAGTAAGCCTGGAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTAcctcc                                                                                                                            <  2:1216075/150‑1 (MQ=255)
                       ggTGTTCTTAGTAAGCCTGGAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTAcctcct                                                                                                                           <  2:1014095/150‑1 (MQ=255)
                               tAGTAAGCCTGGAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCATCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCt                                                                                                                   >  2:510669/1‑150 (MQ=255)
                                aGTAAGCCTGGAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTg                                                                                                                  >  2:474000/1‑150 (MQ=255)
                                aGTAAGCCTGGAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTg                                                                                                                  <  1:951666/150‑1 (MQ=255)
                                  tAAGCCTGGAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGaa                                                                                                                <  1:195978/150‑1 (MQ=255)
                                   aaGCCTGGAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAAt                                                                                                               >  2:1409603/1‑150 (MQ=255)
                                    aGCCTGGAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATg                                                                                                              >  1:1243617/1‑150 (MQ=255)
                                      ccTGGAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCa                                                                                                            >  1:73550/1‑150 (MQ=255)
                                         ggAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTc                                                                                                         >  2:531244/1‑150 (MQ=255)
                                         ggAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTc                                                                                                         <  1:791526/150‑1 (MQ=255)
                                          gAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCt                                                                                                        >  2:358002/1‑150 (MQ=255)
                                            aGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTAt                                                                                                      <  1:546105/150‑1 (MQ=255)
                                                  tCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTaa                                                                                                >  1:417675/1‑150 (MQ=255)
                                                    aCTAAGATTATCACCAGCATGATGACGTGCTTCCTCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAAcg                                                                                              <  2:450436/150‑1 (MQ=255)
                                                    aCTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAAcg                                                                                              >  2:1377259/1‑150 (MQ=255)
                                                        agagTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTg                                                                                          <  2:1243617/150‑1 (MQ=255)
                                                        agagTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTg                                                                                          <  2:1127142/150‑1 (MQ=255)
                                                        agagTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTg                                                                                          <  2:73550/150‑1 (MQ=255)
                                                         gagTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGa                                                                                         <  2:833596/150‑1 (MQ=255)
                                                         gagTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGa                                                                                         >  1:278064/1‑150 (MQ=255)
                                                          agTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGaa                                                                                        <  1:863023/150‑1 (MQ=255)
                                                          agTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGaa                                                                                        <  2:773354/150‑1 (MQ=255)
                                                             aTCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATaa                                                                                     >  2:1219208/1‑150 (MQ=255)
                                                               cACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATaaaa                                                                                   >  2:144924/1‑150 (MQ=255)
                                                               cACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATaaaa                                                                                   >  1:211572/1‑150 (MQ=255)
                                                                 gggttATGATGAAGCACGTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATa                                                                                 >  1:612667/5‑150 (MQ=21)
                                                                       tgatgaCGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTcc                                                                           <  1:1219208/150‑1 (MQ=255)
                                                                          tgaCGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCctct                                                                        >  2:1159678/1‑150 (MQ=255)
                                                                           gaCGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCtg                                                                       >  1:1390538/1‑150 (MQ=255)
                                                                           gaCGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCtg                                                                       >  2:488969/1‑150 (MQ=255)
                                                                           gaCGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCtg                                                                       <  1:488969/150‑1 (MQ=255)
                                                                             cGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCtgtg                                                                     >  1:342198/1‑150 (MQ=255)
                                                                                  ttCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGattat                                                                <  1:117075/150‑1 (MQ=255)
                                                                                   tcatcaTAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGattatt                                                               >  1:186878/1‑150 (MQ=255)
                                                                                      tcatAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGATTATTGCg                                                            <  2:819504/150‑1 (MQ=255)
                                                                                      tcatAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGATTATTGCg                                                            <  1:1409603/150‑1 (MQ=255)
                                                                                      tcatAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGATTATTGCg                                                            <  2:674537/150‑1 (MQ=255)
                                                                                               tttCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCtg                                                                       <  2:1159514/130‑1 (MQ=255)
                                                                                               tttCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCtg                                                                       >  1:1159514/1‑130 (MQ=255)
                                                                                                ttCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATTGTTTCCTCTGTGATTATTGCGTAACGACTCg                                                  <  1:1113408/150‑1 (MQ=255)
                                                                                                 tCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGATTATTGCGTAACGACTCGc                                                 <  1:660506/150‑1 (MQ=255)
                                                                                                 tCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGATTATTGCGTAACGACTCGc                                                 >  2:1342672/1‑150 (MQ=255)
                                                                                                   cTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGATTATTGCGTAACGACTCGCaa                                               >  2:267510/1‑150 (MQ=255)
                                                                                                     tattaAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGATTATTGCGTAACGACTCGCAAAg                                             <  2:612667/150‑1 (MQ=255)
                                                                                                          aaaGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGATTATTGCGTAACGACTCGCAAAGAatat                                        >  1:1144319/1‑150 (MQ=255)
                                                                                                            aGCCCTCTTCTCCGGGAGAGGCATTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGATTATTGCGTAACGACTCGCAAAGAATAttt                                      <  1:763642/150‑1 (MQ=255)
                                                                                                               ccTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGATTATTGCGTAACGACTCGCAAAGAATATTTGCt                                   <  2:1051281/150‑1 (MQ=255)
                                                                                                                  cttctCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGATTATTGCGTAACGACTCGCAAAGAATATTTGCTGGc                                <  1:433143/150‑1 (MQ=255)
                                                                                                                   ttctCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGATTATTGCGTAACGACTCGCAAAGAATATTTGCTGGCa                               <  2:319452/150‑1 (MQ=255)
                                                                                                                   ttctCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGATTATTGCGTAACGACTCGCAAAGAATATTTGCTGGCa                               <  1:378248/150‑1 (MQ=255)
                                                                                                                    tctcCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGATTATTGCGTAACGACTCGCAAAGAATATTTGCTGGCAg                              >  2:695721/1‑150 (MQ=255)
                                                                                                                       ccGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGATTATTGCGt                                                           >  2:1174182/1‑118 (MQ=255)
                                                                                                                       ccGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGATTATTGCGt                                                           <  1:1174182/118‑1 (MQ=255)
                                                                                                                         gaagagggCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGATTATTGCGTAACGACTCGCAAAGAATATTTGCTGGCAGGATCg                         <  1:246720/144‑1 (MQ=255)
                                                                                                                           gagaGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGATTATTGCGTAACGACTCGCAAAGAATATTTGCTGGCAGGATCGCa                       <  2:211572/150‑1 (MQ=255)
                                                                                                                               ggCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGATTATTGCGTAACGACTCGCAAAGAATATTTGCTGGCAGGATCGCAGACt                   >  1:704765/1‑150 (MQ=255)
                                                                                                                                 cTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGATTATTGCGTAACGACTCGCAAAGAATATTTGCTGGCAGGATCGCAGACTAc                 >  1:1231800/1‑150 (MQ=255)
                                                                                                                                     cccGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGATTATTGCGTAACGACTCGCAAAGAATATTTGCTGGCAGGATCGCAGACTACAAAg             >  1:320163/1‑150 (MQ=255)
                                                                                                                                     cccGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGATTATTGCGTAACGACTCGCAAAGAATATTTGCTGGCAGGATCGCAGACTACAAAg             <  1:681520/150‑1 (MQ=255)
                                                                                                                                     atgatttCAGCGGGACGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGATTATTGCGTAACGACTCGCAAAGAATATTTGCTGGCAGGATCGCAGACTACAAAg             >  1:1364371/5‑150 (MQ=255)
                                                                                                                                        atttCAGCGGGACGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGATTATTGCGTAACGACTCGCAAAGAATATTTGCTGGCAGGATCGCAGACTACAAAGCCt          >  2:251542/2‑150 (MQ=255)
                                                                                                                                                gTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGATTATTGCGTAACGACTCGCAAAGAATATTTGCTGGCAGGATCGCAGACTACAAAGCCTGCGGATTg  <  1:440804/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GGGAGAGGTTTCCCCCTCCCCCTGGTGTTCTTAGTAAGCCTGGAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGATTATTGCGTAACGACTCGCAAAGAATATTTGCTGGCAGGATCGCAGACTACAAAGCCTGCGGATTG  >  NC_000913/2153319‑2153612

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: