Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,148,049 C→G 8.4% G90A (GGC→GCC)  yijO ← AraC family putative transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,148,0490CG 8.4% 360.1 / 14.4 114G90A (GGC→GCC) yijOAraC family putative transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (42/59);  new base G (5/5);  total (49/65)
Fisher's exact test for biased strand distribution p-value = 7.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 117 total aligned reads displayed.

GCTGCTCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGTGGA  >  NC_000913/4147900‑4148198
                                                                                                                                                     |                                                                                                                                                     
gctgctCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGc                                                                                                                                                       >  1:830440/1‑150 (MQ=255)
   gctCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCaa                                                                                                                                                    <  2:830440/150‑1 (MQ=255)
    ctCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCaaa                                                                                                                                                   >  1:482629/1‑150 (MQ=255)
     tCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAAc                                                                                                                                                  <  2:609080/150‑1 (MQ=255)
      ccTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACa                                                                                                                                                 <  2:521634/150‑1 (MQ=255)
        tGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATTTTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAg                                                                                                                                                <  2:524982/149‑1 (MQ=255)
          cGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCAc                                                                                                                                             <  1:613396/150‑1 (MQ=255)
           ggCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACg                                                                                                                                            <  2:1254663/150‑1 (MQ=255)
           ggCTGCATTTGCAATTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGGCAAACAGCaac                                                                                                                                            <  2:787247/150‑3 (MQ=255)
            gCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGc                                                                                                                                           <  1:996505/150‑1 (MQ=255)
              tgcatttgcattTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGGCAAACAGCaactcg                                                                                                                                         <  2:777705/150‑6 (MQ=255)
              tgcatttgcattTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTa                                                                                                                                         <  1:1405298/150‑1 (MQ=255)
              ggcatttgcatttCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTa                                                                                                                                         <  1:686541/149‑1 (MQ=255)
               gcatttgcatttCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTaa                                                                                                                                        >  1:1295622/1‑150 (MQ=255)
                catttgcatttCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTaaa                                                                                                                                       >  1:26701/1‑150 (MQ=255)
                 atttgcatttCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAg                                                                                                                                      <  1:189808/150‑1 (MQ=255)
                  tttgcatttCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGAGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAgg                                                                                                                                     >  1:10926/1‑150 (MQ=255)
                   ttgcatttCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGgt                                                                                                                                    <  2:104011/150‑1 (MQ=255)
                   ttgcatttCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGgt                                                                                                                                    <  2:1123487/150‑1 (MQ=255)
                   ttgcatttCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGgt                                                                                                                                    <  2:1596329/150‑1 (MQ=255)
                      catttCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGgtagt                                                                                                                                 <  2:85577/150‑1 (MQ=255)
                      catttCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGgtagt                                                                                                                                 >  2:1267210/1‑150 (MQ=255)
                           cATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCgg                                                                                                                            >  2:1409866/1‑150 (MQ=255)
                            aTTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCggg                                                                                                                           <  1:1450317/150‑1 (MQ=255)
                             ttGAGCGTTTGTAGCAGAAAAGAACCTATGCTTGGGCCGCGTCGGGCGACGGGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTGGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCgggg                                                                                                                          <  1:899617/150‑1 (MQ=255)
                               gAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGc                                                                                                                        <  2:1169030/150‑1 (MQ=255)
                                  cGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCTTCCCATTGCACGTCGCTGACGCCGAGTTGCTGTTTGCCAAACAGCTCGCTAAAGGTAGTAGCGGGGGCtcg                                                                                                                     >  2:1229277/1‑148 (MQ=255)
                                  cGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTg                                                                                                                     <  1:1409866/150‑1 (MQ=255)
                                  cGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTg                                                                                                                     <  1:869244/150‑1 (MQ=255)
                                  cGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTg                                                                                                                     <  2:1686917/150‑1 (MQ=255)
                                   ggttGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGc                                                                                                                    <  2:1114198/148‑1 (MQ=255)
                                   gTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGc                                                                                                                    <  2:1347431/150‑1 (MQ=255)
                                      tGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCAt                                                                                                                 <  1:451892/150‑1 (MQ=255)
                                         aGCAGAAAAGAACCTTTGGGTTGGCCGCGTCGGGCGACGTGTTGCTTCGCAAGATTTTTTTTTTGTTTGCCATCCCATGGCACGACGCTGAAGCCGAGTTGCTGTTTGACACACTGCACGCTAAAGGTAGTAGCGGGTGCTTGCCATtgt                                                                                                              <  2:261417/150‑1 (MQ=255)
                                         aGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATtgt                                                                                                              >  1:1581978/1‑150 (MQ=255)
                                          gCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATtgtg                                                                                                             <  1:1639272/150‑1 (MQ=255)
                                          gCAGAAAAGAACCTATGCGTGGGCCGCGTCAGGCGCCGAGTTTCTTCGCCAGATATTGATATTGTTTGCCATCCCATTGCACGACGCTTAAGCCGAGTTGCTGGTTGACACACAGCGAGGTATAGGTCGTAGCGGGGGCTTGGCAtggat                                                                                                             >  1:261549/1‑146 (MQ=255)
                                           cAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTgg                                                                                                            <  1:1392118/150‑1 (MQ=255)
                                             gAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGaa                                                                                                          >  2:488572/1‑150 (MQ=255)
                                               aaaGAACCTATGAGTGGGCCGCGTCGGTCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCAGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCAAATGGTAGTAGCGGGGTCTTGACATTGTGGaaaa                                                                                                        >  1:1064958/1‑150 (MQ=255)
                                                   aaCCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTcc                                                                                                    <  1:902097/150‑1 (MQ=255)
                                                     ccTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATCTTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAg                                                                                                  <  1:1208514/150‑1 (MQ=255)
                                                      cTATGCGTGGGCCGTGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGc                                                                                                 >  2:671345/1‑150 (MQ=255)
                                                      cTATGCGTGGGCCGTGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGc                                                                                                 >  2:44180/1‑150 (MQ=255)
                                                      cTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGc                                                                                                 <  2:1106850/150‑1 (MQ=255)
                                                       tATGCGTGGGCCGCGTCGGGCGACGTGTTTCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGACAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGcc                                                                                                >  2:162130/1‑150 (MQ=255)
                                                       tATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGTGTTGCTGTTTGGCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGcc                                                                                                <  2:139701/150‑1 (MQ=255)
                                                       tATGCGTGGGCCGCGGCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGcc                                                                                                <  1:1510862/150‑1 (MQ=255)
                                                                 ccGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTAGCCAAACAGCACGCAAAAGGTAGTAGCGGGGGCTTGCCATTCGGGAAAATTCCAGCCACCAGCAGGa                                                                                      >  1:387711/1‑150 (MQ=255)
                                                                    cgTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACg                                                                                   >  2:110906/1‑150 (MQ=255)
                                                                      tCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTa                                                                                 <  1:1001846/150‑1 (MQ=255)
                                                                      gcGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTa                                                                                 <  2:1663536/149‑1 (MQ=255)
                                                                        gggCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTAGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACa                                                                               >  2:1406217/1‑150 (MQ=255)
                                                                            gACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGcac                                                                           <  2:1707581/150‑1 (MQ=255)
                                                                                   tGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCg                                                                    >  1:416896/1‑150 (MQ=255)
                                                                                     cTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCggg                                                                  >  2:916318/1‑150 (MQ=255)
                                                                                     cTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCggg                                                                  >  1:584105/1‑150 (MQ=255)
                                                                                     cTTCGCCAGATGTTGATAATGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCggg                                                                  <  1:636623/150‑1 (MQ=255)
                                                                                     cTTCCCCAGACTTTGATATTTTTAGCCATCCCATTGCAAGCCCCTGAAGCCGAGTTGCTGTTTGCCAAACATCAAGCTTAATGTACTTGCGGGGGCTTGCCATTGTTGTAAATTCCAGCCACCAGGATGAACGTTCAGTACATCGCCggg                                                                  <  2:885355/150‑1 (MQ=255)
                                                                                        cGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAAc                                                               >  1:1583236/1‑150 (MQ=255)
                                                                                         gCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACt                                                              >  1:979709/1‑150 (MQ=255)
                                                                                         gCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACt                                                              <  2:1581978/150‑1 (MQ=255)
                                                                                         gCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACt                                                              >  2:1436653/1‑150 (MQ=255)
                                                                                         gCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCAGTTTGCCAAACAGCACGCTAAAGGTACTAGCGGGGGCTTGCCATAGTCGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACt                                                              >  1:1575320/1‑150 (MQ=255)
                                                                                           cAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTaa                                                            <  2:482629/150‑1 (MQ=255)
                                                                                           cAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCGAAAGGTAGTAGCGGGGGCTTGCCATGGTGGAAAATTCCAGCCCCCAGCAGGAACGTACAGCACATCGCCGGGAACGaa                                                            >  1:311056/1‑150 (MQ=255)
                                                                                            aGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAAc                                                           <  1:1604307/150‑1 (MQ=255)
                                                                                            aGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAAc                                                           >  1:1500408/1‑150 (MQ=255)
                                                                                             gATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGGCAAACAGCAAGCTAAAGGTAGTAGCGGGGGCTTGCGATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACg                                                          >  1:335103/1‑150 (MQ=255)
                                                                                              aTTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGt                                                         <  2:1295622/150‑1 (MQ=255)
                                                                                                tttGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTg                                                       <  1:1267210/150‑1 (MQ=255)
                                                                                                 ttGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGc                                                      >  2:1385375/1‑150 (MQ=255)
                                                                                                    atatTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTcc                                                   <  1:828/150‑1 (MQ=255)
                                                                                                      atTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAg                                                 <  2:1060151/150‑1 (MQ=255)
                                                                                                           ttGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCg                                            <  2:1686801/150‑1 (MQ=255)
                                                                                                              ccATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTa                                         >  2:1234243/1‑150 (MQ=255)
                                                                                                              ccATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTa                                         >  1:1525493/1‑150 (MQ=255)
                                                                                                                aTCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATc                                       <  2:10926/150‑1 (MQ=255)
                                                                                                                 tCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATca                                      <  1:488572/150‑1 (MQ=255)
                                                                                                                   ccATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATcaac                                    >  1:287423/1‑150 (MQ=255)
                                                                                                                    cATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATcaaca                                   >  1:875399/1‑150 (MQ=255)
                                                                                                                       tGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAAc                                <  2:1498712/150‑1 (MQ=255)
                                                                                                                        gCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACt                               >  2:237251/1‑150 (MQ=255)
                                                                                                                            gacgCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCAcc                           <  1:480484/150‑1 (MQ=255)
                                                                                                                            gacgCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCAcc                           >  2:1547770/1‑150 (MQ=255)
                                                                                                                            cacgCTGAAGCCGAGTAGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCAcc                           <  2:627609/149‑1 (MQ=255)
                                                                                                                             acgCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCt                          <  2:1525493/150‑1 (MQ=255)
                                                                                                                                cTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTcc                       >  2:1202626/1‑150 (MQ=255)
                                                                                                                                cTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTcc                       <  2:1458918/150‑1 (MQ=255)
                                                                                                                                 tGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCa                      <  1:110906/150‑1 (MQ=255)
                                                                                                                                  gAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAg                     >  1:537248/1‑150 (MQ=255)
                                                                                                                                    aGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCGTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGca                   <  1:831490/150‑1 (MQ=255)
                                                                                                                                      ttagcgtGCTGTTTGGCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGcaca                 <  2:269351/144‑1 (MQ=255)
                                                                                                                                        agcgtGCTGTTTGGCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAAt               >  1:1523274/5‑150 (MQ=255)
                                                                                                                                         gcgtGCTGTTTGGCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAAtt              >  1:1509790/4‑150 (MQ=255)
                                                                                                                                         gcgtGCTGTTTGGCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAAtt              >  1:537111/4‑150 (MQ=255)
                                                                                                                                         gcgtGCTGTTTGGCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAAtt              >  1:1199772/4‑150 (MQ=255)
                                                                                                                                            tGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTcc           >  1:1552809/1‑150 (MQ=255)
                                                                                                                                              cTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAg         >  1:1528612/1‑150 (MQ=255)
                                                                                                                                                     ccAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGTgca  >  2:197576/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCTGCTCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGTGGA  >  NC_000913/4147900‑4148198

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