Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,631,700 T→C 100% L323P (CTT→CCT)  slt → lytic murein transglycosylase, soluble

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,631,7000TC98.3% 422.1 / ‑3.0 129L323P (CTT→CCT) sltlytic murein transglycosylase, soluble
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (52/74);  total (53/76)
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.71e-01

Only 100 of 129 total aligned reads displayed.

GGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGATTT  >  NC_000913/4631603‑4631796
                                                                                                 |                                                                                                
ggCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTg                                                                                                <  2:225120/100‑1 (MQ=255)
ggCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTg                                                                                                <  2:440376/100‑1 (MQ=255)
ggCAACGATGTCAACGACGAGCATGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAAGCTACTTCCCTTATTGAACGCCGTGTACGAATGGCGCCTg                                                                                                <  2:726234/100‑1 (MQ=255)
 gCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTgg                                                                                               <  2:151129/100‑1 (MQ=255)
    aCGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCAc                                                                                            <  1:1429108/100‑1 (MQ=255)
     cGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCAcc                                                                                           >  1:460167/1‑100 (MQ=255)
     cGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCAcc                                                                                           <  1:1867098/100‑1 (MQ=255)
        tGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGc                                                                                        <  2:1696332/100‑1 (MQ=255)
          tCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGa                                                                                      >  2:2991510/1‑100 (MQ=255)
           cACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGAt                                                                                     <  1:1996053/100‑1 (MQ=255)
           cACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGAt                                                                                     >  2:2793970/1‑100 (MQ=255)
           cACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGAt                                                                                     <  1:2224973/100‑1 (MQ=255)
            aCCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGAt                                                                                     <  1:2941224/99‑1 (MQ=255)
            aCCGACGAGCAGGCGAAATGGCGCGATAAAAAAATAAAACGATCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGAtc                                                                                    <  2:1444549/100‑1 (MQ=255)
             ccGACGAGCATGCGAAATGGCGCGATGACGACATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGGACGAATGGCGCCTGGCACCGGCGAtcg                                                                                   >  2:1883794/1‑100 (MQ=255)
             ccGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGAtcg                                                                                   <  1:2342661/100‑1 (MQ=255)
               gacgaGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGAtcgtc                                                                                 >  2:2342524/1‑100 (MQ=255)
                acgaGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGAtcgtcg                                                                                <  2:2443886/100‑1 (MQ=255)
                   aGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCgg                                                                             >  2:1467404/1‑100 (MQ=255)
                   aGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCgg                                                                             >  1:2444660/1‑100 (MQ=255)
                   aGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCCCCTGGCACCGGCGATCGTcgcg                                                                              >  2:736599/1‑99 (MQ=255)
                     cAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGTCGCCTGGCACCGGCGATCGTCGCGGcc                                                                           >  1:2280683/1‑100 (MQ=255)
                       ggCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCgg                                                                             >  2:1036070/1‑96 (MQ=255)
                       ggCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTg                                                                         >  1:795020/1‑100 (MQ=255)
                       ggCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTg                                                                         <  1:2336287/100‑1 (MQ=255)
                       ggCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCACCTGGCACCGGCGATCGTCGCGGCCTg                                                                         >  1:2508264/1‑100 (MQ=255)
                       ggCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGAGCCTGGCACCGGCGATCGTCGCGGCCTg                                                                         >  2:2547966/1‑100 (MQ=255)
                        gCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTg                                                                         <  2:1490252/99‑1 (MQ=255)
                        gCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGa                                                                        <  2:686795/100‑1 (MQ=255)
                         cGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGaa                                                                       >  1:1453882/1‑100 (MQ=255)
                         cGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGaa                                                                       <  1:975083/100‑1 (MQ=255)
                         cGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGaa                                                                       >  1:622597/1‑100 (MQ=255)
                            aaTGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATAc                                                                    <  2:2295531/100‑1 (MQ=255)
                             aTGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATAcc                                                                   >  1:866746/1‑100 (MQ=255)
                             aTGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATAcc                                                                   <  2:1793928/100‑1 (MQ=255)
                               ggCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACctg                                                                 <  1:2794002/100‑1 (MQ=255)
                               ggCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACctg                                                                 <  1:2677940/100‑1 (MQ=255)
                                gcgcGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACctgg                                                                <  1:1521429/100‑1 (MQ=255)
                                 cgcgATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACctggc                                                               >  1:1916310/1‑100 (MQ=255)
                                 cgcgATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACctggc                                                               <  2:2844858/100‑1 (MQ=255)
                                    gATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACctggctgg                                                            <  1:1449403/100‑1 (MQ=255)
                                      tGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGgcg                                                          >  1:2916423/1‑100 (MQ=255)
                                      tGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGgcg                                                          <  1:1611470/100‑1 (MQ=255)
                                        aCGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGcgcg                                                        >  2:1048867/1‑100 (MQ=255)
                                        aCGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGAGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGgcgc                                                         >  1:1490972/1‑99 (MQ=255)
                                        aCGCCATTATGCGCGCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGcgcg                                                        <  2:2069552/100‑1 (MQ=255)
                                         cGCCCTTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGACCGAAGCGGCCTGAATACCTGGCTGGcgcg                                                        >  2:2076617/1‑99 (MQ=255)
                                          gCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGACGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTc                                                      >  1:1226957/1‑100 (MQ=255)
                                           ccATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCt                                                     <  2:2102568/100‑1 (MQ=255)
                                           ccATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGAGCGTCGCGGCCTGAATACCTGGCTGGCGCGTCt                                                     <  1:2885989/100‑1 (MQ=255)
                                            cATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTg                                                    <  2:51868/100‑1 (MQ=255)
                                             aTTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGc                                                   <  1:1013657/100‑1 (MQ=255)
                                             aTTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGACTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGc                                                   >  2:1844113/1‑100 (MQ=255)
                                              ttATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGcc                                                  <  2:1236479/100‑1 (MQ=255)
                                              ttATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGcc                                                  >  1:421307/1‑100 (MQ=255)
                                              ttATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGCGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGcc                                                  <  2:2281644/100‑1 (MQ=255)
                                                   cgcTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATgg                                             <  1:2342551/100‑1 (MQ=255)
                                                    gcTCGCAATCTACTTCGCTTATTGAACGCTGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGa                                            <  1:809825/100‑1 (MQ=255)
                                                    gcTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGa                                            <  1:2641715/100‑1 (MQ=255)
                                                      tCGTAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAg                                          <  1:328886/100‑1 (MQ=255)
                                                       cGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGACTGGCGCGTCTGCCTATGGAAGc                                         >  2:1333720/1‑100 (MQ=255)
                                                        gCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCg                                        <  1:2544022/100‑1 (MQ=255)
                                                          aaTCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaa                                      <  2:866742/100‑1 (MQ=255)
                                                           aTCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaa                                     <  1:450991/100‑1 (MQ=255)
                                                            tCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaag                                    <  2:1156014/100‑1 (MQ=255)
                                                             cTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaaga                                   <  2:136223/100‑1 (MQ=255)
                                                             cTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACATGGCTGGCGCGTCTGCCTATGGAAGCgaaaga                                   >  1:213220/1‑100 (MQ=255)
                                                             cTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGAATGAAAAACTGGCAGGCGCGTCTGCCTATGGAAGCgaaaga                                   >  1:2647616/1‑100 (MQ=255)
                                                              tACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaagag                                  >  2:206963/1‑100 (MQ=255)
                                                              tACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaaga                                   >  2:2097060/1‑99 (MQ=255)
                                                              tACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCCGCCTATGGAAGCgaaagag                                  >  1:1053521/1‑100 (MQ=255)
                                                               aCTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaagaga                                 >  1:519027/1‑100 (MQ=255)
                                                                 ttCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaagagaaa                               >  2:2227097/1‑100 (MQ=255)
                                                                   cGCTTATTGAACGCCGTGTACTAATGGCGCCTGGCACCTGCGATCGTCGCGGCCTGAATCCCTGGCTGGCGCTTCTGCCTATGGAAGCgaaagagaaaga                             <  2:2657262/100‑1 (MQ=255)
                                                                    gCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGAt                            >  2:1740709/1‑100 (MQ=255)
                                                                        aTTGAACGCCGTGTACGAAAGGCGCCTTGCCCAGTCGTTCGTCCCCATCTGAATACATAACAAACACGTCTGCATATAGAAACAAAAAAGAAAGATaaa                         >  1:1277525/1‑99 (MQ=255)
                                                                           gAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGc                     >  2:1258716/1‑100 (MQ=255)
                                                                            aaCGCCGTGTACGTTTTGCGCCTGTCACCTGCGATCGTCGCGGTCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCg                    <  1:486074/100‑1 (MQ=255)
                                                                             aCGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCg                    <  2:669302/99‑1 (MQ=255)
                                                                              cGCCGTGTACGAATGGCGCCTGGCACCGGCGATTGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGtt                  <  1:346917/100‑1 (MQ=255)
                                                                               gCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTa                 >  1:834329/1‑100 (MQ=255)
                                                                                ccGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTa                 <  1:1423762/99‑1 (MQ=255)
                                                                                ccGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTAc                >  2:1860636/1‑100 (MQ=255)
                                                                                ccGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTAc                >  2:1860592/1‑100 (MQ=255)
                                                                                ccGTGTACGAATGGCGCCGGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTAc                <  2:6368/100‑1 (MQ=255)
                                                                                  gtgtACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTg              <  1:1467419/100‑1 (MQ=255)
                                                                                   tgtACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTg              <  2:2839459/99‑1 (MQ=255)
                                                                                    gtACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCTTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGc            <  1:1882510/100‑1 (MQ=255)
                                                                                    gtACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGc            <  1:601868/100‑1 (MQ=255)
                                                                                    gtACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGc            <  2:382676/100‑1 (MQ=255)
                                                                                    ggaCGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGc            <  2:1903526/98‑1 (MQ=255)
                                                                                      aCGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAg          >  2:1711358/1‑100 (MQ=255)
                                                                                      aCGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAg          <  2:213925/100‑1 (MQ=255)
                                                                                        gAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAACGAAAGAGAAAGATGAATGACGTTACTGGCAGGc        >  1:1118982/1‑100 (MQ=255)
                                                                                         aaTGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCg       <  1:387489/100‑1 (MQ=255)
                                                                                         aaTGGCGCCTGGCACCGGCAATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCa           >  2:1841687/1‑96 (MQ=255)
                                                                                          aTGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCgg      <  1:1651469/100‑1 (MQ=255)
                                                                                             gcgcCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTATTGGCAGGCGGaat   <  1:2058882/100‑3 (MQ=255)
                                                                                             gcgcCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGAtt   <  2:140171/100‑1 (MQ=255)
                                                                                              cgcCTGGCACCGGCGTTCGTCGGGGCCTGAATACCTGGCTGGCTCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGAttt  <  2:1213097/100‑1 (MQ=255)
                                                                                                 |                                                                                                
GGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGATTT  >  NC_000913/4631603‑4631796

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: