Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,814,860 A→G 13.7% noncoding (59/78 nt) micA → sRNA antisense regulator of ompA, lamB, ompX, and phoP; Hfq‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,814,8600AG13.7% 216.8 / 20.9 80noncoding (59/78 nt)micAsRNA antisense regulator of ompA, lamB, ompX, and phoP; Hfq‑dependent
Reads supporting (aligned to +/- strand):  ref base A (32/37);  new base G (7/4);  total (39/41)
Fisher's exact test for biased strand distribution p-value = 3.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

TGAAGCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGCT  >  NC_000913/2814712‑2815007
                                                                                                                                                    |                                                                                                                                                   
tgAAGCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCa                                                                                                                                                     >  3:413817/1‑149 (MQ=255)
 gAAGCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCAc                                                                                                                                                    >  4:384408/1‑149 (MQ=255)
  aaGCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACg                                                                                                                                                   >  4:422010/1‑149 (MQ=255)
   aGCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACg                                                                                                                                                   <  2:686110/148‑1 (MQ=255)
      tATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGt                                                                                                                                                <  1:135553/148‑1 (MQ=255)
       aTCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTgg                                                                                                                                              <  4:879636/149‑1 (MQ=255)
        tCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGTCGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGc                                                                                                                                             >  3:111777/1‑149 (MQ=255)
        tCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTgg                                                                                                                                              <  1:374750/148‑1 (MQ=255)
        tCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGc                                                                                                                                             <  4:790596/149‑1 (MQ=255)
        tCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGc                                                                                                                                             >  3:705542/1‑149 (MQ=255)
            acaacGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCtt                                                                                                                                          >  1:887080/1‑148 (MQ=255)
              aacGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCttttt                                                                                                                                       <  4:151600/149‑1 (MQ=255)
                    aTTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttctttt                                                                                                                                  <  2:854323/148‑1 (MQ=255)
                         gCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTc                                                                                                                             <  1:902174/148‑1 (MQ=255)
                          ccACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAg                                                                                                                           <  4:413817/149‑1 (MQ=255)
                          ccACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAg                                                                                                                           <  3:384408/149‑1 (MQ=255)
                              cTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCg                                                                                                                        >  2:411433/1‑148 (MQ=255)
                                ccGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCgtg                                                                                                                      >  2:867450/1‑148 (MQ=255)
                                  ggTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGttt                                                                                                                   <  4:705542/149‑1 (MQ=255)
                                   gTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTCGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGttt                                                                                                                   <  2:733637/148‑1 (MQ=255)
                                       tttttttAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCa                                                                                                              <  3:422010/149‑1 (MQ=255)
                                       tttttttAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCa                                                                                                              <  4:111777/149‑1 (MQ=255)
                                        ttttttAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCa                                                                                                              <  1:411433/148‑1 (MQ=255)
                                            ttAAAAATTTTCTGAACTCCTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCcac                                                                                                         <  4:62442/149‑1 (MQ=255)
                                                        tGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGt                                                                                             >  3:532471/1‑149 (MQ=255)
                                                         gAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGtt                                                                                            >  4:466684/1‑149 (MQ=255)
                                                          aaCTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGtt                                                                                            <  2:887080/148‑1 (MQ=255)
                                                          aaCTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTc                                                                                           >  4:221451/1‑149 (MQ=255)
                                                          aaCTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTc                                                                                           <  4:532471/149‑1 (MQ=255)
                                                           aCTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTc                                                                                           <  1:876828/148‑1 (MQ=255)
                                                            ctctTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGt                                                                                                                                                <  1:212663/94‑1 (MQ=255)
                                                            ctctTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGt                                                                                                                                                >  2:212663/1‑94 (MQ=255)
                                                            ctctTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttcttttc                                                                                                                                 <  2:771436/109‑1 (MQ=255)
                                                            ctctTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttcttttc                                                                                                                                 >  1:771436/1‑109 (MQ=255)
                                                               tttcttCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTAt                                                                                      >  4:299109/1‑149 (MQ=255)
                                                               tttcttCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTAt                                                                                      >  3:43043/1‑149 (MQ=255)
                                                                ttcttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTAt                                                                                      >  1:853984/1‑148 (MQ=255)
                                                                             aGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCAt                                                                         <  1:867450/148‑1 (MQ=255)
                                                                                cTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTc                                                                      <  2:798282/148‑1 (MQ=255)
                                                                                cTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTc                                                                      >  1:492232/1‑148 (MQ=255)
                                                                                 tGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCct                                                                    >  4:395586/1‑149 (MQ=255)
                                                                                   aGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCctgc                                                                  <  3:221451/149‑1 (MQ=255)
                                                                                        atGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCttttt                                                                                                                                       <  1:555267/75‑1 (MQ=255)
                                                                                        atGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCttttt                                                                                                                                       >  2:555267/1‑75 (MQ=255)
                                                                                                       ttGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTTCTGACGGCGTTCAGACGCCt                                               >  1:575997/1‑148 (MQ=255)
                                                                                                          ttATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCAttt                                                                       <  1:713137/121‑1 (MQ=255)
                                                                                                          ttATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCAttt                                                                       >  2:713137/1‑121 (MQ=255)
                                                                                                           tATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCGTGCTGACGGCGTTCAGACGCCTCgcgc                                          >  3:622108/1‑149 (MQ=255)
                                                                                                              catcatcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCg                                       <  3:679851/149‑1 (MQ=255)
                                                                                                                      cTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACaa                                >  1:347390/1‑148 (MQ=255)
                                                                                                                       tGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACaaaa                              <  4:43043/149‑1 (MQ=255)
                                                                                                                       tGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACaaa                               >  2:357936/1‑148 (MQ=255)
                                                                                                                        gAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACaaaa                              >  1:91923/1‑148 (MQ=255)
                                                                                                                         aaTTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAAt                             <  2:91923/148‑1 (MQ=255)
                                                                                                                         aaTTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATc                            <  3:299109/149‑1 (MQ=255)
                                                                                                                          aTTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATc                            >  1:48036/1‑148 (MQ=255)
                                                                                                                            tCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCct                          <  2:853984/148‑1 (MQ=255)
                                                                                                                             cAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCctc                         >  1:623488/1‑148 (MQ=255)
                                                                                                                                       aTTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTACGCGCAGa               <  2:89426/148‑1 (MQ=255)
                                                                                                                                         aaaggCCACTCGTGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTc                                                                                           <  3:415230/67‑1 (MQ=38)
                                                                                                                                         aaaggCCACTCGTGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTc                                                                                           >  4:415230/4‑70 (MQ=38)
                                                                                                                                          ttGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCCAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAAttt           <  4:622108/149‑1 (MQ=255)
                                                                                                                                          aaggCCACTCGTGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACgg                                                            >  2:768536/3‑100 (MQ=255)
                                                                                                                                          aaggCCACTCGTGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACgg                                                            <  1:768536/98‑1 (MQ=255)
                                                                                                                                          aaggCCACTCGTGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATc                            >  3:811503/3‑132 (MQ=255)
                                                                                                                                          aaggCCACTCGTGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATc                            <  4:811503/130‑1 (MQ=255)
                                                                                                                                          aaggCCACTCGTGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAAttt           >  3:698974/3‑149 (MQ=255)
                                                                                                                                           aggCCACTCGTGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGGGGCTTc                                                                            >  4:351330/2‑83 (MQ=255)
                                                                                                                                           aggCCACTCGTGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTc                                                                            <  3:351330/82‑1 (MQ=255)
                                                                                                                                            ggCCACTCGTGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTc          >  2:787786/1‑148 (MQ=255)
                                                                                                                                             gCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTcc         >  2:119539/1‑148 (MQ=255)
                                                                                                                                             gCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTcc         <  1:357936/148‑1 (MQ=255)
                                                                                                                                             gCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCa        <  3:466684/149‑1 (MQ=255)
                                                                                                                                              ccACTCGTGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCa        >  1:510761/1‑148 (MQ=255)
                                                                                                                                              ccACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGtt                                                                                            >  3:340952/1‑64 (MQ=255)
                                                                                                                                              ccACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGtt                                                                                            <  4:340952/64‑1 (MQ=255)
                                                                                                                                              ccACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAg       >  3:513115/1‑149 (MQ=255)
                                                                                                                                                  tCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCgg    >  2:815247/1‑148 (MQ=255)
                                                                                                                                                  tCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGc   <  3:395586/149‑1 (MQ=255)
                                                                                                                                                   cACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGc   >  1:399504/1‑148 (MQ=255)
                                                                                                                                                   cACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGc   >  2:53727/1‑148 (MQ=255)
                                                                                                                                                    aCGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGct  <  2:575997/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGAAGCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGCT  >  NC_000913/2814712‑2815007

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