Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,153,466 C→A 6.1% intergenic (+12/‑180) 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,4660CA 6.1% 268.4 / 4.3 81intergenic (+12/‑180)sibA/sibBsmall RNA SibA/small RNA SibB
Reads supporting (aligned to +/- strand):  ref base C (34/42);  new base A (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 = 8.56e-01

GGGAGAGGTTTCCCCCTCCCCCTGGTGTTCTTAGTAAGCCTGGAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGATTATTGCGTAACGACTCGCAAAGAATATTTGCTGGCAGGATCGCAGACTACAAAGCCTGCGGATTGACA  >  NC_000913/2153319‑2153615
                                                                                                                                                   |                                                                                                                                                     
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:168326/150‑1 (MQ=255)
              cctccccctGGTGTTCTTAGTAAGCCTGGAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGc                                                                                                                                       <  1:991357/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                                                                                                            <  1:791526/150‑1 (MQ=255)
                                         ggAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTc                                                                                                            >  2:531244/1‑150 (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                                                                                            >  1:278064/1‑150 (MQ=255)
                                                         gagTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGa                                                                                            <  2:833596/150‑1 (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                                                                          >  2:488969/1‑150 (MQ=255)
                                                                           gaCGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCtg                                                                          <  1:488969/150‑1 (MQ=255)
                                                                           gaCGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCtg                                                                          >  1:1390538/1‑150 (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:674537/150‑1 (MQ=255)
                                                                                      tcatAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGATTATTGCg                                                               <  1:1409603/150‑1 (MQ=255)
                                                                                      tcatAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGATTATTGCg                                                               <  2:819504/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                                                    >  2:1342672/1‑150 (MQ=255)
                                                                                                 tCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGATTATTGCGTAACGACTCGc                                                    <  1:660506/150‑1 (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                                                              <  1:1174182/118‑1 (MQ=255)
                                                                                                                       ccGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGATTATTGCGt                                                              >  2:1174182/1‑118 (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)
                                                                                                                                                  cccGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGATTATTGCGTAACGACTCGCAAAGAATATTTGCTGGCAGGATCGCAGACTACAAAGCCTGCGGATTGAc   >  2:1359460/1‑150 (MQ=255)
                                                                                                                                                  cccGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGATTATTGCGTAACGACTCGCAAAGAATATTTGCTGGCAGGATCGCAGACTACAAAGCCTGCGGATTGAc   <  2:1390538/150‑1 (MQ=255)
                                                                                                                                                   ccGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGATTATTGCGTAACGACTCGCAAAGAATATTTGCTGGCAGGATCGCAGACTACAAAGCCTGCGGATTGACa  >  2:1060857/1‑150 (MQ=255)
                                                                                                                                                   ccGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGATTATTGCGTAACGACTCGCAAAGAATATTTGCTGGCAGGATCGCAGACTACAAAGCCTGCGGATTGACa  >  1:1159883/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGGAGAGGTTTCCCCCTCCCCCTGGTGTTCTTAGTAAGCCTGGAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCCTCTGTGATTATTGCGTAACGACTCGCAAAGAATATTTGCTGGCAGGATCGCAGACTACAAAGCCTGCGGATTGACA  >  NC_000913/2153319‑2153615

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: