Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,917,0050TC72.7% 348.1 / 176.2 242M493T (ATG→ACG) b1834hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (35/31);  new base C (96/80);  total (131/111)
Fisher's exact test for biased strand distribution p-value = 8.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 243 total aligned reads displayed.

GGTGGGGCTGGATGGCGTTGAGTTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTTAGCGATTTACCCACCACAACCGTGAGTTTGAGTGCA  >  NC_000913/1916913‑1917096
                                                                                            |                                                                                           
ggTGGGGCTGGATGGCGTTGAGTTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGAt                                                                                             >  2:405478/1‑93 (MQ=255)
    gggCTGGATGGCGTTGAGTTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGGTACGTATTCTGCCGGGCGATAAAGGGGAGATGaaa                                                                                         >  1:61362/1‑93 (MQ=255)
    gggCTGGATGGCGTTGAGTTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGaaa                                                                                         >  2:6740746/1‑93 (MQ=255)
    gggCTGGATGGCGTTGAGTTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGACGaaa                                                                                         <  1:528953/93‑1 (MQ=255)
        tGGATGGCGTTGAGTTTCTCGGTGCTAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCa                                                                                     >  2:5493437/1‑93 (MQ=255)
        tGGATGGCGTTGAGTTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCa                                                                                     >  2:2167130/1‑93 (MQ=255)
        tGGATGGCGTTGAGTTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCa                                                                                     >  1:664800/1‑93 (MQ=255)
        tGGATGGCGTTGAGTTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCa                                                                                     <  1:399758/93‑1 (MQ=255)
        tGGATGGCGTTGAGTTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCa                                                                                     <  1:3032487/93‑1 (MQ=255)
         ggATGGCGTTGAGTTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAg                                                                                    >  2:3074488/1‑93 (MQ=255)
          gATGGCGTTGAGTTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGc                                                                                   <  2:518773/93‑1 (MQ=255)
          gATGGCGTTGAGTTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGc                                                                                   <  1:3634254/93‑1 (MQ=255)
          gATGGCGTTGAGTTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGc                                                                                   >  1:5225071/1‑93 (MQ=255)
          gATGGCGTTGAGTTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGc                                                                                   <  1:3540918/93‑1 (MQ=255)
           aTGGCGTTGAGTTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCt                                                                                  <  2:3305562/93‑1 (MQ=255)
            tGGCGTTGAGTTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTa                                                                                 <  2:4307389/93‑1 (MQ=255)
            tGGCGTTGAGTTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTa                                                                                 <  1:7265967/93‑1 (MQ=255)
              gCGTTGAGTTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATc                                                                               <  2:6108887/93‑1 (MQ=255)
                gTTGAGTTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCa                                                                             <  2:3666975/93‑1 (MQ=255)
                gTTGAGTTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCa                                                                             >  2:3950122/1‑93 (MQ=255)
                  tGAGTTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACt                                                                           >  2:4607287/1‑93 (MQ=255)
                   gAGTTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg                                                                          >  1:4728550/1‑93 (MQ=255)
                    aGTTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGt                                                                         <  2:2342071/93‑1 (MQ=255)
                     gTTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTa                                                                        >  1:2637067/1‑93 (MQ=255)
                      tttCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTAt                                                                       <  1:6103075/93‑1 (MQ=255)
                       ttCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATg                                                                      <  1:6726710/93‑1 (MQ=255)
                        tctcGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGc                                                                     >  2:5150678/1‑93 (MQ=255)
                         ctcGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGcc                                                                    <  2:3976862/93‑1 (MQ=255)
                          tcGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCa                                                                   >  2:2970376/1‑93 (MQ=255)
                             gTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCCGCTATCCACTGTATGCCCAtt                                                                >  1:1295046/1‑93 (MQ=255)
                               gCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTg                                                              >  1:781349/1‑93 (MQ=255)
                               gCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTg                                                              >  2:6711077/1‑93 (MQ=255)
                                ccAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTgg                                                             <  1:2525805/93‑1 (MQ=255)
                                ccAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTgg                                                             <  2:3322991/93‑1 (MQ=255)
                                 cAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGa                                                            <  2:6950784/93‑1 (MQ=255)
                                  aGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGaa                                                           >  1:734463/1‑93 (MQ=255)
                                   gcgcCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGaaa                                                          <  2:2478617/93‑1 (MQ=255)
                                   gcgcCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGaaa                                                          >  1:2042594/1‑93 (MQ=255)
                                       cTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAgc                                                      >  2:4224587/1‑93 (MQ=255)
                                       cTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAgc                                                      >  2:3140400/1‑93 (MQ=255)
                                       cTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAgc                                                      >  2:4664369/1‑93 (MQ=255)
                                        tCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAgcg                                                     <  1:6411307/93‑1 (MQ=255)
                                        tCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAgcg                                                     >  2:720389/1‑93 (MQ=255)
                                         cAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAgcgc                                                    <  2:91956/93‑1 (MQ=255)
                                         cAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAgcgc                                                    <  1:3830847/93‑1 (MQ=255)
                                          aGAATGGATCAATGGCGGGATACGTATTCTGCCGGGAGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCt                                                   >  1:7472219/1‑93 (MQ=255)
                                           gAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTg                                                  >  2:6985460/1‑93 (MQ=255)
                                            aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTgg                                                 >  2:5629975/1‑93 (MQ=255)
                                             aTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGga                                                >  1:5535204/1‑93 (MQ=255)
                                             aTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGga                                                <  2:1535581/93‑1 (MQ=255)
                                             aTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGGGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGga                                                >  1:6174205/1‑93 (MQ=255)
                                              tGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGgag                                               >  1:5746145/1‑93 (MQ=255)
                                               ggATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGgaga                                              >  2:6494866/1‑93 (MQ=255)
                                                 aTCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAAc                                            <  1:2294337/93‑1 (MQ=255)
                                                  tCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACa                                           <  1:1338990/93‑1 (MQ=255)
                                                  tCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACa                                           <  2:3825546/93‑1 (MQ=255)
                                                  tCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACa                                           >  1:2432019/1‑93 (MQ=255)
                                                  tCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACa                                           >  2:2826226/1‑93 (MQ=255)
                                                   cAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGGGAGATGAAAGCCATCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAg                                          >  1:7373567/1‑93 (MQ=255)
                                                    aaTGGCGGGATACGTATTCTGCCGGGCGATAAAGGGGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGc                                         >  1:6869670/1‑93 (MQ=255)
                                                     aTGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGcc                                        <  2:322119/93‑1 (MQ=255)
                                                     aTGGCGGGATAAGTATTCTTCCGGGCGATAAAGGGGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTTGAAAAAGCGCTGGAGAACAGcc                                        >  1:4201498/1‑93 (MQ=255)
                                                      tGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCt                                       <  2:4974765/93‑1 (MQ=255)
                                                      tGGCGGGATACGTATTCTGCCGGGCGATAAAGGGGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCt                                       >  2:1458606/1‑93 (MQ=255)
                                                       ggCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCtt                                      <  1:5894208/93‑1 (MQ=255)
                                                        gCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTTa                                     <  2:7067386/93‑1 (MQ=255)
                                                          gggATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTTAGc                                   <  1:5418488/93‑1 (MQ=255)
                                                          gggATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTTAGc                                   <  1:1464639/93‑1 (MQ=255)
                                                          gggATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTTAGc                                   >  1:4163753/1‑93 (MQ=255)
                                                          gggATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTTAGc                                   >  2:4345109/1‑93 (MQ=255)
                                                           ggATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTTAGCg                                  >  2:1076955/1‑93 (MQ=255)
                                                            gATACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTTAGCGa                                 >  1:1709047/1‑93 (MQ=255)
                                                              tACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTTAGCGAtt                               <  1:210090/93‑1 (MQ=255)
                                                              tACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTTAGCGAtt                               >  1:2589305/1‑93 (MQ=255)
                                                              tACGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTTAGCGAtt                               >  2:487993/1‑93 (MQ=255)
                                                               aCGTATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTTAGCGAttt                              >  1:1382208/1‑93 (MQ=255)
                                                                  tATTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTTAGCGATTTAcc                           >  2:4155457/1‑93 (MQ=255)
                                                                   aTTCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTTAGCGATTTAccc                          <  1:6291606/93‑1 (MQ=255)
                                                                    ttCTGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTTAGCGATTTACcca                         <  2:4589618/93‑1 (MQ=255)
                                                                      cTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTTAGCGATTTACccacc                       <  1:2245000/93‑1 (MQ=255)
                                                                       tGCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTTAGCGATTTACccacca                      >  2:1735727/1‑93 (MQ=255)
                                                                        gCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTTAGCGATTTACCcaccac                     <  2:6370002/93‑1 (MQ=255)
                                                                        gCCGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTTAGCGATTTACCcaccac                     <  2:2614715/93‑1 (MQ=255)
                                                                          cGGGCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTTAGCGATTTACCCACCACaa                   >  2:6280753/1‑93 (MQ=255)
                                                                           gggCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTTAGCGATTTACCCACCACAAc                  >  2:6629649/1‑93 (MQ=255)
                                                                           gggCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTTAGCGATTTACCCACCACAAc                  <  1:973511/93‑1 (MQ=255)
                                                                           gggCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTTAGCGATTTACCCACCACAAc                  >  2:1935276/1‑93 (MQ=255)
                                                                             gCGATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTTAGCGATTTACCCACCACAACCg                >  2:1821350/1‑93 (MQ=255)
                                                                               gATAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTTAGCGATTTACCCACCACAACCGTg              >  1:5695312/1‑93 (MQ=255)
                                                                                aTAAAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTTAGCGATTTACCCACCACAACCGTGa             >  1:407179/1‑93 (MQ=255)
                                                                                aTAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTTAGCGATTTACCCACCACAACCGTGa             >  1:1902739/1‑93 (MQ=255)
                                                                                aTAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTTAGCGATTTACCCACCACAACCGTGa             <  2:5838507/93‑1 (MQ=255)
                                                                                aTAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTTAGCGATTTACCCACCACAACCGTGa             >  1:1145157/1‑93 (MQ=255)
                                                                                 tAAAGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTTAGCGATTTACCCACCACAACCGTGAg            >  1:3195138/1‑93 (MQ=255)
                                                                                    aGGTGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTTAGCGATTTACCCACCACAACCGTGAGttt         >  2:1474356/1‑93 (MQ=255)
                                                                                       tGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTTAGCGATTTACCCACCACAACCGTGAGTTTGAg      >  2:1292641/1‑93 (MQ=255)
                                                                                       tGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTTAGCGATTTACCCACCACAACCGTGAGTTTGAg      <  1:704560/93‑1 (MQ=255)
                                                                                       tGAGACGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTTAGCGATTTACCCACCACAACCGTGAGTTTGAg      <  2:2718256/93‑1 (MQ=255)
                                                                                        gagaCGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTTAGCGATTTACCCACCACAACCGTGAGTTTGAGt     >  1:4438014/1‑93 (MQ=255)
                                                                                        gagaCGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTTAGCGATTTACCCACCACAACCGTGAGTTTGAGt     >  2:7511778/1‑93 (MQ=255)
                                                                                            |                                                                                           
GGTGGGGCTGGATGGCGTTGAGTTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTTAGCGATTTACCCACCACAACCGTGAGTTTGAGTGCA  >  NC_000913/1916913‑1917096

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