Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,770,649 A→C 6.7% L244W (TTG→TGG) ‡ yibH ← inner membrane protein YibH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,770,6490AC 6.7% 813.7 / 40.1 254L244W (TTG→TGG) ‡yibHinner membrane protein YibH
Reads supporting (aligned to +/- strand):  ref base A (110/127);  new base C (14/3);  total (124/130)
Fisher's exact test for biased strand distribution p-value = 4.91e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.60e-01

Only 100 of 257 total aligned reads displayed.

GTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGATTGCTCGCCGTTAACCAT  >  NC_000913/3770500‑3770797
                                                                                                                                                     |                                                                                                                                                    
gTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCa                                                                                                                                                      >  2:530402/1‑150 (MQ=255)
 tttGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCACCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCaa                                                                                                                                                     <  2:1208398/150‑1 (MQ=255)
    gCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCATGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGc                                                                                                                                                  <  2:240761/150‑1 (MQ=255)
    gCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGc                                                                                                                                                  >  1:678228/1‑150 (MQ=255)
     ccGTGAAACACCTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCa                                                                                                                                                 >  2:824051/1‑150 (MQ=255)
      cGTGAAACACCTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAg                                                                                                                                                >  2:143409/1‑150 (MQ=255)
       gTGAAACACCTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGc                                                                                                                                               >  1:661348/1‑150 (MQ=255)
        tGAAACACCTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCt                                                                                                                                              >  1:362159/1‑150 (MQ=255)
          aaaCACCTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGc                                                                                                                                            <  2:1303271/150‑1 (MQ=255)
             cacCTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTa                                                                                                                                         >  2:1086379/1‑150 (MQ=255)
              acCTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCTGCCTTGCCGCTGCGtcc                                                                                                                                        <  1:763695/150‑3 (MQ=38)
               ccTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATg                                                                                                                                       >  1:1106999/1‑150 (MQ=255)
                 tGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTa                                                                                                                                     >  1:119791/1‑150 (MQ=255)
                   cccAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTAcc                                                                                                                                   <  1:1423304/150‑1 (MQ=255)
                     cAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTg                                                                                                                                 <  1:538239/150‑1 (MQ=255)
                     cAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTg                                                                                                                                 >  2:212691/1‑150 (MQ=255)
                      aGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTgg                                                                                                                                <  1:1400048/150‑1 (MQ=255)
                        gTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGc                                                                                                                              >  2:1959430/1‑150 (MQ=255)
                         tAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCg                                                                                                                             >  1:925060/1‑150 (MQ=255)
                            cgcgTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGAt                                                                                                                          <  1:530402/150‑1 (MQ=255)
                             gcgTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATc                                                                                                                         <  1:824051/150‑1 (MQ=255)
                              cgTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCa                                                                                                                        >  1:1693242/1‑150 (MQ=255)
                                     caccacTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTg                                                                                                                 <  2:544957/150‑1 (MQ=255)
                                       ccacTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAg                                                                                                               <  1:1642337/150‑1 (MQ=255)
                                         acTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTa                                                                                                             <  2:49312/150‑1 (MQ=255)
                                            tCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACg                                                                                                          <  2:334362/150‑1 (MQ=255)
                                             ccGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCTGCCTTGCCGCTGCGtc                                                                                                                                         <  1:1680346/118‑2 (MQ=255)
                                             ccGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGt                                                                                                         <  2:1044680/150‑1 (MQ=255)
                                              cGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTTAGTAACGTa                                                                                                        <  2:1863933/150‑1 (MQ=255)
                                                  tcatcGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTa                                                                                                                                         >  1:78452/1‑113 (MQ=255)
                                                  tcatcGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTa                                                                                                                                         <  2:78452/113‑1 (MQ=255)
                                                  tcatcGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCa                                                                                                    >  1:29489/1‑150 (MQ=255)
                                                   catcGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCAt                                                                                                   >  1:458276/1‑150 (MQ=255)
                                                    atcGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCAtt                                                                                                  <  1:1959430/150‑1 (MQ=255)
                                                      cGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGc                                                                                                >  2:709855/1‑150 (MQ=255)
                                                           ggTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGgc                                                                                           >  1:1474216/1‑150 (MQ=255)
                                                            gTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGgcg                                                                                          >  1:964452/1‑150 (MQ=255)
                                                            gTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCAGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGAGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGgcg                                                                                          <  1:1027042/150‑1 (MQ=255)
                                                                cAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGcgcgcg                                                                                      <  2:1604657/150‑1 (MQ=255)
                                                                cAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGcgcgcg                                                                                      >  2:1421299/1‑150 (MQ=255)
                                                                   aCGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAAt                                                                                   <  2:1658724/150‑1 (MQ=255)
                                                                   aCGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAAt                                                                                   <  2:262102/150‑1 (MQ=255)
                                                                   aCGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAAt                                                                                   >  2:1366955/1‑150 (MQ=255)
                                                                    cGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATg                                                                                  <  1:1890558/150‑1 (MQ=255)
                                                                      tAACAGCGACTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGAc                                                                                >  2:1683626/1‑150 (MQ=255)
                                                                       aaCAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACa                                                                               >  1:1809773/1‑150 (MQ=255)
                                                                       aaCAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACa                                                                               >  1:266944/1‑150 (MQ=255)
                                                                        aCAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAg                                                                              <  1:1225734/150‑1 (MQ=255)
                                                                         cAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGt                                                                             >  1:1951165/1‑150 (MQ=255)
                                                                          aGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTc                                                                            >  2:498888/1‑150 (MQ=255)
                                                                             gAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGACGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGc                                                                         >  1:892757/1‑150 (MQ=255)
                                                                              aGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCt                                                                        <  1:875322/150‑1 (MQ=255)
                                                                                ttttGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCa                                                                      <  1:644112/150‑1 (MQ=255)
                                                                                   tgtcGAAATTGGGCGACAATTTGCCGTTTTTTCTCGGGGATGAAGACCATCACCGGACGCTGCTGCAAGGCAGCTGCGTATTTACCTGGGCGGATCAGTACCTTAGTAACGTAGCCTTTGCTCGGCGCGCGAATGACAGTCTGCTCAAgg                                                                   <  2:568311/147‑1 (MQ=255)
                                                                                   tGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAgg                                                                   <  2:119791/150‑1 (MQ=255)
                                                                                    gCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGt                                                                  <  1:1608168/150‑1 (MQ=255)
                                                                                    gCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGt                                                                  >  2:647646/1‑150 (MQ=255)
                                                                                    gCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGt                                                                  <  1:438985/150‑1 (MQ=255)
                                                                                     ccGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGtt                                                                 >  1:1838282/1‑150 (MQ=255)
                                                                                       gAAATTTGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGtttt                                                               <  2:195150/150‑3 (MQ=255)
                                                                                        aaaTTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTtata                                                              <  1:1735824/150‑1 (MQ=255)
                                                                                             gggCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTg                                                         <  1:1013936/150‑1 (MQ=255)
                                                                                             gggCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTg                                                         >  2:178600/1‑150 (MQ=255)
                                                                                              ggCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCTGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGc                                                        <  2:938730/150‑1 (MQ=255)
                                                                                              ggCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAg                                                                                                                       >  1:251476/1‑87 (MQ=255)
                                                                                                cGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCtt                                                      <  1:215005/150‑1 (MQ=255)
                                                                                                cGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCtt                                                      >  2:1263632/1‑150 (MQ=255)
                                                                                                 gACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTc                                                     <  2:1379264/150‑1 (MQ=255)
                                                                                                 gACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTc                                                     <  2:1935845/150‑1 (MQ=255)
                                                                                                   cAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACAGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAg                                                   >  1:1264498/1‑150 (MQ=255)
                                                                                                     aTTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTa                                                 <  2:877073/150‑1 (MQ=255)
                                                                                                     aTTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTa                                                 <  2:458276/150‑1 (MQ=255)
                                                                                                          ccGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTg                                            <  1:560033/150‑1 (MQ=255)
                                                                                                             tttttGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTgcgc                                         <  2:214452/150‑1 (MQ=255)
                                                                                                                  gCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTa                                    <  1:498888/150‑1 (MQ=255)
                                                                                                                    tCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAg                                  >  2:1897916/1‑150 (MQ=255)
                                                                                                                     cGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGc                                 <  1:1086379/150‑1 (MQ=255)
                                                                                                                      ggggATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCt                                <  2:425419/150‑1 (MQ=255)
                                                                                                                        ggATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTca                              <  1:210553/150‑1 (MQ=255)
                                                                                                                        ggATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTca                              >  1:1182612/1‑150 (MQ=255)
                                                                                                                         gATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTcac                             <  1:1771857/150‑1 (MQ=255)
                                                                                                                             aaGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATc                         >  2:621661/1‑150 (MQ=255)
                                                                                                                               gACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTg                       >  2:1771039/1‑150 (MQ=255)
                                                                                                                                 ccATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTCGGAGGAACAGAACCAGAGTAACGTAGACATTGCTCGGCGCGCGAAAGACAGACTGCTCAAGGTTATAATTTGCTTAAATAAGAAGCGCTCTTAAGCTCACAAACAGCg                     >  2:1840888/1‑150 (MQ=255)
                                                                                                                                   aTCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGAt                   <  1:395812/150‑1 (MQ=255)
                                                                                                                                   aTCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGAt                   <  1:1366955/150‑1 (MQ=255)
                                                                                                                                   aTCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGAt                   <  2:1729326/150‑1 (MQ=255)
                                                                                                                                   aTCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGAt                   <  1:178600/150‑1 (MQ=255)
                                                                                                                                     cACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGATTg                 >  2:900146/1‑150 (MQ=255)
                                                                                                                                        cataCGCAGCTGCCTTGCCGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGATTGCTc              >  2:84733/4‑150 (MQ=38)
                                                                                                                                         ataCGCAGCTGCCTTGCCGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGATTGCTCg             >  1:1601966/3‑150 (MQ=38)
                                                                                                                                         ataCGCAGCTGCCTTGCCGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGATCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGATTGCTCg             >  2:1831950/3‑150 (MQ=38)
                                                                                                                                           aCGCAGCTGCCTTGCCGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGATTGCTCGcc           <  1:81548/150‑1 (MQ=39)
                                                                                                                                            cGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGATTGCTCGCCg          >  1:573107/1‑150 (MQ=255)
                                                                                                                                            cGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGATTGCTCGCCg          <  1:621661/150‑1 (MQ=255)
                                                                                                                                               aGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGATTGCTCGCCGTTa       >  2:1224282/1‑150 (MQ=255)
                                                                                                                                                 cGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTc                                                     <  2:364/102‑1 (MQ=255)
                                                                                                                                                 cGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGATTGCTCGCCGTTAAc     <  2:29489/150‑1 (MQ=255)
                                                                                                                                                 cGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGATTGCTCGCCGTTAAc     >  2:1639602/1‑150 (MQ=255)
                                                                                                                                                    cAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGATTGCTCGCCGTTAACCAt  >  2:434545/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGATTGCTCGCCGTTAACCAT  >  NC_000913/3770500‑3770797

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