Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,388,931 T→G 7.4% D147A (GAC→GCC)  mscM ← mechanosensitive channel protein, miniconductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,388,9310TG 7.4% 328.3 / 2.9 142D147A (GAC→GCC) mscMmechanosensitive channel protein, miniconductance
Reads supporting (aligned to +/- strand):  ref base T (52/75);  new base G (2/11);  total (56/86)
Fisher's exact test for biased strand distribution p-value = 8.12e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 142 total aligned reads displayed.

GACAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGA  >  NC_000913/4388784‑4389079
                                                                                                                                                   |                                                                                                                                                    
gACAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCg                                                                                                                                                    >  1:675327/1‑150 (MQ=255)
 aCAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCTTAATACGTGCAGAGGCAGACTGCAACGCGAAAGTTTGTGCCTGATTGAGCGGAGTATTGCCTGTAAGCGTTCCCAGGCGGCGCTCGATATCATTTAACTGACGGCGGGCGTCgg                                                                                                                                                   <  1:769031/150‑1 (MQ=255)
 aCAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGTGGGCGTCgg                                                                                                                                                   <  1:1759241/150‑1 (MQ=255)
 aCAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCgg                                                                                                                                                   <  1:1823841/150‑1 (MQ=255)
  cAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGt                                                                                                                                                  <  1:1746410/150‑1 (MQ=255)
    gCTGCGCCAGCTCCAGTTCATCAACGAGCGCCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGttt                                                                                                                                                >  1:379554/1‑150 (MQ=255)
      tGCGCCAGCTCCAGTTCATCAACGAGCGCCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGc                                                                                                                                              <  2:339540/150‑1 (MQ=255)
            cgCTCCAGTTCATCAACGAGCGCCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGc                                                                                                                                        <  2:79628/149‑1 (MQ=255)
             gCTCCAGTTCATCAACGAGCGCCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCg                                                                                                                                       >  2:203566/1‑150 (MQ=255)
              cTCCAGTTCATCAACGAGCGCCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGGCGGTTTGCTGTTGCgg                                                                                                                                      >  2:129613/1‑150 (MQ=255)
                ccAGTTCATCAACGAGCGCCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCa                                                                                                                                    >  2:759022/1‑150 (MQ=255)
                    ttCATCAACGAGCGCCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTg                                                                                                                                <  2:376281/150‑1 (MQ=255)
                     tcatcaACGAGCGCCTTAAGACGTGCAGAGTCAGCCTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGa                                                                                                                               >  2:851603/1‑150 (MQ=255)
                     tcatcaACGAGCGCCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAGTTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGa                                                                                                                               <  1:257481/150‑1 (MQ=255)
                      catcaACGAGCGCCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGAt                                                                                                                              >  1:1702503/1‑150 (MQ=255)
                      aatcaACGAGCGCCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGAt                                                                                                                              <  2:561090/149‑1 (MQ=255)
                        tcaACGAGCGCCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGACTGATTGAGCGGAGTATTGCCGGTACGCGTTCCCAGGCGGCGCTGGAGCTCATTTAACTGACGGCCGGCGACGGATTGCTGTTGCGGCAGTTGAttg                                                                                                                            >  2:115012/1‑149 (MQ=255)
                         caACGAGCGCCTTAAGACGTGCAGAGTCAGCCTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCAGTTAACTGACGGCGGGCGACGGTTTGCTGTTTCGGCAGTTGATTcc                                                                                                                           >  2:1202575/1‑149 (MQ=255)
                            cGAGCGCCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCg                                                                                                                        >  2:1777436/1‑150 (MQ=255)
                                   cTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGc                                                                                                                 <  2:1449809/150‑1 (MQ=255)
                                    ttAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCa                                                                                                                <  2:929561/150‑1 (MQ=255)
                                       aGACGTGCAGAGGCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAAtc                                                                                                             <  1:1424471/150‑1 (MQ=255)
                                        gACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAAtct                                                                                                            >  2:1180398/1‑150 (MQ=255)
                                        gACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAAtct                                                                                                            >  1:743917/1‑150 (MQ=255)
                                          cGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTcg                                                                                                          >  1:86045/1‑150 (MQ=255)
                                            tGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTcgcg                                                                                                        <  2:86045/150‑1 (MQ=255)
                                            tGCAGAGGCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTcgcg                                                                                                        <  2:1513105/150‑1 (MQ=255)
                                                gaggccgACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCGCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGgcg                                                                                                    <  1:1467051/144‑1 (MQ=255)
                                                gagTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCTGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGgcg                                                                                                    <  2:1063984/150‑1 (MQ=255)
                                                 agTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGgcgc                                                                                                   >  2:554600/1‑150 (MQ=255)
                                                 agTCAGACTGCAACGCGAAATTTTGGGCCTGATTGAGCGGAGTATTGCCGGTACGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGgcgc                                                                                                   <  2:1150442/150‑1 (MQ=255)
                                                  ggcAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCTGCGCTCGATCTCATTTAACTGACGGCGGGCGGCGGTATGCTGGTGCGGCAGTTGATTCCGCGAATCGTCAATATAGCGGgcgcc                                                                                                  >  2:267757/3‑149 (MQ=255)
                                                   tCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGgcgcgc                                                                                                 >  2:1758490/1‑150 (MQ=255)
                                                       aCTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCt                                                                                             >  1:1492746/1‑150 (MQ=255)
                                                        cttcAACGCGAAAGTTTGTGCCTGATTGAGCGGAGTAGTGCCGGTAAGCGTTCCCAGGCGGCGCTCGAGCTCATTTCACTGACGCCGGGCGGCGGTTTGCTGTTGCGGGAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTg                                                                                            <  1:334730/147‑1 (MQ=255)
                                                        cTGCAACGCGAACTGATGTGCCTGCTTGAGCGGAGTATTGCCGGGAACCGCTCCCAGGCGGCGCTCGATCTCATCTAACTGACGGCGGGCGGCGGTGTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGTGCGCGCTCCTg                                                                                            <  1:265067/150‑1 (MQ=255)
                                                           cAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGctc                                                                                         >  1:1153973/1‑150 (MQ=255)
                                                           cAACGCGAAATTTTGTGCCGGTTTGAGCGGAGTATTGCCGGTAAGCGGTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGctc                                                                                         <  2:908697/150‑1 (MQ=255)
                                                            aaCGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGctct                                                                                        <  2:379554/150‑1 (MQ=255)
                                                            aaCGCGAAATTTTGTGCCGGATTGAGCGGAGTATTGCCGGTAAGCGTGCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTGTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGctct                                                                                        <  1:344394/150‑1 (MQ=255)
                                                               gcgAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGc                                                                                     >  2:945805/1‑150 (MQ=255)
                                                                cgAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTTATTCAGCGAATCGGCAATCTCTCGGGCGCGCTCCTGCTCTTGCt                                                                                    >  1:1175692/1‑150 (MQ=255)
                                                                cgAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCt                                                                                    >  2:1721970/1‑150 (MQ=255)
                                                                cgAAATTTTGAGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGTTTTGCTGTTGAGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCt                                                                                    >  1:450787/1‑150 (MQ=255)
                                                                 gAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCACTCTCGCGGGCGCGCTCCTGCTCTTGCTg                                                                                   >  1:325474/1‑150 (MQ=255)
                                                                   aaTTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTggg                                                                                 >  1:465436/1‑150 (MQ=255)
                                                                   aaTTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTggg                                                                                 <  1:756550/150‑1 (MQ=255)
                                                                   aaTTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTggg                                                                                 <  2:1613647/150‑1 (MQ=255)
                                                                    aTTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGc                                                                                <  1:1014831/150‑1 (MQ=255)
                                                                    aTTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGc                                                                                <  2:591096/150‑1 (MQ=255)
                                                                     ttttGTGCCTGATTGAGCGGGGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGcc                                                                               <  1:1495035/150‑1 (MQ=255)
                                                                     ttttGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGGCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGcc                                                                               <  2:1630386/150‑1 (MQ=255)
                                                                     ttttGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGGTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGGCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGcc                                                                               <  2:487102/150‑1 (MQ=255)
                                                                         gtgCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGAc                                                                           <  2:35902/150‑1 (MQ=255)
                                                                          tgCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACg                                                                          <  2:1469783/150‑1 (MQ=255)
                                                                          tgCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACg                                                                          <  2:1376079/150‑1 (MQ=255)
                                                                             cTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGCTCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCt                                                                       >  2:233331/1‑150 (MQ=255)
                                                                              tGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCtt                                                                      >  2:420384/1‑150 (MQ=255)
                                                                              tGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCtt                                                                      <  2:399054/150‑1 (MQ=255)
                                                                               gagtGAGCGGAGTACTGCCGGTCAGCGATCCCCGGCGGCGCGCGATCTCATTTCACTGACGGCGGGCGTTGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGCCGCTCTCCTGCTCTTGCTGGGCCTGACGGCttt                                                                     <  1:274940/147‑1 (MQ=255)
                                                                               gATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCttt                                                                     <  1:537261/150‑1 (MQ=255)
                                                                               gATTGAGCGGAGGATTGCCGGTATGCGTTCCCAGGCGGCGCGCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCttt                                                                     <  1:1202575/150‑1 (MQ=255)
                                                                                agtGAGCGGAGTAATGCCGGTAAGCGGTCCCAGGCGGCGCTCGATCTCATATAACTGACGGCGGGCGGCGGTTTGCGGTTGCGGCAGTTGATTCAGCGAATCGTCAATATCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCtttt                                                                    <  1:1195151/148‑1 (MQ=255)
                                                                                  tGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGCTCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGGTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTatat                                                                  >  2:105888/1‑150 (MQ=255)
                                                                                   gAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGCCGGCTTTTATc                                                                 >  2:904274/1‑150 (MQ=255)
                                                                                    aGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATcc                                                                <  2:332279/150‑1 (MQ=255)
                                                                                     gCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCa                                                               >  1:735132/1‑150 (MQ=255)
                                                                                      cGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAg                                                              <  1:51198/150‑1 (MQ=255)
                                                                                           tATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACt                                                         >  1:853222/1‑150 (MQ=255)
                                                                                           tATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACt                                                         <  1:556268/150‑1 (MQ=255)
                                                                                               gCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGgct                                                     >  1:1802458/1‑150 (MQ=255)
                                                                                                ccGGTAAGCGTTGCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGgctg                                                    <  1:162337/150‑1 (MQ=255)
                                                                                                 cGGTAAGCGTTCCCAGGCGGCGCTCGATCTAAATTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGgctgc                                                   <  2:475687/150‑1 (MQ=255)
                                                                                                   gTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCGGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGctgctg                                                 <  2:672755/150‑1 (MQ=255)
                                                                                                    tAAGCGTTCCCCGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGGCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGa                                                <  1:1060756/150‑1 (MQ=255)
                                                                                                     aaGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGCCGGCTTTTATCCAGCAACTGGCTGCTGAc                                               >  1:1368899/1‑150 (MQ=255)
                                                                                                     aaGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGAc                                               >  1:537707/1‑150 (MQ=255)
                                                                                                      aGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTAGGTTTGCTGTTTCGGCAGTTGATTCAGCGAATCGGCAATCTAGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGAAGGCTTTTATCGAGCAACTGGCTGCTGAcc                                              <  1:1339591/150‑1 (MQ=255)
                                                                                                      aGCGTTCCCAGGCGGCGCTCGATCTAATTTAACTGACGGCGTGCGTCTGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGAcc                                              <  2:1217966/150‑1 (MQ=255)
                                                                                                      aGCGTGCCCAGGCGGCGCTCGATCTCATTTAAGTGACGGCGGGCGTCGGGTTGCTGTGGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCGTTTATCCAGCAACTGGCTGCTGAcc                                              <  1:1381177/150‑1 (MQ=255)
                                                                                                       gCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCt                                             >  1:1686989/1‑150 (MQ=255)
                                                                                                        cGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTg                                            >  1:1249991/1‑150 (MQ=255)
                                                                                                           tCCCAGGCGGCGCTCGATCTCATTTAACGGACGGCGGGCGGCGGTTTGCTGGTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGgag                                         <  2:1125957/150‑1 (MQ=255)
                                                                                                             ccAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGaa                                       >  1:1088846/1‑150 (MQ=255)
                                                                                                                  cggcgCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTcc                                  <  2:470447/150‑1 (MQ=255)
                                                                                                                  cggcgCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGGTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTcc                                  <  1:696889/150‑1 (MQ=255)
                                                                                                                  cggcgCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGGTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTcc                                  <  1:1552962/150‑1 (MQ=255)
                                                                                                                   ggcgCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCt                                 <  1:1666650/150‑1 (MQ=255)
                                                                                                                   gcctcTCGATCTCATTTAACTGTCGGAGGGCGTCGGTTTGCTGTTGCTGCAGTTGATTCAGCGAATGGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACTGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCt                                 <  1:65370/146‑1 (MQ=255)
                                                                                                                    gcgcTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTg                                >  2:244653/1‑150 (MQ=255)
                                                                                                                    gcgcTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTg                                >  1:1807737/1‑150 (MQ=255)
                                                                                                                    gcgcTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTg                                <  2:685114/150‑1 (MQ=255)
                                                                                                                      gcTCGATCTCATTTAACTGACTGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGAt                              <  2:1189532/150‑1 (MQ=255)
                                                                                                                      gcTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGACAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGAt                              <  2:913644/150‑1 (MQ=255)
                                                                                                                       cTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGAtt                             <  1:796641/150‑1 (MQ=255)
                                                                                                                         cGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCa                           >  2:800689/1‑150 (MQ=255)
                                                                                                                           aTCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTGGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAgc                         <  1:173837/150‑1 (MQ=255)
                                                                                                                                      cTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGa              <  2:675327/150‑1 (MQ=255)
                                                                                                                                                 cgTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCg   >  2:1238892/1‑150 (MQ=255)
                                                                                                                                                 cgTCGGGTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCg   <  1:1758490/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GACAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGA  >  NC_000913/4388784‑4389079

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: