Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,152,247 +C intergenic (+10/+29) yghB → / ← yqhC hypothetical protein/putative ARAC‑type regulatory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,152,2441.C99.4% 529.5 / ‑3.8 166intergenic (+7/+32)yghB/yqhChypothetical protein/putative ARAC‑type regulatory protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (106/59);  total (107/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.96e-01

Only 100 of 198 total aligned reads displayed.

TGATGATCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCGAA  >  NC_000913/3152152‑3152332
                                                                                              |                                                                                          
tgatgaTCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGAtttt                                                                                             >  1:3843937/1‑93 (MQ=255)
tgatgaTCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGAtttt                                                                                             >  1:3940926/1‑93 (MQ=255)
 gatgatCCTGCCAATTGCCTTGTTAACCGCTGGCTTGGTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGAGGATTttc                                                                                           >  2:2819139/1‑93 (MQ=255)
   tgatCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTttccc                                                                                         <  1:3620705/93‑1 (MQ=255)
   tgatCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTttccc                                                                                         >  2:4365894/1‑93 (MQ=255)
   tgatCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTttccc                                                                                         >  1:4357498/1‑93 (MQ=255)
    gatCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTttcccc                                                                                        >  1:4402806/1‑92 (MQ=255)
    gatCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTttcccc                                                                                        <  2:4444581/93‑2 (MQ=255)
     atCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTttccccg                                                                                       <  2:3935952/93‑3 (MQ=255)
      tCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTttccccgt                                                                                      >  2:3023515/1‑90 (MQ=255)
         tGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTtcccc‑‑gttcc                                                                                  >  2:3694614/1‑93 (MQ=255)
         tGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTtcccc‑‑gttcc                                                                                  <  1:887328/93‑1 (MQ=255)
             aaTTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGt                                                                              <  2:2235226/93‑1 (MQ=255)
              aTTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGtt                                                                             >  1:283793/1‑93 (MQ=255)
              aTTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGtt                                                                             <  2:1577110/93‑1 (MQ=255)
                tGCGTTGTTAACCGCTGTCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGc                                                                           <  2:468534/93‑1 (MQ=255)
                tGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATAAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGc                                                                           <  2:149318/93‑1 (MQ=255)
                 gCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCt                                                                          >  2:2801959/1‑93 (MQ=255)
                  ccTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTg                                                                         >  1:4152542/1‑93 (MQ=255)
                       ttAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCg                                                                    <  2:1306412/93‑1 (MQ=255)
                        tAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCgg                                                                   >  1:4580227/1‑93 (MQ=255)
                          aCCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGa                                                                 >  1:1017940/1‑93 (MQ=255)
                           ccGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGaa                                                                <  1:2930966/93‑1 (MQ=255)
                            cGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAAc                                                               >  1:1029515/1‑93 (MQ=255)
                            cGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAAc                                                               >  2:3679194/1‑93 (MQ=255)
                             gCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACg                                                              >  1:2998736/1‑93 (MQ=255)
                             gCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACg                                                              >  1:1086210/1‑93 (MQ=255)
                               tGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGCACGtt                                                            >  1:651691/1‑92 (MQ=255)
                               tGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTa                                                            >  2:4176495/1‑93 (MQ=255)
                                ggCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTAt                                                           >  1:356722/1‑93 (MQ=255)
                                ggCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTAt                                                           >  2:2888937/1‑93 (MQ=255)
                                ggCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTACTGTACCGGGAACGTAt                                                           >  1:4125847/1‑93 (MQ=255)
                                 gCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTAtt                                                          <  1:2209332/93‑1 (MQ=255)
                                 gCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTAtt                                                          <  1:670339/93‑1 (MQ=255)
                                 gCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCCGTACCGGGAACGTAtt                                                          >  2:4393187/1‑93 (MQ=255)
                                  cTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTAttt                                                         >  1:199170/1‑93 (MQ=255)
                                    tGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTaa                                                       >  2:1200448/1‑93 (MQ=255)
                                      ttAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAAtt                                                     >  2:4554060/1‑93 (MQ=255)
                                         ggCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTccc                                                  >  2:3665556/1‑93 (MQ=255)
                                         ggCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTccc                                                  >  1:1527744/1‑93 (MQ=255)
                                          gCACGCTGTTTGTTGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTcccc                                                 <  2:4599028/93‑1 (MQ=255)
                                              gCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCa                                             >  2:298069/1‑93 (MQ=255)
                                               gtGTTTGTGGTTATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCAt                                            <  1:180076/92‑1 (MQ=255)
                                               cTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCAt                                            <  2:337384/93‑1 (MQ=255)
                                                tGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATc                                           >  2:2211370/1‑93 (MQ=255)
                                                tGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATc                                           <  1:2620361/93‑1 (MQ=255)
                                                 gTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCg                                          >  1:2715036/1‑93 (MQ=255)
                                                     gtggtgATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCg                                      >  2:1906178/1‑93 (MQ=255)
                                                     gtggtgATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCg                                      >  2:3587641/1‑93 (MQ=255)
                                                      tggtgATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGc                                     >  1:490980/1‑93 (MQ=255)
                                                         tgATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCAtt                                  >  1:2108532/1‑93 (MQ=255)
                                                         tgATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCAtt                                  >  1:3393654/1‑93 (MQ=255)
                                                          gATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTc                                 >  2:4615943/1‑93 (MQ=255)
                                                          gATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTc                                 >  1:2667930/1‑93 (MQ=255)
                                                            ttAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGGACGCATTTAATTCCCCTGCATCGGCCGCATTCtt                               >  2:1835136/1‑93 (MQ=255)
                                                            ttAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCtt                               >  2:598059/1‑93 (MQ=255)
                                                            ttAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GATCCCGGCTTCTTTAGCGGGGACCTATTTATATGTCCTGCCTCGCCCGCACTCtt                               >  2:1655435/1‑93 (MQ=255)
                                                             tAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCTCCTGCATCGCCCGCATTCTTg                              >  2:506189/1‑93 (MQ=255)
                                                             tAAAAAAAAATACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTg                              >  2:2729611/1‑93 (MQ=255)
                                                              aaaaaaaaaTACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGc                             <  1:3223347/93‑1 (MQ=255)
                                                              aaaaaaaaaTACTGTAACGCCTGACGATTTTCCCC‑‑CTTCCCGGTTGCTGTACCGGGAACGCATTTAATTCCCCTGCACCGCCCGCATTCCTGc                             >  2:1106101/1‑93 (MQ=255)
                                                               aaaaaaaaTACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGcc                            >  2:3860136/1‑93 (MQ=255)
                                                               aaaaaaaaTACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGcc                            >  1:1192483/1‑93 (MQ=255)
                                                               aaaaaaaaTACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGcc                            <  2:2022169/93‑1 (MQ=255)
                                                               aaaaaaaaTACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGcc                            >  1:776672/1‑93 (MQ=255)
                                                               aaaaaaaaTACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATGCCCCTGCATCGCCCGCATTCTTGcc                            <  2:3325133/93‑1 (MQ=255)
                                                                aaaaaaaTACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCCGGGACGGATTTAATTTCCCTGCATCGCCCGCCTTCTTGccc                           >  2:821498/1‑92 (MQ=255)
                                                                 aaaaaaTACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGc                          >  1:1109992/1‑93 (MQ=255)
                                                                 aaaaaaTACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGc                          >  1:2754945/1‑93 (MQ=255)
                                                                    aaaTACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATc                       <  2:1300777/93‑1 (MQ=255)
                                                                    aaaTACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATc                       >  1:591535/1‑93 (MQ=255)
                                                                     aaTACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCt                      >  2:2286838/1‑93 (MQ=255)
                                                                     aaTACTGTAACGCCTGACGATTTTCCCCGGTTCCCGGTGGTGTACCGGGGACGTATTTAATTCCCCCGCATCGCCCGCCATCTTTGCGCATCt                      >  2:3964524/1‑93 (MQ=255)
                                                                      aTACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCtt                     >  1:2593441/1‑93 (MQ=255)
                                                                       tACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTc                    >  1:1864246/1‑93 (MQ=255)
                                                                       tACTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTc                    >  1:4235865/1‑93 (MQ=255)
                                                                        aCTGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTCTACCGGGGACGGATTTAATTCCCCTGCATCGCCCGCATTCCTTCCGCCTCtttc                   >  2:2533812/1‑91 (MQ=255)
                                                                          tGTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTcccc                 <  2:4307012/93‑1 (MQ=255)
                                                                           gTAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTccccc                >  2:1465588/1‑93 (MQ=255)
                                                                            tAACGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCg               >  2:792451/1‑93 (MQ=255)
                                                                             aaCGCCTGACGATGTACCCC‑‑GTTCCCGGTTGCTGTACCGGGACTGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCGTCCTCCCCCgg              <  2:4572929/93‑1 (MQ=255)
                                                                              aCGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGc             <  2:2396492/93‑1 (MQ=255)
                                                                               cGCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCg            <  1:2978616/93‑1 (MQ=255)
                                                                                gCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGGACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCcgcgc           >  2:4076215/1‑92 (MQ=255)
                                                                                gCCTGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGt           >  1:4346732/1‑93 (MQ=255)
                                                                                  cTGACGATTTTCCCC‑‑GTTCCCCGGTGCTGTACCGGGGACGTATTTTATTTCCCTGCATCGCCCGCCTTCTTTCCGCATCTTCCCCCCGCGTcc         >  2:848660/1‑92 (MQ=255)
                                                                                   tGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTcac        >  2:141400/1‑93 (MQ=255)
                                                                                   tGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTcac        >  2:2785870/1‑93 (MQ=255)
                                                                                   tGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTcac        >  2:640548/1‑93 (MQ=255)
                                                                                   tGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTcac        <  2:3814349/93‑1 (MQ=255)
                                                                                   tGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCCGCGTcac        >  1:1969375/1‑93 (MQ=255)
                                                                                   tGACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCCGCGTcac        >  1:1397459/1‑93 (MQ=255)
                                                                                    gACGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTcaca       <  1:1716957/93‑1 (MQ=255)
                                                                                     aCGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCacac      >  1:3709861/1‑93 (MQ=255)
                                                                                      gattttcCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGGATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCg    >  2:2835143/4‑93 (MQ=255)
                                                                                      cGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACAcc     <  1:1798213/93‑1 (MQ=255)
                                                                                      cGATTTTCCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACAcc     <  2:3214701/93‑1 (MQ=255)
                                                                                       attttcCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCGa   <  2:603828/91‑1 (MQ=255)
                                                                                       attttcCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCGa   <  1:1267806/91‑1 (MQ=255)
                                                                                        aTTTTCCCC‑‑GTTCCCCGTTGCTGTACCGGGGACGTATTTAATTCCCCTGCATCGCCCGCCTTCTTGCCGCATCTTCCCCCCGCGTCCCACCGa   >  2:4142358/1‑93 (MQ=255)
                                                                                              |                                                                                          
TGATGATCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCC‑‑GTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCGAA  >  NC_000913/3152152‑3152332

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