Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,196,356 G→C 5.2% A151A (GCC→GCG rsd ← stationary phase protein, binds sigma 70 RNA polymerase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,196,3560GC 5.2% 368.8 / 7.5 110A151A (GCC→GCGrsdstationary phase protein, binds sigma 70 RNA polymerase subunit
Reads supporting (aligned to +/- strand):  ref base G (62/42);  new base C (4/2);  total (66/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 115 total aligned reads displayed.

TATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTG  >  NC_000913/4196209‑4196505
                                                                                                                                                   |                                                                                                                                                     
tATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGCCGCGggcggc                                                                                                                                                     >  1:154263/1‑150 (MQ=255)
    aGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCa                                                                                                                                                 >  2:1746510/1‑150 (MQ=255)
    aGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCAGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCa                                                                                                                                                 >  2:357212/1‑150 (MQ=255)
    aGCCTGATCAGGTTCCGCGGATCCCGAATAAACGATCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCa                                                                                                                                                 >  1:1532937/1‑150 (MQ=255)
     gCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTcta                                                                                                                                                >  2:1637129/1‑148 (MQ=255)
       cTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGTCGGCGACAAGc                                                                                                                                              >  1:64299/1‑150 (MQ=255)
       cTGATCAGGTTCCGCGGATCCCGAATAAACGGGCTCCGCCAGCTAAGGCACTCCGACAAGAAATAAGCCCCGCAAATGGGGCATTGAATGTAACTTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGc                                                                                                                                              <  1:1681914/150‑1 (MQ=255)
        tGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCa                                                                                                                                             >  1:715688/1‑150 (MQ=255)
        tGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCa                                                                                                                                             >  1:1725736/1‑150 (MQ=255)
         gATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCCGCCCGCGTCAAACat                                                                                                                                            <  1:736649/150‑2 (MQ=255)
           tCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCCAGCAGGATGTTTGACGCTGGCGGCGTCAAGCACca                                                                                                                                          >  2:1404061/1‑150 (MQ=255)
            cAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACcag                                                                                                                                         <  1:202536/150‑1 (MQ=255)
             aGGCTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCCGCCCGCGTCAAACatcctg                                                                                                                                        <  1:1345700/150‑6 (MQ=255)
                    gcgGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCa                                                                                                                                 <  2:1452303/150‑1 (MQ=255)
                    gcgGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCa                                                                                                                                 <  2:1459691/150‑1 (MQ=255)
                      ggATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGc                                                                                                                               <  1:1013998/150‑1 (MQ=255)
                         tCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTa                                                                                                                            >  1:1070055/1‑150 (MQ=255)
                           ccGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATc                                                                                                                          >  1:262531/1‑150 (MQ=255)
                           ccGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATc                                                                                                                          >  1:1149965/1‑150 (MQ=255)
                           ccGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGAAGCGGGCGGCGTCAAGCACAAGCAGAATCCGCTTATc                                                                                                                          >  1:1395841/1‑150 (MQ=255)
                            cGAATAAACGGTCTCAGCCAGTTAAGGCCCTCCGACAAGAAATTAGACGCGCAAATGGGGCATTGAATGTACAGTACCCGGTAACAGCGCAGAACTCCAGCCGGAGGTTTGACGGGGGCCGCGGCACGCACCAGCAGAATCAGCTTATca                                                                                                                         >  1:300845/1‑149 (MQ=255)
                            cGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCt                                                                                                                         >  1:727117/1‑150 (MQ=255)
                            cGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCt                                                                                                                         >  2:1300302/1‑150 (MQ=255)
                                 aaaCGGTCTCAGCCAGGTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAg                                                                                                                    <  1:1747343/150‑1 (MQ=255)
                                  aaCGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGc                                                                                                                   >  2:1490448/1‑150 (MQ=255)
                                      gTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACaa                                                                                                               >  1:1753057/1‑150 (MQ=255)
                                        ctcAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAg                                                                                                             >  1:429825/1‑150 (MQ=255)
                                          cAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCCGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAcgcg                                                                                                           >  2:1349538/1‑150 (MQ=255)
                                          cAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCCAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAgcg                                                                                                           >  2:1185295/1‑150 (MQ=255)
                                            gCCAGTTAATGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCATGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAgcgcg                                                                                                         <  1:112719/150‑1 (MQ=255)
                                            gCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAgcgcg                                                                                                         >  1:1726244/1‑150 (MQ=255)
                                            gCCAGTTAAGGCACGCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAgcgcg                                                                                                         <  2:1014010/150‑1 (MQ=255)
                                             ccAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAgcgcgc                                                                                                        <  1:1160249/150‑1 (MQ=255)
                                             ccAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAgcgcgc                                                                                                        >  1:462992/1‑150 (MQ=255)
                                              cAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCt                                                                                                       <  2:1507105/150‑1 (MQ=255)
                                                 ttAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCATGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTcc                                                                                                    <  2:189649/150‑1 (MQ=255)
                                                 ttAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCGTATCTTCCCGCACAAAGCGCGCTTcc                                                                                                    >  2:890048/1‑150 (MQ=255)
                                                    aGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAgc                                                                                                 <  2:1222909/150‑1 (MQ=255)
                                                       cACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGATTGTAAATTACGCGTTAACAGCGCAGAACTCAATCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCTCAAAGAGTTCTTCCAGCGcc                                                                                              >  2:272046/1‑150 (MQ=255)
                                                       cACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGcc                                                                                              <  1:1051839/150‑1 (MQ=255)
                                                       cACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGcc                                                                                              <  1:335290/150‑1 (MQ=255)
                                                           ccGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGc                                                                                          >  2:664994/1‑150 (MQ=255)
                                                              aCAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGa                                                                                       <  2:926284/150‑1 (MQ=255)
                                                               cAAGAAATTAGCCCCGCAAATGGGGCAATGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGAt                                                                                      >  1:1213331/1‑150 (MQ=255)
                                                                aaGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATg                                                                                     <  1:1467986/150‑1 (MQ=255)
                                                                 aGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGCATCAGATTCTCTTCCAGCACAAAGCGCGCTGCCAGCGCCTCGCCgcggg                                                                                    >  1:1025235/1‑146 (MQ=255)
                                                                  gAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTc                                                                                   >  1:1172738/1‑150 (MQ=255)
                                                                    aaTTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAg                                                                                 <  1:513896/150‑1 (MQ=255)
                                                                    aaTTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAg                                                                                 >  2:1761118/1‑150 (MQ=255)
                                                                      ttAGCCCCGCAACTGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGAt                                                                               <  2:884819/150‑1 (MQ=255)
                                                                           cccGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAAc                                                                          >  1:130467/1‑150 (MQ=255)
                                                                           cccGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAAc                                                                          >  1:1780699/1‑150 (MQ=255)
                                                                           cccGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGAGAAAAc                                                                          >  2:667323/1‑150 (MQ=255)
                                                                            ccGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAAcc                                                                         >  2:1467935/1‑150 (MQ=255)
                                                                             cGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCGGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGAGGTCAGATAAAGCCt                                                                        >  2:336883/1‑150 (MQ=255)
                                                                              gCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTg                                                                       <  1:1065936/150‑1 (MQ=255)
                                                                                aaaTGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGtt                                                                     >  2:735759/1‑150 (MQ=255)
                                                                                 aaTGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATGAAACCTGTTg                                                                    <  2:1725736/150‑1 (MQ=255)
                                                                                  aTGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGa                                                                   <  2:602507/150‑1 (MQ=255)
                                                                                     gggCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAAc                                                                <  2:1028562/150‑1 (MQ=255)
                                                                                         aTTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCaa                                                            >  2:399810/1‑150 (MQ=255)
                                                                                         aTTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCaa                                                            >  1:1611254/1‑150 (MQ=255)
                                                                                           tGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAgg                                                          >  1:956423/1‑150 (MQ=255)
                                                                                           tGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAgg                                                          <  1:1259/150‑1 (MQ=255)
                                                                                            gAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGt                                                         >  2:960202/1‑150 (MQ=255)
                                                                                                 tAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTa                                                    >  2:1149919/1‑150 (MQ=255)
                                                                                                    aTTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCa                                                 >  1:1515407/1‑150 (MQ=255)
                                                                                                     ttACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGCCGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCGTCCAGCGCCTCGCAGATGTCAGATAAAACCTGTGGACACTCAAGGTAATGATCag                                                >  1:267484/1‑149 (MQ=255)
                                                                                                     ttACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCAACAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCAt                                                >  1:267489/1‑150 (MQ=255)
                                                                                                        cgcgTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGAt                                             >  1:1021782/1‑150 (MQ=255)
                                                                                                        cgcgTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGAt                                             >  1:832355/1‑150 (MQ=255)
                                                                                                         gcgTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCTCGCTTCCAGCGCCTCGCCGATGGCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATc                                            <  1:364417/150‑1 (MQ=255)
                                                                                                         gcgTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATc                                            >  2:1716624/1‑150 (MQ=255)
                                                                                                         gcgTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATc                                            <  1:389094/150‑1 (MQ=255)
                                                                                                         gcgTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATc                                            <  2:262531/150‑1 (MQ=255)
                                                                                                               aCAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAAACTGTTGAAACTCAAGGTAATTATCATGATCGATGGc                                      >  2:1304072/1‑150 (MQ=255)
                                                                                                                cAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCGTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTGGAAACTCACGGTAAGTATCAGGAGCGATGTCg                                     >  1:275319/1‑150 (MQ=255)
                                                                                                                 aGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCgg                                    >  2:1743547/1‑150 (MQ=255)
                                                                                                                    gcAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGttt                                 >  1:1758781/1‑150 (MQ=255)
                                                                                                                      aGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTcc                               <  1:1793550/150‑1 (MQ=255)
                                                                                                                      aGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCACCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTcc                               <  2:1149965/150‑1 (MQ=255)
                                                                                                                      aGAAATCAAGCAGGATGTTTGCCGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCAAAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAACCCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGGTTcc                               <  2:899554/150‑1 (MQ=255)
                                                                                                                       gAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTGCCCGCGCCACGCCGAGGTCAGATAAAAGCTGTTGAAACTGAAGGTAATTATAATGTTCGATGGCGGTTTACa                              >  2:267147/1‑150 (MQ=255)
                                                                                                                            cAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTTTCATGATCGATGGCGGTTTCCAGACTg                         >  1:1571435/1‑150 (MQ=255)
                                                                                                                            cAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTAGCATGATCGATGGCGGTTTCCAGACTg                         >  1:1262666/1‑150 (MQ=255)
                                                                                                                             aaGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTgg                        <  2:1379047/150‑1 (MQ=255)
                                                                                                                             aaGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTgg                        <  2:220888/150‑1 (MQ=255)
                                                                                                                               gCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGaa                      >  1:1859348/1‑150 (MQ=255)
                                                                                                                                 aGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATc                    >  1:875415/1‑150 (MQ=255)
                                                                                                                                   ggtgcttGACGCCGCCCGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCgt                  >  2:1558256/6‑150 (MQ=255)
                                                                                                                                       tttGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCgtagta              >  2:1860305/1‑150 (MQ=255)
                                                                                                                                       cttGACGCCGCCCGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCCGCGCCTCGCCGATGTCAGATTAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCgtagta              >  1:354070/2‑150 (MQ=255)
                                                                                                                                       cttGACGCCGCCCGCGTCAAGCACCAGCAGAATCAGCTTATCTACCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAACACCTGTTGAAACTCAAGGTAATTATAATGATCGATGGCGGTTTCAAGACTGGAATCgtagta              >  1:600717/2‑150 (MQ=255)
                                                                                                                                        ttGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTaa             >  1:1819484/1‑150 (MQ=255)
                                                                                                                                          gACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGAAGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATTGTAGTAATc           <  2:1726244/150‑1 (MQ=255)
                                                                                                                                            cgcgGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCa         >  1:1416218/1‑150 (MQ=255)
                                                                                                                                                ggcggcGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGAGCCTCGCCGATGTCAGATAAAACCTGTAGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAAt     >  2:971852/1‑150 (MQ=255)
                                                                                                                                                  cggcgTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCt   >  1:706456/1‑150 (MQ=255)
                                                                                                                                                   ggcgTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTg  <  2:1532937/150‑1 (MQ=255)
                                                                                                                                                   ggcgTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTg  <  1:1149919/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTG  >  NC_000913/4196209‑4196505

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