Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,643,737 A→T 5.7% L190I (TTA→ATA) ‡ ydfU ← Qin prophage; protein YdfU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,643,7370AT 5.7% 277.6 / 5.3 87L190I (TTA→ATA) ‡ydfUQin prophage; protein YdfU
Reads supporting (aligned to +/- strand):  ref base A (41/41);  new base T (3/2);  total (44/43)
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.14e-01

CGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTG  >  NC_000913/1643588‑1643880
                                                                                                                                                     |                                                                                                                                               
cGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTaa                                                                                                                                                 <  1:221621/150‑1 (MQ=255)
       tCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGc                                                                                                                                          <  2:1380091/150‑1 (MQ=255)
       tCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTGGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGc                                                                                                                                          <  2:1492436/150‑1 (MQ=255)
        ccTTTGGTTATTTTTGTTGCTGGTTTTCGTAATGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCg                                                                                                                                         <  2:1563415/150‑1 (MQ=255)
        ccTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCg                                                                                                                                         >  1:58124/1‑150 (MQ=255)
        ccTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCg                                                                                                                                         <  2:127047/150‑1 (MQ=255)
           gtggtgATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAac                                                                                                                                      <  1:3831/150‑1 (MQ=255)
           gtggtgATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGGATTACCAGCggaa                                                                                                                                      <  1:222715/150‑4 (MQ=255)
           gtggtgATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGGATTACCAGCggaa                                                                                                                                      <  2:655884/150‑4 (MQ=255)
              gtgATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAacgac                                                                                                                                   >  1:857869/1‑150 (MQ=255)
                gATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGAcgc                                                                                                                                 >  1:975159/1‑150 (MQ=255)
                      tgttgCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACa                                                                                                                           <  2:1655382/150‑1 (MQ=255)
                      tgttgCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACa                                                                                                                           <  1:988660/150‑1 (MQ=255)
                      tgttgCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACa                                                                                                                           <  2:58124/150‑1 (MQ=255)
                      tgttgCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACa                                                                                                                           <  1:208265/150‑1 (MQ=255)
                       gttgCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAg                                                                                                                          >  1:700102/1‑150 (MQ=255)
                          gCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCt                                                                                                                       >  1:1376682/1‑150 (MQ=255)
                           cTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTc                                                                                                                      <  2:304354/150‑1 (MQ=255)
                            tggtggGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCa                                                                                                                     <  2:297076/150‑1 (MQ=255)
                            tggtggGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCa                                                                                                                     <  1:1456303/150‑1 (MQ=255)
                              gtggGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAgg                                                                                                                   <  1:673359/150‑1 (MQ=255)
                               tggGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGt                                                                                                                  >  2:190036/1‑150 (MQ=255)
                               tggGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGt                                                                                                                  <  1:1124555/150‑1 (MQ=255)
                                gggCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTa                                                                                                                 >  1:1085939/1‑150 (MQ=255)
                                   cGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTg                                                                                                              >  1:235477/1‑150 (MQ=255)
                                    gTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGc                                                                                                             <  2:1134146/150‑1 (MQ=255)
                                     tAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCa                                                                                                            >  1:661131/1‑150 (MQ=255)
                                     tAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCa                                                                                                            <  2:241278/150‑1 (MQ=255)
                                       aGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATc                                                                                                          >  2:358918/1‑150 (MQ=255)
                                        ggCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAg                                                                                                                          >  1:176483/1‑133 (MQ=255)
                                        ggCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAg                                                                                                                          <  2:176483/133‑1 (MQ=255)
                                        ggCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATcc                                                                                                         <  2:1402547/150‑1 (MQ=255)
                                        ggCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATcc                                                                                                         >  1:26880/1‑150 (MQ=255)
                                        ggCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATcc                                                                                                         >  1:20713/1‑150 (MQ=255)
                                        ggCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATcc                                                                                                         >  1:988405/1‑150 (MQ=255)
                                        ggCAACCGAGCAACATTTTTTCTGCGTTGTATCAATATGCTGGAGGCGGTCAGCTCTCCCGGTACGATGTCGCTATCGCGGTACACCGAGCAGATATTATCCCCTGGTAATCCCAGCGAACGACGAGATACCTCCTCAGGTAGTGCaaaa                                                                                                         >  2:155167/1‑147 (MQ=255)
                                             cGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCcaccac                                                                                                    <  2:327389/150‑1 (MQ=255)
                                               gggCAAGATTTTTTGTGCGTTGGTTCAATATGCTGGTGGCGGTCAGCGCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCt                                                                                                  <  2:796547/148‑1 (MQ=255)
                                                tGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATcc                                                                                                         >  2:1577384/1‑142 (MQ=255)
                                                tGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATcc                                                                                                         <  1:1577384/142‑1 (MQ=255)
                                                  ccaGATTTTTTGTGCGTTGTTTCAATCTGCTGGTGGCGGTCAGCGCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGAt                                                                                               <  2:308574/148‑1 (MQ=255)
                                                  cAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGAt                                                                                               <  1:453555/150‑1 (MQ=255)
                                                    aGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTGGTGGATGCACTACCTGaggc                                                                                             <  2:1617665/150‑4 (MQ=255)
                                                     gATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGc                                                                                            <  1:1294507/150‑1 (MQ=255)
                                                     gATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGc                                                                                            >  2:1242938/1‑150 (MQ=255)
                                                     gATTTTTTGTGCGTTGTTTCAATAAGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGc                                                                                            >  1:705745/1‑150 (MQ=255)
                                                      aTTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCa                                                                                           >  2:1239864/1‑150 (MQ=255)
                                                              tgCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGAtt                                                                                              <  2:281432/139‑1 (MQ=255)
                                                              tgCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGAtt                                                                                              >  1:281432/1‑139 (MQ=255)
                                                                 gttgttTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCccacc                                                                                >  2:1266013/1‑150 (MQ=255)
                                                                    gttTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGCCGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAgg                                                                             >  1:1611761/1‑150 (MQ=255)
                                                                          atatGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAAtt                                                                       >  1:571669/1‑150 (MQ=255)
                                                                           tatGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTc                                                                      >  2:1369257/1‑150 (MQ=255)
                                                                            atGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCa                                                                     <  2:522748/150‑1 (MQ=255)
                                                                             tGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAg                                                                    >  1:722029/1‑150 (MQ=255)
                                                                             tGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAg                                                                    >  1:1323132/1‑150 (MQ=255)
                                                                             tGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAg                                                                    >  1:664243/1‑150 (MQ=255)
                                                                                tggtggCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCa                                                                 >  1:565119/1‑150 (MQ=255)
                                                                                          cAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTc                                                       >  1:1351900/1‑150 (MQ=255)
                                                                                            gCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTTGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGGATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTccc                                                     <  2:407292/150‑1 (MQ=255)
                                                                                            gCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTccc                                                     >  1:1228863/1‑150 (MQ=255)
                                                                                                 cccGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTcc                                                >  2:1367466/1‑150 (MQ=255)
                                                                                                     gTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCg                                            >  1:269931/1‑150 (MQ=255)
                                                                                                     gTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCg                                            <  2:661131/150‑1 (MQ=255)
                                                                                                     gTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCg                                            <  2:269931/150‑1 (MQ=255)
                                                                                                       aCGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTa                                          >  2:754175/1‑150 (MQ=255)
                                                                                                        cGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTAc                                         <  2:700102/150‑1 (MQ=255)
                                                                                                        cGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTAc                                         <  1:781566/150‑1 (MQ=255)
                                                                                                        cGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTAc                                         >  1:202932/1‑150 (MQ=255)
                                                                                                         gATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTAcc                                        >  2:1276593/1‑150 (MQ=255)
                                                                                                               gCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTAt                                  >  2:381783/1‑150 (MQ=255)
                                                                                                                   tCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCg                              >  2:142498/1‑150 (MQ=255)
                                                                                                                   tCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCg                              <  1:1239864/150‑1 (MQ=255)
                                                                                                                   tCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCg                              <  1:143654/150‑1 (MQ=255)
                                                                                                                    ctcgGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTTATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGt                             <  2:1127461/148‑1 (MQ=255)
                                                                                                                        gTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACg                         <  2:1569424/150‑1 (MQ=255)
                                                                                                                        gTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACg                         <  2:1376682/150‑1 (MQ=255)
                                                                                                                           cacCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGAt                      >  2:531564/1‑150 (MQ=255)
                                                                                                                            acCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATg                     <  1:1521849/150‑1 (MQ=255)
                                                                                                                               gAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTATTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACg                  <  1:567143/150‑1 (MQ=255)
                                                                                                                               gAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACg                  <  2:975159/150‑1 (MQ=255)
                                                                                                                                  cAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCa               >  2:197672/1‑150 (MQ=255)
                                                                                                                                      ttttttCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCa           >  1:950017/1‑150 (MQ=255)
                                                                                                                                          gttcGCTGGGATTACCAGCGGACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCATACCGCTTATTGCGTGACGGATGACGTCAATCACCCa       >  1:1062851/4‑150 (MQ=255)
                                                                                                                                          gttcGCTGGGATTACCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCa       >  2:1320609/4‑150 (MQ=255)
                                                                                                                                          gttcGCTGGGATTACCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCa       >  1:598157/4‑150 (MQ=255)
                                                                                                                                            ccGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATg     >  2:116694/1‑150 (MQ=255)
                                                                                                                                               cTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTg  <  1:1242938/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
CGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTG  >  NC_000913/1643588‑1643880

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