Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,194,875 T→C 11.8% intergenic (‑10/+224) sibD ← / ← sibE sRNA antisense regulator of toxic IbsD protein/sRNA antisense regulator of toxic IbsE protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,194,8750TC11.8% 125.0 / 8.9 54intergenic (‑10/+224)sibD/sibEsRNA antisense regulator of toxic IbsD protein/sRNA antisense regulator of toxic IbsE protein
Reads supporting (aligned to +/- strand):  ref base T (20/25);  new base C (3/3);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.62e-01

CCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCA  >  NC_000913/3194728‑3195022
                                                                                                                                                   |                                                                                                                                                   
cccTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCctttacttt                                                                                                                                                    <  1:98855/149‑1 (MQ=255)
cccTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCctataattt                                                                                                                                                    >  1:208608/1‑149 (MQ=255)
 ccTCCCGGAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCctttacttta                                                                                                                                                   <  2:754969/149‑1 (MQ=255)
    cccGTAGAAGGGGGCCTTTCTGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGt                                                                                                                                                <  1:458870/149‑1 (MQ=255)
    cccGAAGAAGTGGGCCTTACAGAAGGAAAGGGTTATGAAGAAGCTCGTAATCATACTGATTGTGTTGTTACTCGTAAGTTTCTCAGCTTATTAACAGACAATCAGAGGGGGGAGAAATCCTCCCTCCCACTTGTTCCATTACTTTAgtt                                                                                                                                                >  2:964491/1‑147 (MQ=18)
              ggggCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGaaaaaaaa                                                                                                                                      <  2:838392/149‑1 (MQ=255)
              ggggCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGaaaaaaaa                                                                                                                                      <  1:853505/149‑1 (MQ=255)
                ggCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATc                                                                                                                                    >  2:107630/1‑149 (MQ=255)
                 gCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCa                                                                                                                                   >  2:988746/1‑149 (MQ=255)
                 gCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCa                                                                                                                                   <  2:621879/149‑1 (MQ=255)
                 gCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCa                                                                                                                                   >  2:591559/1‑149 (MQ=255)
                 gCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGGTGTTACTCGTAAGGTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCGCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCa                                                                                                                                   <  2:42074/149‑1 (MQ=255)
                   cctACAGCCGGAAAGGGTTAGTCTGAACCTCGTCATCATACTGATTGTGGAGTGACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAACTCCTCCCGCACCCTTGTGCCTTTACTTTAGGTTGAAAAAAAATCAGc                                                                                                                                 <  2:686079/147‑1 (MQ=255)
                          aaGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTGTAGGTTGGAAAAAAATCAGCGTCaaga                                                                                                                          >  2:968811/1‑149 (MQ=255)
                             gAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCTTAAGTTTCGCTTCTTATTAACAACAAATCAGAGGGGGGAGAAATCCTCCCTAAACCTTGTTCATTTACTATAGGTTGAAAAAACAACAGCGTCaataagc                                                                                                                       >  1:601617/1‑148 (MQ=11)
                               aaGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTAACCATATTTCCTGTACATAAGTATGAAAAAAACTCAGCGGCAAGAAGCtt                                                                                                                     >  1:397387/1‑149 (MQ=18)
                                   gTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCAt                                                                                                                 >  1:392938/1‑149 (MQ=255)
                                        gatgatGGTCGTCATACTACTGTTTGTTTTGTTACTCGTAAGTTTCGCTGCTTATTAACAGCGAATAAGAGGGGGGAGACATCCTCCCGCACGCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGTTGGTTTCATTAaca                                                                                                            <  2:122593/149‑1 (MQ=255)
                                             aGCTCGTCATCATACTGATTGTGTTGTTACTCGTTAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTACCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACatat                                                                                                       <  2:736762/149‑1 (MQ=255)
                                               cTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATcc                                                                                                     >  2:219098/1‑149 (MQ=255)
                                                   tcatcaTACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTgg                                                                                                 >  2:360458/1‑149 (MQ=255)
                                                    catcatACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGGc                                                                                                <  2:354706/149‑1 (MQ=255)
                                                       catACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGGCCCa                                                                                             >  1:411497/1‑149 (MQ=255)
                                                              aTTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGGCCCAAAAAGcc                                                                                      >  2:842250/1‑149 (MQ=255)
                                                              aTTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCACCCTCACCCTTGTTCATTTAATATAGGTTGAAAAAAAAACAGCGACAAGAAGGATTCATTAAAACATATCCCTGACACAAAAAGcc                                                                                      >  1:657609/1‑149 (MQ=255)
                                                               ttGGGTAGTTACTCGTAAGATTCGCAGCTTATTAACAGCCAATCAGAGGGGGTAGAAATCCTCCAGCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCATCGTCAAGAAGGTTTCATTAACACATATCCCTGGCCCAAAAAGCCg                                                                                     <  2:684394/149‑1 (MQ=255)
                                                                       ttACTCGTAAGTTTCGCAGCTAATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGGCCCAAAAAGCCGGGGAtttt                                                                             <  1:501392/149‑1 (MQ=255)
                                                                          cTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGGCCCAAAAAGCCGGGGATTTTTta                                                                          <  1:591559/149‑1 (MQ=255)
                                                                                 gTTTCGCAGCTTATTAACAGCCTATGAGATGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGc                                                                   <  2:929660/149‑1 (MQ=255)
                                                                                 gTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGc                                                                   <  1:253891/149‑1 (MQ=255)
                                                                                  tttCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCg                                                                  >  1:788212/1‑149 (MQ=255)
                                                                                     cGCAGCTTATTAACATCCATTCAGAGGGGGGATAAATCATCCCTCACGCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGACCCAAATAGCCGGGGATTTTTTATATCTGCGTTc                                                               <  2:657609/149‑1 (MQ=39)
                                                                                        aGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGc                                                            <  1:286102/149‑1 (MQ=255)
                                                                                          cTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTa                                                          >  2:374472/1‑149 (MQ=255)
                                                                                           ttattaACAGCCAATCAGAGGGGGGAGAAATCCTCCCTTACCCTTGTTCCTTTACTCTAGGTTGaa                                                                                                                                            <  2:610587/66‑1 (MQ=255)
                                                                                           ttattaACAGCCAATCAGAGGGGGGAGAAATCCTCCCTAACCCTTGTTCCTTTACTCTAGGTTGaa                                                                                                                                            >  1:610587/1‑66 (MQ=255)
                                                                                                aaCAGCCAATCAGAGGGGGGAGAAATCCTCCCTTACCCTTGTTCCTTTACTCTAgg                                                                                                                                                 <  2:754186/56‑1 (MQ=11)
                                                                                                aaCAGCCAATCAGAGGGGGGAGAAATCCTCCCTTACCCTTGTTCCTTTACTCTAGGt                                                                                                                                                <  2:885687/57‑1 (MQ=11)
                                                                                                aaCAGCCAATCAGAGGGGGGAGAAATCCTCCCTGACCCTTGTTCCTTTACTCTAgg                                                                                                                                                 >  1:754186/1‑56 (MQ=14)
                                                                                                aaCAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTCTAGGt                                                                                                                                                >  1:885687/1‑57 (MQ=255)
                                                                                                  cAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGc                                                  <  1:964491/149‑1 (MQ=255)
                                                                                                   aGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCa                                                 >  1:63333/1‑149 (MQ=255)
                                                                                                    gCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCaa                                                >  1:67978/1‑149 (MQ=255)
                                                                                                      cAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAAt                                              >  2:958980/1‑149 (MQ=255)
                                                                                                         tCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGcg                                           >  1:290872/1‑148 (MQ=255)
                                                                                                                   tagaAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTg                                 <  2:891590/148‑1 (MQ=255)
                                                                                                                   gagaAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGTTTTCATTAACACATATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCTTTCCGCTAAAAGGTGAAAATGCTCAGGCCGTTg                                 >  2:240426/1‑149 (MQ=255)
                                                                                                                     taaaTCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCa                               <  1:432486/148‑1 (MQ=255)
                                                                                                                     gaAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCa                               <  1:438900/149‑1 (MQ=255)
                                                                                                                     gaAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCa                               <  2:739122/149‑1 (MQ=255)
                                                                                                                         tcctccCTCACCCTTGTTCCTATACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGGCCCAAAAAGCCGGGGAGTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGt                           >  1:780230/1‑149 (MQ=255)
                                                                                                                          cctccctcACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGtt                          >  1:653079/1‑149 (MQ=255)
                                                                                                                           ctccctcACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGttt                         <  2:67978/149‑1 (MQ=255)
                                                                                                                            tccctcACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGtttt                        <  2:411497/149‑1 (MQ=255)
                                                                                                                               ctcACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCg                     <  1:89815/149‑1 (MQ=255)
                                                                                                                                gcACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGt                    <  2:1015175/148‑1 (MQ=255)
                                                                                                                                      ttGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACTTTGCAGCGTTTTGCGTGACCGc              >  1:474821/1‑149 (MQ=255)
                                                                                                                                            ctttactttaGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCgggg        >  1:565672/1‑149 (MQ=255)
                                                                                                                                            ctttactttaGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCGGCGTTCCGCTAAAAGGTGCAAATGCTGAGGACGTTGCAGCGTTTTGCGTGACCGCTCgggg        >  2:901267/1‑149 (MQ=255)
                                                                                                                                              ttactttaGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGaa      <  1:360458/149‑1 (MQ=255)
                                                                                                                                                 ctttaGGTTGAAAAAAAATCAGCGTCTAGAAGGTTTCATTAACACATATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCGAGGGGAAGGc   >  1:416925/1‑149 (MQ=255)
                                                                                                                                                  tttaGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCa  <  2:392938/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCA  >  NC_000913/3194728‑3195022

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: