Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,636,465 T→C intergenic (‑296/+14) b1551 ← / ← cspI hypothetical protein/cold shock‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,636,4650TC99.4% 581.9 / ‑4.4 167intergenic (‑296/+14)b1551/cspIhypothetical protein/cold shock‑like protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (90/76);  total (90/77)
Fisher's exact test for biased strand distribution p-value = 4.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

Only 100 of 169 total aligned reads displayed.

CCCCGAAAGGATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGTAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCA  >  NC_000913/1636375‑1636554
                                                                                          |                                                                                         
ccccGAAAGGATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCaa                                                                                         >  1:1803757/1‑93 (MQ=255)
ccccGAAAGGATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCaa                                                                                         <  2:5084230/93‑1 (MQ=255)
ccccGAAAGGATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCaa                                                                                         >  2:1240384/1‑93 (MQ=255)
 cccGAAAGGATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAat                                                                                        >  1:3876678/1‑93 (MQ=255)
 cccGAAAGGATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAat                                                                                        >  1:3709377/1‑93 (MQ=255)
  ccGAAAGGATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAata                                                                                       >  2:6110537/1‑93 (MQ=255)
   cGAAAGGATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAatat                                                                                      >  2:1676529/1‑93 (MQ=255)
   cGAAAGGATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAatat                                                                                      >  2:4418947/1‑93 (MQ=255)
     aaaGGATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTg                                                                                    >  1:3374768/1‑93 (MQ=255)
     aaaGGATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTg                                                                                    >  1:4116316/1‑93 (MQ=255)
     aaaGGATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTg                                                                                    >  1:1345134/1‑93 (MQ=255)
       aGGATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTc                                                                                  <  2:5749326/93‑1 (MQ=255)
        ggATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCt                                                                                 <  2:6236194/93‑1 (MQ=255)
        ggATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCt                                                                                 <  1:5302112/93‑1 (MQ=255)
        ggATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCt                                                                                 <  1:6181011/93‑1 (MQ=255)
        ggATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCt                                                                                 <  2:103652/93‑1 (MQ=255)
         gATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTa                                                                                >  1:4385338/1‑93 (MQ=255)
          aTTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTAc                                                                               >  1:3766733/1‑93 (MQ=255)
           ttCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTAcc                                                                              <  2:3764899/93‑1 (MQ=255)
            tCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCt                                                                             >  1:4062583/1‑93 (MQ=255)
            tCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCt                                                                             >  2:1631966/1‑93 (MQ=255)
             cATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTc                                                                            >  1:2687859/1‑93 (MQ=255)
                gTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCaaa                                                                         <  1:5410952/93‑1 (MQ=255)
                gTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCaaa                                                                         >  1:5474701/1‑93 (MQ=255)
                 tCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAg                                                                        <  2:478352/93‑1 (MQ=255)
                   aGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAgcg                                                                      >  2:682546/1‑93 (MQ=255)
                    gAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAgcgc                                                                     >  2:4134611/1‑93 (MQ=255)
                    gAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAgcgc                                                                     >  1:2815441/1‑93 (MQ=255)
                       aCATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCAc                                                                  >  1:2037049/1‑93 (MQ=255)
                       aCATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCAc                                                                  >  2:31966/1‑93 (MQ=255)
                       aCATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCAc                                                                  <  2:3297382/93‑1 (MQ=255)
                       aCATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCAc                                                                  >  1:3354649/1‑93 (MQ=255)
                       aCATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCAc                                                                  <  1:3048536/93‑1 (MQ=255)
                        cATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACt                                                                 <  2:1267566/93‑1 (MQ=255)
                         aTCAAAAGATGATGTTCTATACTGGGAACCATGATAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTa                                                                >  2:1077684/1‑93 (MQ=255)
                         aTCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTa                                                                >  1:4208301/1‑93 (MQ=255)
                           cAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACa                                                              <  1:2482790/93‑1 (MQ=255)
                            aaaaGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTGCCTCAAAGCGACACTACAt                                                             <  1:6319786/93‑1 (MQ=255)
                             aaaGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATg                                                            >  2:3059402/1‑93 (MQ=255)
                              aaGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGa                                                           >  1:4287910/1‑93 (MQ=255)
                              aaGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGa                                                           >  1:5290195/1‑93 (MQ=255)
                              aaGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGa                                                           >  1:2463745/1‑93 (MQ=255)
                               aGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGaa                                                          <  1:2325631/93‑1 (MQ=255)
                               aGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGaa                                                          >  1:6459377/1‑93 (MQ=255)
                               aGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGaa                                                          >  2:5348271/1‑93 (MQ=255)
                                gatgatGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGaaa                                                         >  1:100043/1‑92 (MQ=255)
                                gatgatGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAAc                                                         >  1:595796/1‑93 (MQ=255)
                                 atgatgTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACa                                                        <  1:4810781/93‑1 (MQ=255)
                                   gatgTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAgg                                                      >  2:2042491/1‑92 (MQ=255)
                                       ttCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCa                                                  >  2:5457526/1‑93 (MQ=255)
                                        tCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAg                                                 <  1:5000839/93‑1 (MQ=255)
                                        tCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAg                                                 <  2:2654568/93‑1 (MQ=255)
                                          tataCTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGa                                               >  1:2685660/1‑93 (MQ=255)
                                          tataCTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGa                                               >  1:3560500/1‑93 (MQ=255)
                                          tataCTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGa                                               >  2:4437165/1‑93 (MQ=255)
                                          tataCTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGa                                               >  1:485767/1‑93 (MQ=255)
                                            taCTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGAcc                                             >  2:893544/1‑93 (MQ=255)
                                               tGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCttt                                          >  2:6359328/1‑93 (MQ=255)
                                                gggAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTa                                         >  1:696090/1‑93 (MQ=255)
                                                 ggAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAg                                        <  1:871910/93‑1 (MQ=255)
                                                 ggAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAg                                        >  1:1716775/1‑93 (MQ=255)
                                                 ggAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAg                                        >  1:3117707/1‑93 (MQ=255)
                                                 ggAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAg                                        >  1:4502600/1‑93 (MQ=255)
                                                  gAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCCTTAgg                                       >  1:3886881/1‑93 (MQ=255)
                                                   aaCCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGt                                      <  2:4797264/93‑1 (MQ=255)
                                                   aaCCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGt                                      <  2:413090/93‑1 (MQ=255)
                                                   aaCCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGt                                      <  2:3720/93‑1 (MQ=255)
                                                    aCCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAGCAGCGGCAGGACCTTTAGGTc                                     <  2:1189252/93‑1 (MQ=255)
                                                     ccATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGGCCTTTAAGTcc                                    >  1:6389717/1‑93 (MQ=255)
                                                      cATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGAGCTTTAGGTCCg                                   >  2:3871710/1‑93 (MQ=255)
                                                       aTGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGt                                  <  2:1594187/93‑1 (MQ=255)
                                                          acaacaCGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCCTTAGGTCCGGtct                               >  2:393843/1‑93 (MQ=255)
                                                             acacGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCaa                            >  1:3220726/1‑93 (MQ=255)
                                                               acGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATa                          >  2:2466969/1‑93 (MQ=255)
                                                                 gggCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCCATAcc                        >  1:2547021/1‑93 (MQ=255)
                                                                 gggCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATAcc                        >  1:1050764/1‑93 (MQ=255)
                                                                  ggCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCa                       >  1:5527000/1‑93 (MQ=255)
                                                                   gCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCaa                      <  2:295220/93‑1 (MQ=255)
                                                                     aTCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAAt                    >  1:1632931/1‑93 (MQ=255)
                                                                       cTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTc                  <  2:1999403/93‑1 (MQ=255)
                                                                       cTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCCTTAGGTCCGGTCTCCATACCAAATTc                  >  1:3319002/1‑93 (MQ=255)
                                                                         aaaGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCaa                <  2:3998793/93‑1 (MQ=255)
                                                                          aaGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAAc               <  2:1886082/93‑1 (MQ=255)
                                                                           aGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACt              <  2:4417408/93‑1 (MQ=255)
                                                                           aGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACt              >  2:4705687/1‑93 (MQ=255)
                                                                            gTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACtt             <  2:2305135/93‑1 (MQ=255)
                                                                            gTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACtt             <  2:1469125/93‑1 (MQ=255)
                                                                             tGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTc            <  2:389709/93‑1 (MQ=255)
                                                                             tGAATATGCTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTc            <  1:2383822/93‑1 (MQ=255)
                                                                               aaTATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCt          <  2:2787035/93‑1 (MQ=255)
                                                                               aaTATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCt          <  1:3615113/93‑1 (MQ=255)
                                                                               aaTATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCt          <  2:4300101/93‑1 (MQ=255)
                                                                               aaTATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCCTTAGGTCCGTTCTCAATAACAAATTCAACTTCCt          >  2:6199693/1‑93 (MQ=255)
                                                                                atatGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTg         >  1:4701899/1‑93 (MQ=255)
                                                                                 tatGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTgt        >  2:5329317/1‑92 (MQ=255)
                                                                                    agTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATtc     <  2:2572004/92‑1 (MQ=255)
                                                                                     gTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGGCCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATtct    >  2:3492541/1‑93 (MQ=255)
                                                                                      tttGCAATATTGTCTACCTCAAAGCGCCACTACATGACCACCGGCAGGACCTTTAGGCCCGTTTTCAATACCAAATTCAACTTCCTGATtccc   >  2:2880225/1‑91 (MQ=255)
                                                                                      tttGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATtctc   <  1:4831857/93‑1 (MQ=255)
                                                                                       ttGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCa  <  1:6365598/93‑1 (MQ=255)
                                                                                          |                                                                                         
CCCCGAAAGGATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGTAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCA  >  NC_000913/1636375‑1636554

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