Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,446,3620A.25.0% 341.3 / 119.0 132intergenic (+156/+16)ydbL/feaRDUF1318 family protein/transcriptional activator for tynA and feaB
Reads supporting (aligned to +/- strand):  ref base A (45/54);  new base . (20/13);  total (65/67)
Fisher's exact test for biased strand distribution p-value = 1.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.95e-01

Only 100 of 137 total aligned reads displayed.

CTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAA  >  NC_000913/1446263‑1446461
                                                                                                     |                                                                                                   
cTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACa                                                                                                     >  1:806273/1‑100 (MQ=255)
cTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACa                                                                                                     >  2:363624/1‑100 (MQ=255)
cTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACa                                                                                                     >  2:52563/1‑100 (MQ=255)
 tGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACAt                                                                                                    <  1:1814329/100‑1 (MQ=255)
  gAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATg                                                                                                   <  1:2385598/100‑1 (MQ=255)
   aaaCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGc                                                                                                  <  1:2101150/100‑1 (MQ=255)
   aaaCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGc                                                                                                  <  1:3784030/100‑1 (MQ=255)
        aGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTTCAAAATCCATACTCAAATAGCTGTTTTATCTAATAGACGATGC‑A‑Gtacataca                                                                                                 >  2:767388/1‑96 (MQ=255)
        aGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑ATGTAC‑TGCAACGc                                                                                             <  2:951158/100‑1 (MQ=38)
        aGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑ATGTAC‑TGCAACGc                                                                                             <  1:941701/100‑1 (MQ=38)
        aGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑ATGTAC‑TGCAACGc                                                                                             <  1:74787/100‑1 (MQ=38)
         ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑ATGTAC‑TGCAACGcc                                                                                            <  2:1474368/100‑1 (MQ=38)
         ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑ATGTAC‑TGCAACGcc                                                                                            <  1:1539872/100‑1 (MQ=38)
         ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGcc                                                                                            >  1:782502/1‑100 (MQ=255)
         ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGcc                                                                                            <  2:229918/100‑1 (MQ=255)
          gTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑ATGTAC‑TGCAACGCCa                                                                                           <  1:2224621/100‑1 (MQ=38)
          gTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCa                                                                                           <  1:350649/100‑1 (MQ=255)
            gAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑ATGTAC‑TGCAACGCCAAt                                                                                         <  1:1861483/100‑1 (MQ=38)
             aaTATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑ATGTAC‑TGCAACGCCAAtt                                                                                        <  2:3250433/100‑1 (MQ=38)
             aaTATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAAtt                                                                                        <  1:2334426/100‑1 (MQ=255)
             aaTATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAAtt                                                                                        <  1:2565121/100‑1 (MQ=255)
             aaTATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAAtt                                                                                        <  2:806271/100‑1 (MQ=255)
              atatTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAAtta                                                                                       <  2:3549618/100‑1 (MQ=255)
               natTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAAttag                                                                                      >  1:352057/2‑100 (MQ=255)
                 ttCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAAttagtt                                                                                    <  2:3393799/100‑1 (MQ=255)
                 ttCAAAAACTCCTGTCAAATTGCATTATGAACAGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAAttagtt                                                                                    <  2:3528955/100‑1 (MQ=38)
                  tCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAAttagtta                                                                                   >  1:3411919/1‑100 (MQ=255)
                   cAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAAttagttag                                                                                  >  1:2782880/1‑100 (MQ=255)
                    aaaaaCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAAAGCCAAttagtta                                                                                   >  2:1707637/1‑98 (MQ=255)
                     aaaaCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCg                                                                                <  1:2108639/100‑1 (MQ=255)
                     aaaaCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCg                                                                                >  1:2980395/1‑100 (MQ=255)
                      aaaCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCgg                                                                               <  2:1984965/100‑1 (MQ=255)
                            cTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTa                                                                         >  1:1478104/1‑100 (MQ=255)
                             tGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATACAACGCCAATTAGTTAGCGGAATTTAc                                                                        >  2:3408300/1‑100 (MQ=255)
                              gTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACg                                                                       <  1:332854/100‑1 (MQ=255)
                              gTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACg                                                                       <  1:3700656/100‑1 (MQ=255)
                              gTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACg                                                                       <  2:737894/100‑1 (MQ=255)
                               tCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGt                                                                      <  1:3718958/100‑1 (MQ=255)
                                 aaaTTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCg                                                                    <  1:1530586/100‑1 (MQ=255)
                                 aaaTTGCC‑TTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGa                                                                   <  1:2303158/100‑1 (MQ=255)
                                   aTTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGAt                                                                  <  2:1158253/100‑1 (MQ=255)
                                      gCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACt                                                               >  1:1991078/1‑100 (MQ=255)
                                       ccTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTc                                                              >  2:2177383/1‑100 (MQ=255)
                                        cTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTACTTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCg                                                             <  2:2862337/100‑1 (MQ=255)
                                        cTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCg                                                             <  2:2823285/100‑1 (MQ=255)
                                          tttGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGcc                                                           <  2:2487474/100‑1 (MQ=255)
                                            tGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCtn                                                         <  2:3815259/100‑2 (MQ=255)
                                            tGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTg                                                         <  2:3660386/100‑1 (MQ=255)
                                               ccTGAAAAATTCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCg                                                      <  2:535937/100‑1 (MQ=255)
                                               ccTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCg                                                      >  2:3705187/1‑100 (MQ=255)
                                                cTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGt                                                     >  2:260613/1‑100 (MQ=255)
                                                cTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGt                                                     <  1:3788853/100‑1 (MQ=255)
                                                 tGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTc                                                    >  1:3106281/1‑100 (MQ=255)
                                                    aaaaTGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATc                                                 >  2:359050/1‑100 (MQ=255)
                                                         gCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCaaa                                            >  2:3183917/1‑100 (MQ=255)
                                                           aTAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGc                                          <  2:2746295/100‑1 (MQ=255)
                                                            tAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCg                                         >  1:1622931/1‑100 (MQ=255)
                                                                 aaaTATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGc                                     <  2:2706646/99‑1 (MQ=255)
                                                                       cTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATa                              <  1:2165420/100‑1 (MQ=255)
                                                                         gTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATAcc                            >  2:692161/1‑100 (MQ=255)
                                                                         gTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATAcc                            >  1:1979332/1‑100 (MQ=255)
                                                                          ngctAACTAATTGGCGTTGC‑ATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCg                           >  1:2198600/4‑100 (MQ=25)
                                                                           gctAACTAATTGGCGTTGC‑ATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGt                          >  2:1707091/3‑100 (MQ=25)
                                                                            ttAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGtt                         >  1:306986/1‑100 (MQ=255)
                                                                                cTAATTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaa                     <  2:1487194/100‑1 (MQ=255)
                                                                                cTAATTGGCGTTGC‑A‑GTACATGCAACGCCAATAAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaa                     >  1:869244/1‑100 (MQ=255)
                                                                                 tAATTGGCGTTGC‑ATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaaa                    >  1:3752296/1‑100 (MQ=38)
                                                                                  aaTTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaaaa                   >  1:1194026/1‑100 (MQ=255)
                                                                                  aaTTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaaaa                   >  1:324563/1‑100 (MQ=255)
                                                                                  aaTTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaaaa                   >  1:381772/1‑100 (MQ=255)
                                                                                   aTTGGCGTTGC‑ATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt                  >  2:185702/1‑100 (MQ=38)
                                                                                   aTTGGCGTTGC‑ATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt                  >  1:2802442/1‑100 (MQ=38)
                                                                                   aTTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt                  >  1:3529815/1‑100 (MQ=255)
                                                                                    ttGGCGTTGCGT‑GTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATg                 >  1:3087707/1‑100 (MQ=38)
                                                                                    ttGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATg                 <  1:369820/100‑1 (MQ=255)
                                                                                     tGGCGTTGC‑ATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGa                >  1:1820055/1‑100 (MQ=38)
                                                                                     tGGCGTTGC‑ATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGa                >  1:3156316/1‑100 (MQ=38)
                                                                                     tGGCGTTGC‑ATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGa                >  1:638124/1‑100 (MQ=38)
                                                                                     tGGCGTTGC‑ATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGa                >  1:931138/1‑100 (MQ=38)
                                                                                     tGGCGTTGC‑ATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGa                >  1:674461/1‑100 (MQ=38)
                                                                                     tGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaa                     >  2:1253325/1‑95 (MQ=255)
                                                                                      ggCGTTGC‑ATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc               >  1:100267/1‑100 (MQ=38)
                                                                                      ggCGTTGC‑ATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc               >  1:1241159/1‑100 (MQ=38)
                                                                                      ggCGTTGC‑ATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc               >  1:3175886/1‑100 (MQ=38)
                                                                                      ggCGTTGC‑ATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc               >  2:1189871/1‑100 (MQ=38)
                                                                                      ggCGTTGC‑ATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc               >  2:2517887/1‑100 (MQ=38)
                                                                                      ggCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc               >  2:3845001/1‑100 (MQ=255)
                                                                                      ggCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc               <  2:1979318/100‑1 (MQ=255)
                                                                                       gCGTTGC‑ATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAct              >  2:1047457/1‑100 (MQ=38)
                                                                                           tGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTgg          <  2:3028601/100‑1 (MQ=255)
                                                                                           tGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTgg          >  2:184883/1‑100 (MQ=255)
                                                                                           tGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTgg          >  1:1409126/1‑100 (MQ=255)
                                                                                           tGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTgg          <  2:782500/100‑1 (MQ=255)
                                                                                            gC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACTTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGTCTCTGGt         >  2:2082154/1‑100 (MQ=255)
                                                                                            gC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGt         >  2:2353921/1‑100 (MQ=255)
                                                                                               a‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCa       >  2:3210812/1‑100 (MQ=255)
                                                                                               a‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCa       <  2:1651594/100‑1 (MQ=255)
                                                                                               a‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAAGACTCTg           >  2:1003314/1‑96 (MQ=255)
                                                                                                 gTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAg      >  2:1786549/1‑100 (MQ=255)
                                                                                                     aTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGaaaa  <  1:1408836/100‑1 (MQ=255)
                                                                                                     |                                                                                                   
CTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGC‑A‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAA  >  NC_000913/1446263‑1446461

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