Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,814,861 C→T 13.7% noncoding (60/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,8610CT13.7% 222.3 / 21.0 80noncoding (60/78 nt)micAsRNA antisense regulator of ompA, lamB, ompX, and phoP; Hfq‑dependent
Reads supporting (aligned to +/- strand):  ref base C (31/38);  new base T (7/4);  total (38/42)
Fisher's exact test for biased strand distribution p-value = 3.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGCT  >  NC_000913/2814713‑2815007
                                                                                                                                                    |                                                                                                                                                  
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                                                                                                                                             >  3:705542/1‑149 (MQ=255)
       tCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGc                                                                                                                                             <  4:790596/149‑1 (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:532471/149‑1 (MQ=255)
                                                         aaCTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTc                                                                                           >  4:221451/1‑149 (MQ=255)
                                                          aCTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTc                                                                                           <  1:876828/148‑1 (MQ=255)
                                                           ctctTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGt                                                                                                                                                <  1:212663/94‑1 (MQ=255)
                                                           ctctTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGt                                                                                                                                                >  2:212663/1‑94 (MQ=255)
                                                           ctctTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttcttttc                                                                                                                                 >  1:771436/1‑109 (MQ=255)
                                                           ctctTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttcttttc                                                                                                                                 <  2:771436/109‑1 (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                                                                                                                                       >  2:555267/1‑75 (MQ=255)
                                                                                       atGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCttttt                                                                                                                                       <  1:555267/75‑1 (MQ=255)
                                                                                                      ttGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTTCTGACGGCGTTCAGACGCCt                                               >  1:575997/1‑148 (MQ=255)
                                                                                                         ttATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCAttt                                                                       >  2:713137/1‑121 (MQ=255)
                                                                                                         ttATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCAttt                                                                       <  1:713137/121‑1 (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                                                                                           >  4:415230/4‑70 (MQ=38)
                                                                                                                                        aaaggCCACTCGTGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTc                                                                                           <  3:415230/67‑1 (MQ=38)
                                                                                                                                         ttGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCCAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAAttt           <  4:622108/149‑1 (MQ=255)
                                                                                                                                         aaggCCACTCGTGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACgg                                                            <  1:768536/98‑1 (MQ=255)
                                                                                                                                         aaggCCACTCGTGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACgg                                                            >  2:768536/3‑100 (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         <  1:357936/148‑1 (MQ=255)
                                                                                                                                            gCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTcc         >  2:119539/1‑148 (MQ=255)
                                                                                                                                            gCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCa        <  3:466684/149‑1 (MQ=255)
                                                                                                                                             ccACTCGTGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCa        >  1:510761/1‑148 (MQ=255)
                                                                                                                                             ccACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGtt                                                                                            <  4:340952/64‑1 (MQ=255)
                                                                                                                                             ccACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGtt                                                                                            >  3:340952/1‑64 (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)
                                                                                                                                                    |                                                                                                                                                  
GAAGCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGCT  >  NC_000913/2814713‑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: