Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,992,0520TA68.9% 160.2 / 110.1 132D300E (GAT→GAAcyaAadenylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (19/22);  new base A (47/44);  total (66/66)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.37e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 133 total aligned reads displayed.

TCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCGTAGGCTGGCGTCGC  >  NC_000913/3991904‑3992196
                                                                                                                                                    |                                                                                                                                                
tctACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAt                                                                                                                                                  >  1:340495/1‑149 (MQ=255)
  tACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATc                                                                                                                                                 >  8:283889/1‑148 (MQ=255)
     aaGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAACCAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATa                                                                                                                                             >  6:111436/1‑149 (MQ=255)
     aaGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATa                                                                                                                                             <  6:363324/149‑1 (MQ=255)
     aaGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATa                                                                                                                                             <  1:127729/149‑1 (MQ=255)
          tATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCa                                                                                                                                        >  2:169033/1‑149 (MQ=255)
          tATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCa                                                                                                                                        <  4:307618/149‑1 (MQ=255)
          tATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCa                                                                                                                                        >  6:172054/1‑149 (MQ=255)
             cGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCatga                                                                                                                                     >  1:7734/1‑149 (MQ=255)
               aTTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCatgatg                                                                                                                                   >  3:260753/1‑149 (MQ=255)
                   cccATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTg                                                                                                                                >  3:352059/1‑148 (MQ=255)
                   cccATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTg                                                                                                                                >  6:126013/1‑148 (MQ=255)
                   cccATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTg                                                                                                                                >  4:88244/1‑148 (MQ=255)
                   cccATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTg                                                                                                                                >  5:341211/1‑148 (MQ=255)
                    ccATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGa                                                                                                                              <  6:337246/149‑1 (MQ=255)
                    ccATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAAACCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTg                                                                                                                                >  4:85042/1‑147 (MQ=255)
                     cATACAGAGCGGTACTGAAACCACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGATCCCCCGTCTGCTGGCGAAAGATATCACACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAg                                                                                                                             >  6:159981/1‑149 (MQ=255)
                     cATACAGAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGATCCCACGTCTGCTGGCGAAAGATATCAACCAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAg                                                                                                                             >  6:165909/1‑149 (MQ=255)
                      aTACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGc                                                                                                                            <  3:88244/149‑1 (MQ=255)
                         cAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCgtg                                                                                                                         <  5:147544/149‑1 (MQ=255)
                             gCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTAc                                                                                                                     <  2:69418/149‑1 (MQ=255)
                              cGGTACTGAAAACACTGCTGCTGGAAGGCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGGATGATGCTGGAGCGTGTTACt                                                                                                                    <  6:36105/149‑1 (MQ=255)
                                 tACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGACCCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGCCGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGaa                                                                                                                 >  2:37776/1‑149 (MQ=255)
                                 tACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGaa                                                                                                                 >  1:140397/1‑149 (MQ=255)
                                 tACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGaa                                                                                                                 <  7:357042/149‑1 (MQ=255)
                                 tACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGaa                                                                                                                 >  5:95733/1‑149 (MQ=255)
                                  aCTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGaa                                                                                                                 >  1:124649/1‑148 (MQ=255)
                                   cTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATa                                                                                                                                             <  5:195676/119‑1 (MQ=255)
                                   cTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGTGTGTTACTGAAt                                                                                                                <  3:212784/148‑1 (MQ=255)
                                   cTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATa                                                                                                               <  7:352690/149‑1 (MQ=255)
                                    tGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGACTAc                                                                                                              >  5:384183/1‑149 (MQ=255)
                                    tGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATAc                                                                                                              >  3:278656/1‑149 (MQ=255)
                                      aaaaCACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCt                                                                                                            >  3:129223/1‑149 (MQ=255)
                                       aaaCACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCt                                                                                                            >  1:170317/1‑148 (MQ=255)
                                        aaCACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGa                                                                                                          <  7:26318/149‑1 (MQ=255)
                                           acTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACgg                                                                                                       <  1:169033/149‑1 (MQ=255)
                                           acTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACgg                                                                                                       >  3:214062/1‑149 (MQ=255)
                                               ctgctgGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGa                                                                                                    <  1:7997/148‑1 (MQ=255)
                                                     gAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGa                                                                                             <  1:84422/149‑1 (MQ=255)
                                                     gAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTg                                                                                                 >  1:320111/1‑145 (MQ=255)
                                                       aGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGAtt                                                                                           <  3:6248/149‑1 (MQ=255)
                                                        gCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGAttt                                                                                          <  5:111436/149‑1 (MQ=255)
                                                         ccTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGACGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGAttt                                                                                          >  2:101025/1‑148 (MQ=255)
                                                           tATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTAc                                                                                       >  4:309344/1‑149 (MQ=255)
                                                           tATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGCTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTAc                                                                                       >  3:34550/1‑149 (MQ=255)
                                                             ttCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCAGGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTAc                                                                                       <  5:172054/147‑1 (MQ=255)
                                                             ttCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTAccc                                                                                     >  8:360213/1‑149 (MQ=255)
                                                             ttCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTAcc                                                                                      <  4:260753/148‑1 (MQ=255)
                                                                cTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTc                                                                                  >  4:154316/1‑149 (MQ=255)
                                                                cTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTc                                                                                  >  1:200180/1‑149 (MQ=255)
                                                                   ggAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCt                                                                                 <  2:340495/147‑1 (MQ=255)
                                                                       tACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGAttt                                                                           <  2:140397/149‑1 (MQ=255)
                                                                         cccGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTa                                                                          <  2:170317/148‑1 (MQ=255)
                                                                         cccGAACCCACGTCTGCCGGCGAAAGATATCAAACAGCGCTTGCACGACGGCGAGATTGTATCTGTTGGTCTCGAACCATACTGTATGATGCTGGAGCGTGCTACTGAATACCTGTCGGCGATTGAAGATTTTACACGTCTGGATTTa                                                                          <  2:222218/148‑1 (MQ=255)
                                                                          ccGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGCACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGt                                                                        >  7:312115/1‑149 (MQ=255)
                                                                          ccGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAg                                                                         >  6:365909/1‑148 (MQ=255)
                                                                          ccGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAg                                                                         >  6:401525/1‑148 (MQ=255)
                                                                               cccACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTc                                                                   <  7:283889/149‑1 (MQ=255)
                                                                                ccACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCg                                                                  >  1:234430/1‑149 (MQ=255)
                                                                                ccACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCg                                                                  >  4:246647/1‑149 (MQ=255)
                                                                                 cACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGa                                                                                                                  <  6:256826/100‑1 (MQ=255)
                                                                                  aagTCTGCTGGCGAAAGATATTTAACAGCGTTTTCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCgct                                                                <  2:94938/147‑1 (MQ=255)
                                                                                  aCGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCgct                                                                >  4:129920/1‑149 (MQ=255)
                                                                                  aCGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATAACTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCgct                                                                >  6:402941/1‑149 (MQ=255)
                                                                                          tgGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTAc                                                                                       >  4:201226/1‑118 (MQ=255)
                                                                                          tgGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTAt                                                        <  2:200180/149‑1 (MQ=255)
                                                                                           ggCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTAtt                                                       <  1:37776/149‑1 (MQ=255)
                                                                                                aaGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTaaaa                                                  <  3:129920/149‑1 (MQ=255)
                                                                                                aaGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTaaaa                                                  <  7:360213/149‑1 (MQ=255)
                                                                                                aaGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTaaaa                                                  <  6:82473/149‑1 (MQ=255)
                                                                                                aaGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTaaaa                                                  <  5:402941/149‑1 (MQ=255)
                                                                                                  gATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAgt                                                >  2:55844/1‑149 (MQ=255)
                                                                                                  gATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGACGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAgt                                                <  1:101025/149‑1 (MQ=255)
                                                                                                       cAAACATCGTTTGCACGACGGCGATATTGTATCGTTTGTTCTCGAAGCATACTGCATGATGCTGGAGCGTGTTACTGAATACATTACTTCGATTGAAGAttt                                                                                          <  3:235671/102‑1 (MQ=255)
                                                                                                       cAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGc                                            >  7:92774/1‑148 (MQ=255)
                                                                                                        aaaCAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCg                                           >  4:356915/1‑148 (MQ=255)
                                                                                                          aCAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGaa                                         >  4:355346/1‑148 (MQ=255)
                                                                                                          aCAGCGTTTGCACGACGGCGAGATTGTATCGTATGGTCTCGAACCATACTGCATGATGCTCGAGCGTGTTACTGAATAACTGCCGGCGATTGAAGATTTTACCCGTCTGTATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGaa                                         >  7:419580/1‑148 (MQ=255)
                                                                                                           cAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGaa                                         >  4:21991/1‑147 (MQ=255)
                                                                                                           cAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGaaaa                                       <  4:35922/149‑1 (MQ=255)
                                                                                                               gTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTc                                   <  2:124649/149‑1 (MQ=255)
                                                                                                                tttGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCa                                  <  3:85042/149‑1 (MQ=255)
                                                                                                                    cACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTAt                                                        <  2:34418/123‑1 (MQ=255)
                                                                                                                    cACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCg                              >  8:111111/1‑149 (MQ=255)
                                                                                                                           gCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAAcgc                       <  4:34550/149‑1 (MQ=255)
                                                                                                                               gaTTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATGGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCt                   <  6:95733/149‑1 (MQ=255)
                                                                                                                                  tGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGc                 <  2:53495/148‑1 (MQ=255)
                                                                                                                                   gTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCg                >  5:384163/1‑148 (MQ=255)
                                                                                                                                    tATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGt                             <  6:29990/134‑1 (MQ=255)
                                                                                                                                    tATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGt                             >  5:29990/1‑134 (MQ=255)
                                                                                                                                    tATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCGt               <  2:234430/148‑1 (MQ=255)
                                                                                                                                    tATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCg                >  2:317603/1‑147 (MQ=255)
                                                                                                                                    tATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCGTa              <  7:111111/149‑1 (MQ=255)
                                                                                                                                    tATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCGTa              <  3:309344/149‑1 (MQ=255)
                                                                                                                                      tCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCGTAg             <  3:356915/148‑1 (MQ=255)
                                                                                                                                        gTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTCCCCGTCTGGATTTAGTACGTCGCGGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCGTAGGc           <  8:312115/148‑1 (MQ=255)
                                                                                                                                           tGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCGTAGGCTgg        <  5:365909/148‑1 (MQ=255)
                                                                                                                                            ggTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCGTAGGCTGGc       <  6:384163/148‑1 (MQ=255)
                                                                                                                                            ggTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCGTAGGCTGGc       <  1:46327/148‑1 (MQ=255)
                                                                                                                                               ctcGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCGTAGGCTGGCGTcg   >  1:42447/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
TCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCGTAGGCTGGCGTCGC  >  NC_000913/3991904‑3992196

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: