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% 314.8 / 7.0 96A151A (GCC→GCGrsdstationary phase protein, binds sigma 70 RNA polymerase subunit
Reads supporting (aligned to +/- strand):  ref base G (41/50);  new base C (3/2);  total (44/52)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.33e-01

Only 100 of 101 total aligned reads displayed.

GTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATA  >  NC_000913/4196208‑4196500
                                                                                                                                                    |                                                                                                                                                
gTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAGTGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACgcgggcgg                                                                                                                                                 <  1:935492/150‑1 (MQ=255)
gTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACgcgggcgg                                                                                                                                                 <  1:1507897/150‑1 (MQ=255)
  aTTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGgcggcg                                                                                                                                               <  2:1065817/150‑1 (MQ=255)
   ttAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGt                                                                                                                                              <  1:969541/150‑1 (MQ=255)
        cTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGc                                                                                                                                         <  1:1502637/150‑1 (MQ=255)
        cTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAACTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCCGCCCGCGTCaaac                                                                                                                                         <  2:1655794/150‑3 (MQ=255)
         tGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCACCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCa                                                                                                                                        <  1:913149/150‑1 (MQ=255)
         tGATCACGTTCCGCGGATCCCGAATAAACGGGCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGATTGACGCCGCCCGCGTCAAACa                                                                                                                                        <  2:1210541/150‑1 (MQ=255)
           aTCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCAcc                                                                                                                                      <  1:359530/150‑1 (MQ=255)
           aTCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCAcc                                                                                                                                      <  2:428853/150‑1 (MQ=255)
           aTCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCAcc                                                                                                                                      <  2:1642382/150‑1 (MQ=255)
            tCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACca                                                                                                                                     <  1:405082/150‑1 (MQ=255)
            tCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGGTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACca                                                                                                                                     <  2:1263515/150‑1 (MQ=255)
             cAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACcag                                                                                                                                    <  1:1097400/150‑1 (MQ=255)
              aGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACcagc                                                                                                                                   <  2:1589761/150‑1 (MQ=255)
               ggatCCGCGGATCCCGAATAATCGGTCTCAGCCAGTTAACGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGACTGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCTCcagca                                                                                                                                  <  1:288353/147‑1 (MQ=255)
                gTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACcagcag                                                                                                                                 <  2:1183041/150‑1 (MQ=255)
                gTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACcagcag                                                                                                                                 <  1:510427/150‑1 (MQ=255)
                gTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACcagcag                                                                                                                                 <  1:349529/150‑1 (MQ=255)
                    cgcgGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTAAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATc                                                                                                                             <  2:879142/150‑1 (MQ=255)
                     gcgGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCa                                                                                                                            >  1:1188459/1‑150 (MQ=255)
                           cccGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTAt                                                                                                                      >  1:648678/1‑150 (MQ=255)
                           cccGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTAt                                                                                                                      >  2:91716/1‑150 (MQ=255)
                            ccGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATc                                                                                                                     <  1:109676/150‑1 (MQ=255)
                            ccGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATc                                                                                                                     >  2:515148/1‑150 (MQ=255)
                                aTAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTcc                                                                                                                 <  1:1254298/150‑1 (MQ=255)
                                   aaCGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACATCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGc                                                                                                              <  2:1182601/150‑1 (MQ=255)
                                   aaCGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGc                                                                                                              >  1:1437595/1‑150 (MQ=255)
                                         ctcAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAg                                                                                                        >  1:832792/1‑150 (MQ=255)
                                             cccAGTTAAGGCACGCCTACAAGAACTTAGCCCCTCAAATGGGGCATTGAATGTAAAGTACGCGTTAACAGCGCAGAAATCAAGCAGGATGGTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAgcgcg                                                                                                    <  1:872787/149‑1 (MQ=255)
                                              ccAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGACCTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAgcgcgc                                                                                                   >  2:1074545/1‑150 (MQ=255)
                                               cAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCt                                                                                                  >  1:655526/1‑150 (MQ=255)
                                                 gTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTc                                                                                                >  2:872511/1‑150 (MQ=255)
                                                  ttAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTcc                                                                                               <  1:1103746/150‑1 (MQ=255)
                                                      ggCACTCCTACAAGAAATTAGCCCCTCATAGGGGGATTTGAATGTAAATTACGCGTTAACAGCGGAGAACTCAAGCAGGATGTTTGACGCGGGCGGCTTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAgcg                                                                                           <  1:1248306/150‑1 (MQ=255)
                                                        cACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGTCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACACAGCGCGCTTCCAGCGcc                                                                                         >  1:1684180/1‑150 (MQ=255)
                                                            ccGACAAGAACTTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACCGCGCAGACCTCAAGCAGGATGTTGGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGc                                                                                     <  1:132128/150‑1 (MQ=255)
                                                             cGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACACAGCGCGCTTCCAGCGCCTCGcc                                                                                    >  2:142318/1‑150 (MQ=255)
                                                              gACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCg                                                                                   <  1:1622194/150‑1 (MQ=255)
                                                               aCAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGa                                                                                  <  1:840074/150‑1 (MQ=255)
                                                               aCAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGa                                                                                  >  1:934191/1‑150 (MQ=255)
                                                                cAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGAt                                                                                 <  2:1188459/150‑1 (MQ=255)
                                                                cAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGAt                                                                                 >  2:1459268/1‑150 (MQ=255)
                                                                cAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGAt                                                                                 <  2:1374275/150‑1 (MQ=255)
                                                                cAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGAt                                                                                 >  2:1364739/1‑150 (MQ=255)
                                                                 aaGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATg                                                                                >  1:338505/1‑150 (MQ=255)
                                                                  aGAAATTAGCCCCGCAAATGGGTCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGt                                                                               >  2:1239243/1‑150 (MQ=255)
                                                                  aGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCCGCTTATCTTCCAGCACAAAGCGCGCTTCCCGCGCCTCGCCgcggt                                                                               >  1:613053/1‑146 (MQ=255)
                                                                   gAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCTGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTc                                                                              >  1:1488630/1‑150 (MQ=255)
                                                                   gAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTc                                                                              <  2:688281/150‑1 (MQ=255)
                                                                    aaaTTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCa                                                                             <  1:918689/150‑1 (MQ=255)
                                                                     aaTTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAACGCGCGCTTCCAGCGCCTCGCCGATGTCAg                                                                            >  2:729082/1‑150 (MQ=255)
                                                                     aaTTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAg                                                                            >  2:917336/1‑150 (MQ=255)
                                                                     aaTTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAg                                                                            <  2:783511/150‑1 (MQ=255)
                                                                     aaTTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAg                                                                            >  2:82255/1‑150 (MQ=255)
                                                                      aTTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGa                                                                           <  1:470354/150‑1 (MQ=255)
                                                                        tAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATa                                                                         >  2:189770/1‑150 (MQ=255)
                                                                             ccGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAACGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAAcc                                                                    >  1:331783/1‑150 (MQ=255)
                                                                                 aaaTGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGtt                                                                <  1:800141/150‑1 (MQ=255)
                                                                                   aTGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGa                                                              <  1:917336/150‑1 (MQ=255)
                                                                                   aTGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGa                                                              <  1:379506/150‑1 (MQ=255)
                                                                                   aTGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGa                                                              <  1:1593366/150‑1 (MQ=255)
                                                                                    tGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGaa                                                             >  1:703934/1‑150 (MQ=255)
                                                                                     ggggCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGaaa                                                            >  1:1267443/1‑150 (MQ=255)
                                                                                       ggCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACt                                                          >  2:1390726/1‑150 (MQ=255)
                                                                                       ggCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCGCGCCGATGTCAGATAAAACCTGTTGAAACt                                                          <  1:91716/150‑1 (MQ=255)
                                                                                         cATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCa                                                        >  1:6862/1‑150 (MQ=255)
                                                                                          aTTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCaa                                                       <  2:1401109/150‑1 (MQ=255)
                                                                                           ttGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAg                                                      <  1:1482551/150‑1 (MQ=255)
                                                                                            tGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCTGGCGGCGGCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAgg                                                     <  1:268650/150‑1 (MQ=255)
                                                                                              aaTGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAgcgcg                                                                                                    <  2:760734/101‑1 (MQ=255)
                                                                                              aaTGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAgcgcg                                                                                                    >  1:760734/1‑101 (MQ=255)
                                                                                               aTGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTaa                                                  >  2:441813/1‑150 (MQ=255)
                                                                                                  tAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTa                                               >  1:37858/1‑150 (MQ=255)
                                                                                                     aTTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCa                                            >  1:1620158/1‑150 (MQ=255)
                                                                                                          gcgTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATc                                       <  1:1364739/150‑1 (MQ=255)
                                                                                                          gcgTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATc                                       <  1:1459268/150‑1 (MQ=255)
                                                                                                           cgTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCg                                      >  1:1164869/1‑150 (MQ=255)
                                                                                                            gTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATc                                       >  2:374290/1‑148 (MQ=255)
                                                                                                            gTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATc                                       <  1:374290/148‑1 (MQ=255)
                                                                                                                aCAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGc                                 >  1:1670671/1‑150 (MQ=255)
                                                                                                                aCAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGc                                 <  1:515148/150‑1 (MQ=255)
                                                                                                                   gcgcAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGt                              >  2:1682218/1‑150 (MQ=255)
                                                                                                                      cAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTc                           >  2:792209/1‑150 (MQ=255)
                                                                                                                              atgCAGGTTGTTTGACGCGGGCGGCGTCAAGCACCAGAAGACTCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCTCAGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTgg                   <  2:1418814/148‑1 (MQ=255)
                                                                                                                                  aGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATc               <  2:1437595/150‑1 (MQ=255)
                                                                                                                                   tggtgcttGACGCCGCCCGCGTCAAACACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCg              >  1:1642176/7‑150 (MQ=255)
                                                                                                                                   ggctGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCg              <  1:1670059/147‑1 (MQ=255)
                                                                                                                                     gtgcttGACGCCGCCCGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCgta            >  1:1259325/5‑150 (MQ=255)
                                                                                                                                       gcttGACGCCGCCCGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGCTGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCgtagt          >  1:340831/3‑150 (MQ=255)
                                                                                                                                          tGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAAt       >  1:1645208/1‑150 (MQ=255)
                                                                                                                                          tGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAAt       <  1:792209/150‑1 (MQ=255)
                                                                                                                                           gACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATc      >  2:240771/1‑150 (MQ=255)
                                                                                                                                           gACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATc      >  2:1682095/1‑150 (MQ=255)
                                                                                                                                           gACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATc      <  1:1239243/150‑1 (MQ=255)
                                                                                                                                           gACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATc      <  1:1074545/150‑1 (MQ=255)
                                                                                                                                           gACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACGCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATc      >  2:67477/1‑150 (MQ=255)
                                                                                                                                             cgcgGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCa    >  1:1031390/1‑150 (MQ=255)
                                                                                                                                               cgggcggCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATa  >  2:534121/1‑150 (MQ=255)
                                                                                                                                               cgggcggCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGAGGTTTCCAGACTGGAATCGTAGTAATCCATa  >  2:1306101/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATA  >  NC_000913/4196208‑4196500

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