Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,905,091 T→C intergenic (‑7/+159) yebH ← / ← cspC hypothetical protein/cold shock protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,905,0910TC85.4% 574.4 / 84.9 247intergenic (‑7/+159)yebH/cspChypothetical protein/cold shock protein
Reads supporting (aligned to +/- strand):  ref base T (15/21);  new base C (116/95);  total (131/116)
Fisher's exact test for biased strand distribution p-value = 1.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 250 total aligned reads displayed.

TATCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACTCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGTAAAAA  >  NC_000913/1905000‑1905182
                                                                                           |                                                                                           
tatCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACTc                                                                                            >  2:6750971/1‑93 (MQ=255)
 atCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAAcccc                                                                                           <  1:4259210/93‑1 (MQ=255)
 atCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAAcccc                                                                                           >  1:3343935/1‑93 (MQ=255)
 atCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACTcc                                                                                           <  1:1365841/93‑1 (MQ=255)
  tCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCg                                                                                          >  2:1247722/1‑93 (MQ=255)
    aaaCTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACTCCGGc                                                                                        <  1:3824416/93‑1 (MQ=255)
     aaCTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCa                                                                                       >  1:6285792/1‑93 (MQ=255)
      aCTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACTCCGGCAg                                                                                      <  1:6927999/93‑1 (MQ=255)
      aCTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAg                                                                                      <  2:3830118/93‑1 (MQ=255)
        tGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGa                                                                                    >  1:3052668/1‑93 (MQ=255)
         gATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACTCCGGCAGGAt                                                                                   <  2:3445995/93‑1 (MQ=255)
          aTGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATa                                                                                  >  2:4723400/1‑93 (MQ=255)
           tGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATAc                                                                                 <  1:7366071/93‑1 (MQ=255)
             tCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTa                                                                               >  2:3871279/1‑93 (MQ=255)
              cGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTaa                                                                              <  1:1616394/93‑1 (MQ=255)
               gCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACTCCGGCAGGATACTaaa                                                                             >  2:4949705/1‑93 (MQ=255)
                  ggggACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACTCCGGCAGGATACTAAAGGc                                                                          >  2:1345624/1‑93 (MQ=255)
                        aGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAg                                                                    <  1:4011050/93‑1 (MQ=255)
                        aGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAg                                                                    >  1:2274444/1‑93 (MQ=255)
                         gATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAgt                                                                   <  2:2400200/93‑1 (MQ=255)
                         gATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAgt                                                                   >  2:3377756/1‑93 (MQ=255)
                         gATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAgt                                                                   >  1:6476509/1‑93 (MQ=255)
                         gATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAgt                                                                   >  1:6088823/1‑93 (MQ=255)
                           tataGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACTCCGGCAGGATACTAAAGGCCGCAAgtgt                                                                 >  1:7255109/1‑93 (MQ=255)
                           tataGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAgtgt                                                                 >  2:4559825/1‑93 (MQ=255)
                             taGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACTCCGGCAGGATACTAAAGGCCGCAAGTGTaa                                                               <  1:729393/93‑1 (MQ=255)
                             taGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTaa                                                               >  1:4839051/1‑93 (MQ=255)
                             taGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTaa                                                               <  2:3018112/93‑1 (MQ=255)
                             taGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTaa                                                               <  1:2514241/93‑1 (MQ=255)
                             taGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTaa                                                               <  1:1041897/93‑1 (MQ=255)
                               gCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACTCCGGCAGGATACTAAAGGCCGCAAGTGTAAcg                                                             >  2:3510210/1‑93 (MQ=255)
                                cTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAAcgc                                                            >  1:7126516/1‑93 (MQ=255)
                                   tttttttCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGaa                                                         >  2:5112605/1‑93 (MQ=255)
                                   tttttttCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGaa                                                         >  1:1118666/1‑93 (MQ=255)
                                   tttttttCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGaa                                                         <  1:2982572/93‑1 (MQ=255)
                                     tttttCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATc                                                       <  2:3192703/93‑1 (MQ=255)
                                      ttttCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCa                                                      <  2:4136167/93‑1 (MQ=255)
                                       tttCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCAc                                                     <  1:4143492/93‑1 (MQ=255)
                                       tttCACGCGAAAGAAGCTGGTGGCAAAGTGGACAGGAAAAAGACATGAACACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATccc                                                     <  1:2385789/93‑3 (MQ=255)
                                        ttCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGTGAATCAcg                                                    >  2:5349923/1‑93 (MQ=255)
                                         tCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCAcgc                                                   <  2:4670005/93‑1 (MQ=255)
                                          cACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCAcgcg                                                  >  1:4532393/1‑93 (MQ=255)
                                           aCGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACgcgg                                                 >  2:6680125/1‑93 (MQ=255)
                                           aCGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACgcgg                                                 >  2:1791729/1‑93 (MQ=255)
                                            cgcgAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACTCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACgcggc                                                <  2:5976995/93‑1 (MQ=255)
                                            cgcgAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACgcggc                                                >  1:2718819/1‑93 (MQ=255)
                                             gcgAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACgcggcg                                               <  1:4718683/93‑1 (MQ=255)
                                              cgAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGcggcgg                                              <  1:6628584/93‑1 (MQ=255)
                                               gAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACTCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCggg                                             >  1:5949988/1‑93 (MQ=255)
                                               gAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACTCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCggg                                             <  2:6166574/93‑1 (MQ=255)
                                               gAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCggg                                             <  2:7519927/93‑1 (MQ=255)
                                                aaaGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACTCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGc                                            >  1:6859964/1‑93 (MQ=255)
                                                aaaGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGc                                            >  2:5597031/1‑93 (MQ=255)
                                                aaaGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGc                                            >  1:5664413/1‑93 (MQ=255)
                                                  agagGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCaa                                          >  1:4985382/1‑93 (MQ=255)
                                                  agagGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCaa                                          >  2:2072161/1‑93 (MQ=255)
                                                  agagGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCaa                                          >  1:3393282/1‑93 (MQ=255)
                                                   gagGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACTCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCaaa                                         <  2:2627828/93‑1 (MQ=255)
                                                    agGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAat                                        >  2:4299482/1‑93 (MQ=255)
                                                     ggCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACTCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAata                                       <  1:3485985/93‑1 (MQ=255)
                                                     ggCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACTCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAata                                       <  2:2264782/93‑1 (MQ=255)
                                                       cTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATc                                     <  1:3110908/93‑1 (MQ=255)
                                                       cTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATc                                     >  2:84881/1‑93 (MQ=255)
                                                         gATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCg                                   >  2:1714422/1‑93 (MQ=255)
                                                         gATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCg                                   >  1:148240/1‑93 (MQ=255)
                                                          aTGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGa                                  >  1:261276/1‑93 (MQ=255)
                                                          aTGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGa                                  <  2:3565173/93‑1 (MQ=255)
                                                             gCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACg                               >  2:6975803/1‑93 (MQ=255)
                                                             gCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACg                               >  2:1558852/1‑93 (MQ=255)
                                                              cAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGc                              >  2:2809231/1‑93 (MQ=255)
                                                                 aGTGGACAGGAAAAAGACATGACAACTCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAg                           <  2:6318600/93‑1 (MQ=255)
                                                                 aGTGGACAGGAAAAAGACATGACAACTCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAg                           <  2:5392111/93‑1 (MQ=255)
                                                                 aGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAg                           <  1:5212019/93‑1 (MQ=255)
                                                                 aGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAg                           <  1:6257872/93‑1 (MQ=255)
                                                                     gACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGc                       <  2:3106658/93‑1 (MQ=255)
                                                                      aCAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCg                      <  1:3872626/93‑1 (MQ=255)
                                                                        aNGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGaa                    <  2:7202813/93‑1 (MQ=255)
                                                                        aGGAAAAAGACATGACAACTCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGaa                    <  2:1576290/93‑1 (MQ=255)
                                                                        aGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGaa                    <  2:1947880/93‑1 (MQ=255)
                                                                        aGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGaa                    >  2:4062941/1‑93 (MQ=255)
                                                                        aGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGaa                    <  2:2712593/93‑1 (MQ=255)
                                                                        aGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGaa                    >  2:5691814/1‑93 (MQ=255)
                                                                         ggAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAg                   >  1:128625/1‑93 (MQ=255)
                                                                          gAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGa                  <  1:912199/93‑1 (MQ=255)
                                                                          gAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGa                  >  2:4202078/1‑93 (MQ=255)
                                                                           aaaaaGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGAt                 >  1:6995587/1‑93 (MQ=255)
                                                                           aaaaaGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGAt                 >  1:4054444/1‑93 (MQ=255)
                                                                            aaaaGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATa                >  1:3559845/1‑93 (MQ=255)
                                                                              aaGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATaaa              >  1:1870383/1‑93 (MQ=255)
                                                                                  cATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCAt          <  1:5560937/93‑1 (MQ=255)
                                                                                  cATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCAt          <  2:4469700/93‑1 (MQ=255)
                                                                                   aTGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATa         >  1:2733564/1‑93 (MQ=255)
                                                                                   aTGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATa         <  1:2030544/93‑1 (MQ=255)
                                                                                    tGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAg        <  1:1878709/93‑1 (MQ=255)
                                                                                    tGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAg        >  2:6881337/1‑93 (MQ=255)
                                                                                    tGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAg        <  1:5382907/93‑1 (MQ=255)
                                                                                    tGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAg        >  1:2586605/1‑93 (MQ=255)
                                                                                     gACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGt       >  1:4582054/1‑93 (MQ=255)
                                                                                       cAACTCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGTaa     <  1:5753843/93‑1 (MQ=255)
                                                                                       cAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGTaa     >  1:3221669/1‑93 (MQ=255)
                                                                                           |                                                                                           
TATCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACTCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGTAAAAA  >  NC_000913/1905000‑1905182

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