Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 572,454 C→T C29C (TGC→TGT ylcH → uncharacterized protein, DLP12 prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913572,4540CT98.8% 269.9 / NA 84C29C (TGC→TGTylcHuncharacterized protein, DLP12 prophage
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (47/36);  minor base G (1/0);  total (48/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTGCTTAAACCGTGGCGTGAGA  >  NC_000913/572320‑572589
                                                                                                                                      |                                                                                                                                        
ccGGAGCCTGAACAACTCGGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAg                                                                                                                                     >  1:1018620/1‑139 (MQ=255)
         gAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGt                                                                                                                            <  2:1555101/139‑1 (MQ=255)
              acTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAAcc                                                                                                                       >  2:1609715/1‑139 (MQ=255)
               cTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACct                                                                                                                      >  1:1254951/1‑139 (MQ=255)
               cTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACct                                                                                                                      >  2:675619/1‑139 (MQ=255)
                        cTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGa                                                                                                             <  2:1254951/139‑1 (MQ=255)
                              cgcTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCAt                                                                                                       >  2:1033411/1‑139 (MQ=255)
                                        ggACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGaa                                                                                                                         >  2:1781597/1‑111 (MQ=255)
                                        ggACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGaa                                                                                                                         <  1:1781597/111‑1 (MQ=255)
                                              ttATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGc                                                                                       <  2:518066/139‑1 (MQ=255)
                                                        aTACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGa                                                                                                             >  2:597155/1‑107 (MQ=255)
                                                        aTACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGa                                                                                                             <  1:597155/107‑1 (MQ=255)
                                                        aTACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCg                                                                             >  1:1504188/1‑139 (MQ=255)
                                                                  aaCTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTAtcggtcg                                                                   >  2:1779767/1‑139 (MQ=255)
                                                                  aaCTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTAtcggtcg                                                                   >  2:1749053/1‑139 (MQ=255)
                                                                  aaCTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTAtcggtcg                                                                   >  1:991501/1‑139 (MQ=255)
                                                                   aCTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTAtcggtcgg                                                                  <  2:991501/139‑1 (MQ=255)
                                                                   aCTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTAtcggtcgg                                                                  <  2:1504188/139‑1 (MQ=255)
                                                                      ctctTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAggg                                                                                                                                   <  1:1733541/71‑1 (MQ=255)
                                                                      ctctTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAggg                                                                                                                                   >  2:1733541/1‑71 (MQ=255)
                                                                        cttctCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCgg                                                                                        >  2:173298/1‑112 (MQ=255)
                                                                        cttctCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCgg                                                                                        <  1:173298/112‑1 (MQ=255)
                                                                        cttctCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGAt                                                             >  1:220446/1‑139 (MQ=255)
                                                                        cttctCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGAt                                                             <  2:1018620/139‑1 (MQ=255)
                                                                          tctcCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATcc                                                           >  2:21451/1‑139 (MQ=255)
                                                                            tcCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGc                                                                                       >  1:1801984/1‑109 (MQ=255)
                                                                            tcCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGc                                                                                       <  2:1801984/109‑1 (MQ=255)
                                                                             tccctCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGc                                                        >  2:802865/5‑139 (MQ=255)
                                                                                 tCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTa                                                    >  2:1642380/1‑139 (MQ=255)
                                                                                 tCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTa                                                    >  1:1616214/1‑139 (MQ=255)
                                                                                   acaGATGCAGAGATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTAt                                                                         >  1:1064881/1‑116 (MQ=255)
                                                                                   acaGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTAt                                                                         <  2:1064881/116‑1 (MQ=255)
                                                                                     aGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGcc                                                       >  1:797315/1‑132 (MQ=255)
                                                                                     aGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGcc                                                       <  2:797315/132‑1 (MQ=255)
                                                                                       aTGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCCACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTgga                                              >  2:1781809/1‑139 (MQ=255)
                                                                                       aTGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGCTCGGCAGATCCAGCCTTATCTgga                                              >  1:749333/1‑139 (MQ=255)
                                                                                        tGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTAtcggtcgg                                                                  <  2:1092069/118‑1 (MQ=255)
                                                                                        tGCAGAAATGCACCTGCAATTCTTTGCATCTACCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTAtcggtcgg                                                                  >  1:1092069/1‑118 (MQ=255)
                                                                                          cAGAAATGCACCTGCAATTCTTTTCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGcc                                                       <  2:1107468/127‑1 (MQ=255)
                                                                                          cAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCa                                                                                                        <  1:1278836/78‑1 (MQ=255)
                                                                                          cAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCa                                                                                                        >  2:1278836/1‑78 (MQ=255)
                                                                                          cAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTa                                                                          <  1:1636723/108‑1 (MQ=255)
                                                                                          cAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTa                                                                          >  2:1636723/1‑108 (MQ=255)
                                                                                          cAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGcc                                                       <  1:1358508/127‑1 (MQ=255)
                                                                                          cAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGcc                                                       >  1:1107468/1‑127 (MQ=255)
                                                                                          cAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGcc                                                       >  2:1358508/1‑127 (MQ=255)
                                                                                          cAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTggagga                                           >  2:101780/1‑139 (MQ=255)
                                                                                             aaaTGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGaggagg                                         >  2:30986/1‑139 (MQ=255)
                                                                                             aaaTGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCAc                                                                               >  2:213666/1‑100 (MQ=255)
                                                                                             aaaTGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCAc                                                                               <  1:213666/100‑1 (MQ=255)
                                                                                                 gCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCaaa                                    >  2:1823973/1‑139 (MQ=255)
                                                                                                  cACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGAt                                                             <  1:1693891/113‑1 (MQ=255)
                                                                                                  cACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGAt                                                             >  2:1693891/1‑113 (MQ=255)
                                                                                                    ccTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCa                                                                                                        >  2:1512388/1‑68 (MQ=255)
                                                                                                    ccTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCa                                                                                                        <  1:1512388/68‑1 (MQ=255)
                                                                                                     cTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCt                                >  1:700650/1‑139 (MQ=255)
                                                                                                       gCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCa                                                                                                        >  2:332277/1‑65 (MQ=255)
                                                                                                       gCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCa                                                                                                        <  1:332277/65‑1 (MQ=255)
                                                                                                       gCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCttt                              >  1:71029/1‑139 (MQ=255)
                                                                                                       gCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCttt                              >  2:1074233/1‑139 (MQ=255)
                                                                                                        cAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTc                             <  1:1609715/139‑1 (MQ=255)
                                                                                                           ttCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAAcc                                                                                                                       >  1:1083677/1‑46 (MQ=255)
                                                                                                           ttCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAAcc                                                                                                                       <  2:1083677/46‑1 (MQ=255)
                                                                                                            tCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGTAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCAt                         <  1:1749053/139‑1 (MQ=255)
                                                                                                            tCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCAt                         <  1:1074233/139‑1 (MQ=255)
                                                                                                            tCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCAt                         <  2:1396409/139‑1 (MQ=255)
                                                                                                            tCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCAt                         <  1:1033411/139‑1 (MQ=255)
                                                                                                            tCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTGATCTGGAGGAGGGCAAATGCTTTCGCAt                         <  2:749333/139‑1 (MQ=255)
                                                                                                             cTTTGCATCTAGCGTTTGGCCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATg                        <  2:1616214/139‑1 (MQ=255)
                                                                                                             cTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATg                        <  2:71029/139‑1 (MQ=255)
                                                                                                               ttGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGAt                                                             >  2:964145/1‑100 (MQ=255)
                                                                                                               ttGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGAt                                                             <  1:964145/100‑1 (MQ=255)
                                                                                                                     cTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTGCTTaa                >  1:473974/1‑139 (MQ=255)
                                                                                                                       aGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATc                                                            >  2:951134/1‑93 (MQ=255)
                                                                                                                       aGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATc                                                            <  1:951134/93‑1 (MQ=255)
                                                                                                                       aGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGcc                                                       <  1:511777/98‑1 (MQ=255)
                                                                                                                       aGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGcc                                                       >  2:511777/1‑98 (MQ=255)
                                                                                                                       aGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTggag                                             >  1:805002/1‑108 (MQ=255)
                                                                                                                       aGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTggag                                             <  2:805002/108‑1 (MQ=255)
                                                                                                                       aGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGTGCTTAAAcc             >  1:1717689/1‑139 (MQ=255)
                                                                                                                        gCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATc                                                            <  1:1280396/92‑1 (MQ=255)
                                                                                                                        gCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATc                                                            >  2:1280396/1‑92 (MQ=255)
                                                                                                                               aCCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTGCTTAAACCGTGGCGt      >  1:1695751/1‑139 (MQ=255)
                                                                                                                                   ctGGGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGAAGGAGGGCAAATGCTTTCGCATGGTGCTTAAACCGTGGCGTgaga  >  1:1179041/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                        
CCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTGCTTAAACCGTGGCGTGAGA  >  NC_000913/572320‑572589

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: