Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,350,682:1 +A 5.1% noncoding (106/121 nt) arcZ → sRNA positive antisense regulator of rpoS; binds Hfq

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,350,6821.A 5.1% 668.0 / 30.4 156noncoding (106/121 nt)arcZsRNA positive antisense regulator of rpoS; binds Hfq
Reads supporting (aligned to +/- strand):  ref base . (70/78);  new base A (2/6);  total (72/84)
Fisher's exact test for biased strand distribution p-value = 2.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-01

Only 100 of 158 total aligned reads displayed.

TTCATGTAACAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGT  >  NC_000913/3350537‑3350830
                                                                                                                                                  |                                                                                                                                                    
 tCATGTAACAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCTCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCgg                                                                                                                                                  <  3:131621/149‑1 (MQ=255)
 tCATGTAACAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCgg                                                                                                                                                  <  1:538435/149‑1 (MQ=255)
   aTGTAACAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCgggg                                                                                                                                                <  4:412518/149‑1 (MQ=255)
    tGTAACAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGG‑CTAGACCGGGGt                                                                                                                                               <  2:900731/149‑1 (MQ=255)
     gTAACAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTc                                                                                                                                              >  1:161984/1‑149 (MQ=255)
        aCAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGG‑CTAGACCGGGGt                                                                                                                                               >  1:1474133/1‑145 (MQ=255)
         cAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCAttt                                                                                                                                          >  1:1419841/1‑149 (MQ=255)
           aaTCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCAttttt                                                                                                                                        >  2:1461683/1‑149 (MQ=255)
           aaTCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCAttttt                                                                                                                                        >  2:1568018/1‑149 (MQ=255)
           aaTCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCAttttt                                                                                                                                        >  1:701267/1‑149 (MQ=255)
             tCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCAtttt                                                                                                                                         >  2:1393438/1‑146 (MQ=255)
             tCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTa                                                                                                                                      <  4:149463/149‑1 (MQ=255)
             tCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTa                                                                                                                                      <  3:175158/149‑1 (MQ=255)
             tCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTa                                                                                                                                      >  4:1386739/1‑149 (MQ=255)
             tCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTa                                                                                                                                      >  2:95841/1‑149 (MQ=255)
               aTTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGt                                                                                                                                    >  2:809691/1‑149 (MQ=255)
               aTTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGt                                                                                                                                    >  3:396849/1‑149 (MQ=255)
                tttAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCgggg                                                                                                                                                <  4:100167/136‑1 (MQ=255)
                tttAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCgggg                                                                                                                                                >  3:100167/1‑136 (MQ=255)
                tttAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTg                                                                                                                                   >  1:1332373/1‑149 (MQ=255)
                    aGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCtt                                                                                                                               <  1:658353/149‑1 (MQ=255)
                     gTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCttt                                                                                                                              >  1:1418979/1‑149 (MQ=255)
                          gCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑TCGGGGTCATTTTTTAGTGGCTTTTGcca                                                                                                                         <  4:1043290/149‑1 (MQ=255)
                             aTCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGccaccc                                                                                                                      <  3:932435/149‑1 (MQ=255)
                               cTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGccacccac                                                                                                                    <  3:469341/149‑1 (MQ=255)
                                 tAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGc                                                                                                                  >  1:736541/1‑149 (MQ=255)
                                    cTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCttt                                                                                                               <  2:604640/149‑1 (MQ=255)
                                    cTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCACCATAATAATTTACGGCGCAGCCAAGATTTCTCTGTTGTTGGCGCAGTATTCCCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCttt                                                                                                               >  3:818845/1‑149 (MQ=255)
                                     tgcgtgcgGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGccaccca                                                                                                                     >  1:599651/1‑142 (MQ=255)
                                     tgcgtgcgGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGccaccca                                                                                                                     <  2:599651/142‑1 (MQ=255)
                                     tgcgtgcgGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTc                                                                                                              >  2:1217886/1‑149 (MQ=255)
                                      gcgtgcgGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCa                                                                                                             <  3:1132560/149‑1 (MQ=255)
                                          gcgGCCTGAAAAACAGTGATGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTACCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACCCTTTCAGCAc                                                                                                         <  2:1099964/149‑1 (MQ=255)
                                                tGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTAc                                                                                                   >  3:1267462/1‑149 (MQ=255)
                                                tGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTAc                                                                                                   <  2:161984/149‑1 (MQ=255)
                                                tGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTAc                                                                                                   >  4:1329771/1‑149 (MQ=255)
                                                 gAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTAcg                                                                                                  >  2:1251027/1‑149 (MQ=255)
                                                   aaaaCAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTAcgtc                                                                                                <  3:1442725/149‑1 (MQ=255)
                                                   aaaaCAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTAcgtc                                                                                                <  3:905535/149‑1 (MQ=255)
                                                     aaCAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCtttt                                                                                                                             >  1:913312/1‑118 (MQ=255)
                                                     aaCAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTAcgtcgt                                                                                              <  1:1108104/149‑1 (MQ=255)
                                                     aaCAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTAcgacgt                                                                                              >  2:502735/1‑149 (MQ=255)
                                                       cactGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGCGa                                                                                            <  3:1386739/146‑1 (MQ=255)
                                                            cTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCa                                                                                       >  3:868439/1‑149 (MQ=255)
                                                             tgtgCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCAc                                                                                      <  4:935142/149‑1 (MQ=255)
                                                             tgtgCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCAc                                                                                      >  1:458790/1‑149 (MQ=255)
                                                                  ccTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACtcttc                                                                                 <  1:405204/149‑1 (MQ=255)
                                                                   cTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACtcttct                                                                                <  2:1419841/149‑1 (MQ=255)
                                                                    ttGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACtcttct                                                                                <  1:272057/148‑1 (MQ=255)
                                                                    ttGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACtcttct                                                                                >  2:272057/1‑148 (MQ=255)
                                                                     tGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCttt                                                                              <  1:328439/149‑1 (MQ=255)
                                                                      gTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTc                                                                             <  1:1217886/149‑1 (MQ=255)
                                                                        aaCTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCAt                                                                                                                                            >  2:368948/1‑84 (MQ=255)
                                                                             atcatAATAATTTACGGCTCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCtt                                                                      >  1:1374271/1‑149 (MQ=255)
                                                                             atcatAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCtt                                                                      <  1:1085977/149‑1 (MQ=255)
                                                                             atcatAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCtt                                                                      >  1:26011/1‑149 (MQ=255)
                                                                             atcatAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCtt                                                                      <  1:1276418/149‑1 (MQ=255)
                                                                              tcatAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCtt                                                                      >  3:498824/1‑148 (MQ=255)
                                                                              tcatAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCtt                                                                      <  4:498824/148‑1 (MQ=255)
                                                                              tcatAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTc                                                                     >  3:193280/1‑149 (MQ=255)
                                                                                ataataATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCaa                                                                   >  1:262153/1‑149 (MQ=255)
                                                                                 taataaTTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCttt                                                                              >  4:809379/1‑137 (MQ=255)
                                                                                 taataaTTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAAt                                                                  >  3:5083/1‑149 (MQ=255)
                                                                                 taataaTTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCACTTCAAt                                                                  >  1:1531347/1‑149 (MQ=255)
                                                                                  aataatTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTcttctt                                                                               <  2:625394/135‑1 (MQ=255)
                                                                                  aataatTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATc                                                                 >  2:402002/1‑149 (MQ=255)
                                                                                      atTTACGGCGCAGCCGAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCCGTCTAG‑CCGGGGTCATTTTTTAATGGCTTTTGCCACCCACGCTTTCAGCACTTCTACCTCGTGAAGCCACTCTCCTTTCATCTCTTCAATCCAtt                                                             >  4:164552/1‑149 (MQ=255)
                                                                                                cAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGtt                                                   <  4:495342/149‑1 (MQ=255)
                                                                                                   ccAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATc                                                >  3:1426725/1‑149 (MQ=255)
                                                                                                    cAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCtt                                                                      >  3:63049/1‑126 (MQ=255)
                                                                                                    cAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTc                                                            >  3:814263/1‑136 (MQ=255)
                                                                                                    cAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTc                                                            <  4:814263/136‑1 (MQ=255)
                                                                                                    cAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCt                                               >  3:403406/1‑149 (MQ=255)
                                                                                                     aaGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAg                                                                                                            <  3:670380/87‑1 (MQ=255)
                                                                                                     aaGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAg                                                                                                            >  4:670380/1‑87 (MQ=255)
                                                                                                     aaGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCtt                                              >  2:112452/1‑149 (MQ=255)
                                                                                                                   gtTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAgg                                >  1:495595/1‑149 (MQ=255)
                                                                                                                    ttGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGt                               <  1:402002/149‑1 (MQ=255)
                                                                                                                    ttGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGt                               >  1:1371329/1‑149 (MQ=255)
                                                                                                                     tGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTc                              <  4:1404490/149‑1 (MQ=255)
                                                                                                                      ggCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCa                             <  4:396849/149‑1 (MQ=255)
                                                                                                                       gcgcAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGa                                                                                            >  1:614646/1‑85 (MQ=255)
                                                                                                                       gcgcAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGa                                                                                            <  2:614646/85‑1 (MQ=255)
                                                                                                                           aGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCACTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGgaga                        <  2:1531347/149‑1 (MQ=255)
                                                                                                                            gTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGAc                       <  1:913661/149‑1 (MQ=255)
                                                                                                                              aTTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTg                     >  2:924158/1‑149 (MQ=255)
                                                                                                                               ttCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGa                    >  1:1224618/1‑149 (MQ=255)
                                                                                                                               ttCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGa                    >  1:647334/1‑149 (MQ=255)
                                                                                                                                 cgcgCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAAt                  <  4:5083/149‑1 (MQ=255)
                                                                                                                                  gcgcACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAAtt                 >  2:744559/1‑149 (MQ=255)
                                                                                                                                  gcgcACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAAtt                 >  4:1261285/1‑149 (MQ=255)
                                                                                                                                  gcgcACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAAtt                 <  4:724954/149‑1 (MQ=255)
                                                                                                                                     cACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGc              <  1:1251027/149‑1 (MQ=255)
                                                                                                                                     cACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGc              <  2:194713/149‑1 (MQ=255)
                                                                                                                                     cACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGc              <  2:1341961/149‑1 (MQ=255)
                                                                                                                                      aCCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCt             <  3:335455/149‑1 (MQ=255)
                                                                                                                                      aCCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCt             <  4:1267462/149‑1 (MQ=255)
                                                                                                                                        cccGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGa           >  1:400137/1‑149 (MQ=255)
                                                                                                                                                 g‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGt  <  1:972378/149‑1 (MQ=255)
                                                                                                                                                 g‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGt  <  2:1371329/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TTCATGTAACAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGT  >  NC_000913/3350537‑3350830

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