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,0910TC99.6% 781.2 / ‑5.6 230intergenic (‑7/+159)yebH/cspChypothetical protein/cold shock protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (121/107);  total (123/107)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 234 total aligned reads displayed.

ATCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACTCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGTAAA  >  NC_000913/1905001‑1905180
                                                                                          |                                                                                         
atCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAAcccc                                                                                         <  2:4577686/93‑1 (MQ=255)
 tCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCg                                                                                        >  1:6193157/1‑93 (MQ=255)
 tCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCg                                                                                        <  1:6359359/93‑1 (MQ=255)
  cAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCgg                                                                                       <  1:364637/93‑1 (MQ=255)
    aaCTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCa                                                                                     >  2:4036912/1‑93 (MQ=255)
    aaCTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCa                                                                                     >  1:1413287/1‑93 (MQ=255)
      cTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAgg                                                                                   <  2:6160683/93‑1 (MQ=255)
        gATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGAt                                                                                 >  2:598891/1‑93 (MQ=255)
        gATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGAt                                                                                 >  2:1452036/1‑93 (MQ=255)
        gATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGAt                                                                                 <  2:1168578/93‑1 (MQ=255)
          tGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCACGGCAGGATAc                                                                               >  1:3544993/1‑93 (MQ=255)
           gTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACt                                                                              >  1:1351261/1‑93 (MQ=255)
               cTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAg                                                                          <  2:5947854/93‑1 (MQ=255)
               cTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAg                                                                          <  1:5613975/93‑1 (MQ=255)
               atGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAg                                                                          <  1:1482122/92‑1 (MQ=255)
                tGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACATCCCCGGCAGGATACTAAAgg                                                                         <  1:4285309/93‑1 (MQ=255)
                tGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAgg                                                                         >  1:249766/1‑93 (MQ=255)
                tGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAgg                                                                         <  1:2653236/93‑1 (MQ=255)
                tGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAgg                                                                         >  1:545238/1‑93 (MQ=255)
                gggggACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAgg                                                                         <  1:3355417/92‑1 (MQ=255)
                  gggACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGcc                                                                       >  2:1718801/1‑93 (MQ=255)
                  gggACAGATATAGCTGGTTTTTTCACGCGAAAGAGGCTGATGGTAAAGTGTAAAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGcc                                                                       <  1:6425174/93‑1 (MQ=255)
                        gATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAgt                                                                 >  2:6305831/1‑93 (MQ=255)
                        gATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAgt                                                                 >  2:732391/1‑93 (MQ=255)
                        gATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAgg                                                                 >  2:3472494/1‑92 (MQ=255)
                          tataGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAgtgt                                                               >  1:1058942/1‑93 (MQ=255)
                          tataGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAgtga                                                               >  2:6475781/1‑92 (MQ=255)
                           ataGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTa                                                              >  2:5427027/1‑93 (MQ=255)
                           ataGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTa                                                              <  2:2204867/93‑1 (MQ=255)
                            taGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTaa                                                             <  2:4967018/93‑1 (MQ=255)
                            taGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTaa                                                             <  1:3040012/93‑1 (MQ=255)
                             aGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAAc                                                            <  2:3965926/93‑1 (MQ=255)
                              gCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAAcg                                                           >  2:1952116/1‑93 (MQ=255)
                               cTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAAcgc                                                          >  1:3408791/1‑93 (MQ=255)
                               cTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAAcgc                                                          <  1:2969663/93‑1 (MQ=255)
                                tGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAAcgcg                                                         >  1:2553225/1‑93 (MQ=255)
                                 gTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGa                                                        >  1:3680224/1‑93 (MQ=255)
                                 gTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGa                                                        >  1:1135575/1‑93 (MQ=255)
                                  tttttttCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGaa                                                       <  2:4470723/93‑1 (MQ=255)
                                  tgtttttGACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGaa                                                       <  1:4312679/91‑1 (MQ=255)
                                      tttCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCAc                                                   <  1:1024655/93‑1 (MQ=255)
                                         catgcgAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCAcgcg                                                >  1:3940858/4‑93 (MQ=255)
                                          aCGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACgcgg                                               >  2:1080865/1‑93 (MQ=255)
                                            gcgAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACgcggcg                                             <  1:4638908/93‑1 (MQ=255)
                                             cgAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGcggcgg                                            <  1:71844/93‑1 (MQ=255)
                                             cgAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGcggcgg                                            <  2:3511198/93‑1 (MQ=255)
                                               aaaGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGc                                          >  2:4159246/1‑93 (MQ=255)
                                                aaGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCa                                         >  2:5265903/1‑93 (MQ=255)
                                                aaGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCa                                         >  1:1554795/1‑93 (MQ=255)
                                                aaGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCa                                         >  1:3740201/1‑93 (MQ=255)
                                                aaGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCa                                         <  1:2055279/93‑1 (MQ=255)
                                                aaGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACACCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCa                                         >  2:4534178/1‑93 (MQ=255)
                                                 agagGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCaa                                        >  2:4603064/1‑93 (MQ=255)
                                                 agagGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCaa                                        >  2:623357/1‑93 (MQ=255)
                                                 agagGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATAACGCGGCGGGCaa                                        >  1:3005759/1‑93 (MQ=255)
                                                    ggCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAata                                     <  1:501561/93‑1 (MQ=255)
                                                    ggCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAata                                     <  1:5019852/93‑1 (MQ=255)
                                                     gCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAatat                                    <  2:1975057/93‑1 (MQ=255)
                                                     gCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAatat                                    <  1:1745091/93‑1 (MQ=255)
                                                     gCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAatat                                    >  1:3117934/1‑93 (MQ=255)
                                                      cTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATc                                   <  1:4930084/93‑1 (MQ=255)
                                                      cTGATGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATc                                   <  2:3200942/93‑1 (MQ=255)
                                                         aTGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGa                                >  1:2200751/1‑93 (MQ=255)
                                                         aTGGCAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGa                                >  1:1942949/1‑93 (MQ=255)
                                                             cAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAATTATCCGAACGc                            >  2:5546768/1‑93 (MQ=255)
                                                             cAAAGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGc                            <  2:2514698/93‑1 (MQ=255)
                                                              aaaGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGcc                           >  1:3103030/1‑93 (MQ=255)
                                                               aaGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCa                          >  2:6013095/1‑93 (MQ=255)
                                                               aaGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCa                          >  1:3845229/1‑93 (MQ=255)
                                                                aGTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAg                         <  1:5530714/93‑1 (MQ=255)
                                                                 gTGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAgg                        >  2:244459/1‑93 (MQ=255)
                                                                  tGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGa                       <  2:5049588/93‑1 (MQ=255)
                                                                  tGGACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGa                       <  2:5670902/93‑1 (MQ=255)
                                                                   ggACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAg                      <  2:182795/93‑1 (MQ=255)
                                                                   ggACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAg                      >  1:3055060/1‑93 (MQ=255)
                                                                    gACAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGc                     >  1:6038456/1‑93 (MQ=255)
                                                                     aCAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCg                    >  2:215356/1‑93 (MQ=255)
                                                                     aCAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCg                    >  1:6097077/1‑93 (MQ=255)
                                                                      cAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGa                   <  2:1304302/93‑1 (MQ=255)
                                                                      cAGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGa                   >  2:4798320/1‑93 (MQ=255)
                                                                       aGGAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGaa                  >  2:3656019/1‑93 (MQ=255)
                                                                        ggAAAAAGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAg                 >  1:4654315/1‑93 (MQ=255)
                                                                          aaaaaGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGAt               >  2:3951976/1‑93 (MQ=255)
                                                                          aaaaaGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGAt               <  1:774737/93‑1 (MQ=255)
                                                                           aaaaGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATa              <  2:2747741/93‑1 (MQ=255)
                                                                           aaaaGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATa              >  2:2848588/1‑93 (MQ=255)
                                                                             aaGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATaaa            <  1:4521532/93‑1 (MQ=255)
                                                                             aaGACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATaaa            <  2:4750610/93‑1 (MQ=255)
                                                                               gACATGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGc          >  1:5806987/1‑93 (MQ=255)
                                                                                  aTGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATa       >  2:1399189/1‑93 (MQ=255)
                                                                                   tGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAg      >  1:4003454/1‑93 (MQ=255)
                                                                                   tGACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAg      >  1:32984/1‑93 (MQ=255)
                                                                                    gACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGt     >  2:3220993/1‑93 (MQ=255)
                                                                                    gACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGt     >  2:3094002/1‑93 (MQ=255)
                                                                                    gACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGt     >  2:1478627/1‑93 (MQ=255)
                                                                                    gACAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCAGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGt     >  1:815730/1‑93 (MQ=255)
                                                                                     aCAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGTa    >  2:5314682/1‑93 (MQ=255)
                                                                                     aCAACCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGTa    <  1:965136/93‑1 (MQ=255)
                                                                                       aaCCCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGTaaa  >  1:1886233/1‑93 (MQ=255)
                                                                                       aaCCCCGGCAGGATACTAAAGGCCGCAAGTGCAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGTAGCGAAGATAAAGCATAGTaaa  >  2:5639641/1‑93 (MQ=255)
                                                                                          |                                                                                         
ATCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACTCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGTAAA  >  NC_000913/1905001‑1905180

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