Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,631,642 G→C 5.6% D304H (GAC→CAC) ‡ slt → lytic murein transglycosylase, soluble

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,631,6420GC 5.6% 478.6 / 13.8 143D304H (GAC→CAC) ‡sltlytic murein transglycosylase, soluble
Reads supporting (aligned to +/- strand):  ref base G (69/66);  new base C (4/4);  total (73/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 146 total aligned reads displayed.

AGGATGCTGAGAATGCACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCG  >  NC_000913/4631496‑4631791
                                                                                                                                                  |                                                                                                                                                     
aGGATGCTGAGAATGCACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACg                                                                                                                                                    <  1:650259/150‑1 (MQ=255)
aGGATGCTGAGAATGCACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACg                                                                                                                                                    <  1:1071263/150‑1 (MQ=255)
       tGAGAATGCACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTAt                                                                                                                                             >  1:1089652/1‑150 (MQ=255)
        gagaATGCCCGGCTGATGAGCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATg                                                                                                                                            <  1:1591744/150‑1 (MQ=255)
        gagaATGCACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGTCATCGCGCCAtttcg                                                                                                                                            <  2:84944/150‑4 (MQ=255)
        gagaATGCACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGTCATCGCGCCAtttcg                                                                                                                                            <  2:514709/150‑4 (MQ=255)
        gagaATGCACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATg                                                                                                                                            >  2:236911/1‑150 (MQ=255)
           aaTGCACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGTCATCGCGCCAtttcgcct                                                                                                                                         <  1:24426/150‑7 (MQ=255)
             tGCACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCCCCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTc                                                                                                                                       >  2:150193/1‑150 (MQ=255)
              gCACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCg                                                                                                                                      <  2:1149128/150‑1 (MQ=255)
              gCACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCg                                                                                                                                      >  2:1251920/1‑150 (MQ=255)
              gCACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCg                                                                                                                                      <  1:315898/150‑1 (MQ=255)
              ccACGCCTGAAGAACCCATCCCTCGCCCAGGCGCAGCAGCTAAATGAAGAGCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATCGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCg                                                                                                                                      <  1:872025/149‑1 (MQ=255)
               cACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGc                                                                                                                                     >  2:1602736/1‑150 (MQ=255)
               cACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGGCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGc                                                                                                                                     <  1:95253/150‑1 (MQ=255)
                    cTGATGATACCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGAGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCt                                                                                                                                <  2:1211632/150‑1 (MQ=255)
                         gatCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCCCCGCCGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACttg                                                                                                                           >  2:1632567/1‑149 (MQ=255)
                         gatCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCCCCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTc                                                                                                                           >  2:699901/1‑150 (MQ=255)
                            cccATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCt                                                                                                                        <  1:183429/150‑1 (MQ=255)
                            cccATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCt                                                                                                                        <  2:1921166/150‑1 (MQ=255)
                                tCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTAtt                                                                                                                    <  2:724157/150‑1 (MQ=255)
                                 cGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTg                                                                                                                   <  1:549882/150‑1 (MQ=255)
                                  gCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTgc                                                                                                                  >  1:323839/1‑149 (MQ=255)
                                       cccAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGcc                                                                                                             >  2:294065/1‑150 (MQ=255)
                                        ccAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCg                                                                                                            <  2:477029/150‑1 (MQ=255)
                                           ggCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCgtgt                                                                                                         >  2:1734215/1‑150 (MQ=255)
                                            gcgcAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTa                                                                                                        <  2:1467498/150‑1 (MQ=255)
                                             cgcAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTAc                                                                                                       >  2:788591/1‑150 (MQ=255)
                                             cgcAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTACGCTTATTGAACGCCGTGTAc                                                                                                       >  1:70227/1‑150 (MQ=255)
                                               cagcagCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGa                                                                                                     >  2:667788/1‑150 (MQ=255)
                                                agcagcTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGaa                                                                                                    <  1:1186706/150‑1 (MQ=255)
                                                    gcTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGCACGCCGTGTACGAATGgc                                                                                                >  2:1614419/1‑150 (MQ=255)
                                                    gcTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGgc                                                                                                >  1:1556322/1‑150 (MQ=255)
                                                     cTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGgcg                                                                                               >  2:352561/1‑150 (MQ=255)
                                                       tAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCt                                                                                             >  2:1043650/1‑150 (MQ=255)
                                                        aaTGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCtt                                                                                            <  1:946347/150‑1 (MQ=255)
                                                          tGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTgg                                                                                          >  2:1658280/1‑150 (MQ=255)
                                                          tGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTgg                                                                                          >  1:1106107/1‑150 (MQ=255)
                                                           gAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGc                                                                                         <  2:876770/150‑1 (MQ=255)
                                                              gATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCAcc                                                                                      >  1:1863835/1‑150 (MQ=255)
                                                              gATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCAcc                                                                                      >  2:64345/1‑150 (MQ=255)
                                                               aTCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCg                                                                                     >  1:853709/1‑150 (MQ=255)
                                                                 aaGAATCAGGATCTTCGCTATATCTTAGCCTGGCGTTTGATGCGAAACGATGTCACCGACGAGCAGGCTAAATTGCGAGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCgtc                                                                                   <  2:1246590/149‑3 (MQ=255)
                                                                  aGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCTCGATGAAGCCATTATGCGCTCGCAATCTACGTCGCTTATTGACCGCCTTGTACTAATGGCGCTTGGCACCGGCg                                                                                  >  1:253473/1‑150 (MQ=255)
                                                                  aGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCg                                                                                  <  2:1652138/150‑1 (MQ=255)
                                                                   gATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGa                                                                                 <  1:1454240/150‑1 (MQ=255)
                                                                    aTTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGAt                                                                                <  1:888372/150‑1 (MQ=255)
                                                                    aTTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGAt                                                                                <  1:1320302/150‑1 (MQ=255)
                                                                     ttCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTTGCACCGGCGAtc                                                                               <  2:887733/150‑1 (MQ=255)
                                                                      tCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGAtcg                                                                              <  1:601997/150‑1 (MQ=255)
                                                                        aGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGTCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGAtcgtc                                                                            <  2:1503389/150‑1 (MQ=255)
                                                                         ggAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCGCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGAtcgtcg                                                                           <  2:92265/150‑1 (MQ=255)
                                                                            gcggcgaGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCTCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCgg                                                                        <  2:525429/147‑1 (MQ=255)
                                                                               gcgcGATATCGTCGCCTGGCGTTTGAGGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCTCCGGCGATCGTCGCGGCCt                                                                     >  1:1987655/1‑150 (MQ=255)
                                                                                cgcgATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTg                                                                    <  2:1438960/150‑1 (MQ=255)
                                                                                 gcgATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGa                                                                   >  2:1461926/1‑150 (MQ=255)
                                                                                 gcgATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGa                                                                   >  2:1375818/1‑150 (MQ=255)
                                                                                    atatCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATa                                                                <  2:996169/150‑1 (MQ=255)
                                                                                      atCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATAcc                                                              <  1:294065/150‑1 (MQ=255)
                                                                                       tcgtcgCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACct                                                             <  2:1723647/150‑1 (MQ=255)
                                                                                        cgtcgCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACctg                                                            <  1:1932235/150‑1 (MQ=255)
                                                                                        cgtcgCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACATCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACctg                                                            <  1:1251920/150‑1 (MQ=255)
                                                                                             ccTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACctggctgg                                                       >  2:1913324/1‑150 (MQ=255)
                                                                                             ccTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACctggctgg                                                       <  2:1089652/150‑1 (MQ=255)
                                                                                                ggCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGgcgc                                                    >  1:1582736/1‑150 (MQ=255)
                                                                                                 gCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGcgcg                                                   <  1:855847/150‑1 (MQ=255)
                                                                                                 gCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGcgcg                                                   >  1:592985/1‑150 (MQ=255)
                                                                                                  cGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGt                                                  >  1:361208/1‑150 (MQ=255)
                                                                                                    tttGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCt                                                <  2:1356832/150‑1 (MQ=255)
                                                                                                     ttGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGGGTACGAATGGCGCTTGGCACCAGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTg                                               <  1:556313/150‑1 (MQ=255)
                                                                                                       gATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGcc                                             <  1:381728/150‑1 (MQ=255)
                                                                                                       gATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGcc                                             <  2:948814/150‑1 (MQ=255)
                                                                                                        aTGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCt                                            >  2:869924/1‑150 (MQ=255)
                                                                                                           ggCAACGATGTCTCCGACGAGCAGGCGAACTGGCGCGATGACCCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCGATg                                         <  2:1878350/150‑1 (MQ=255)
                                                                                                           ggCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATg                                         >  1:455635/1‑150 (MQ=255)
                                                                                                           ggCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATg                                         >  2:502498/1‑150 (MQ=255)
                                                                                                            gCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATgg                                        >  2:383078/1‑150 (MQ=255)
                                                                                                                 gATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCg                                   <  2:1863835/150‑1 (MQ=255)
                                                                                                                     tCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaag                               >  2:742424/1‑150 (MQ=255)
                                                                                                                      cacaGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCGTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGACGCGGCCTGAATACCAGGCTGGCGCGTCTGCCAATGGAAGCgaaaga                              >  2:1263441/1‑150 (MQ=255)
                                                                                                                        ccGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaagaga                            <  1:1168055/150‑1 (MQ=255)
                                                                                                                             gaGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGAt                       <  2:846874/150‑1 (MQ=255)
                                                                                                                              aGCAGGCGAAATGGCGCGATGCCGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATg                      >  1:1651132/1‑150 (MQ=255)
                                                                                                                              aGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATg                      <  1:502498/150‑1 (MQ=255)
                                                                                                                              aGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATg                      >  1:1824364/1‑150 (MQ=255)
                                                                                                                                cAGGCGAAATGGCGCGATGACGCCATTATGCGCGCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGaa                    <  1:1380228/150‑1 (MQ=255)
                                                                                                                                 aGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCGGCCTATGGAAGCGAAAGAGAAAGATGAAt                   >  2:1428068/1‑150 (MQ=255)
                                                                                                                                  ggCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATg                  >  2:829163/1‑150 (MQ=255)
                                                                                                                                   gcataaTGGCGTCATCGCGCCATTATGCGCTCGCAATCTACTTCGCGTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATgg                 >  1:1980991/5‑150 (MQ=255)
                                                                                                                                    cGAAATGGCGCGATGCCGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGGACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGc                >  1:1611066/1‑150 (MQ=255)
                                                                                                                                    cGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGc                >  2:303576/1‑150 (MQ=255)
                                                                                                                                      taaTGGCGTCATCGCGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGt              >  1:823982/2‑150 (MQ=255)
                                                                                                                                      aaaTGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGt              <  2:894898/150‑1 (MQ=255)
                                                                                                                                         tGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTAc           >  1:147008/1‑150 (MQ=255)
                                                                                                                                           gcgcGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTg         >  2:1979879/1‑150 (MQ=255)
                                                                                                                                            cgcgATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTgg        <  2:455635/150‑1 (MQ=255)
                                                                                                                                              cgATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCa      <  1:788591/150‑1 (MQ=255)
                                                                                                                                               gATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAg     <  2:853709/150‑1 (MQ=255)
                                                                                                                                                  gACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGTGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCg  >  1:816028/1‑150 (MQ=255)
                                                                                                                                                  gACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCg  <  1:1270462/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
AGGATGCTGAGAATGCACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCG  >  NC_000913/4631496‑4631791

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: