Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,521,374 G→C 8.4% A216G (GCG→GGG)  yjhF ← KpLE2 phage‑like element; putative transporter YjhF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,521,3740GC 8.4% 271.4 / 13.3 85A216G (GCG→GGG) yjhFKpLE2 phage‑like element; putative transporter YjhF
Reads supporting (aligned to +/- strand):  ref base G (38/39);  new base C (5/2);  total (43/42)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

TCACGTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACATCAGCGCCA  >  NC_000913/4521225‑4521520
                                                                                                                                                     |                                                                                                                                                  
tCACGTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTcg                                                                                                                                                    <  1:1429672/150‑1 (MQ=255)
     ttGTAAAATAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGccc                                                                                                                                               >  1:589571/1‑149 (MQ=255)
     ttGTAAAATAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCg                                                                                                                                               <  1:1250087/150‑1 (MQ=255)
       gTAAAATAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCgag                                                                                                                                             <  2:1391433/150‑1 (MQ=255)
       gTAAAATAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCgag                                                                                                                                             >  2:357542/1‑150 (MQ=255)
          aaaTAAGCCTGAATCCGCGGGACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGCGACCCGCCgattca                                                                                                                                          <  1:386442/150‑5 (MQ=255)
           aaTAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGCGACTCGCCgattcag                                                                                                                                         <  1:1080236/150‑6 (MQ=255)
            aTAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGCGACCCGCCgattcagc                                                                                                                                        <  1:1146579/150‑7 (MQ=255)
             tAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGc                                                                                                                                       <  2:402434/150‑1 (MQ=255)
             tAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGGGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGc                                                                                                                                       <  1:307106/150‑1 (MQ=255)
                gCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCtat                                                                                                                                    >  1:507623/1‑150 (MQ=255)
                 ccTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCtata                                                                                                                                   <  2:55431/150‑1 (MQ=255)
                             gCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGGATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCtttt                                                                                                                       >  2:685409/1‑150 (MQ=255)
                               aCTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCtttttt                                                                                                                     <  1:1006255/150‑1 (MQ=255)
                               aCTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCtttttt                                                                                                                     <  1:1421093/150‑1 (MQ=255)
                                 ttttGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTaa                                                                                                                   <  2:1298283/150‑1 (MQ=255)
                                     ggCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGc                                                                                                               >  1:727274/1‑150 (MQ=255)
                                     ggCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCATTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAgcc                                                                                                               >  2:683673/1‑148 (MQ=255)
                                           aGTCATCTGGATGATGGTGACAACCGGCATCAGGATCACCGGCAGCATTTCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCTTCCgg                                                                                                         <  2:212637/150‑1 (MQ=255)
                                           aGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCgg                                                                                                         <  2:507623/150‑1 (MQ=255)
                                            gTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGa                                                                                                        >  2:1522209/1‑150 (MQ=255)
                                            gTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGa                                                                                                        <  1:357542/150‑1 (MQ=255)
                                              cATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGaat                                                                                                      >  2:1139709/1‑150 (MQ=255)
                                              cATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGaat                                                                                                      <  2:1025129/150‑1 (MQ=255)
                                                 cTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGaataat                                                                                                   >  2:208442/1‑150 (MQ=255)
                                                  tGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATACAAAAGACTGGCATTTCGCTTTCGCTGAATAGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATg                                                                                                  >  2:1590629/1‑150 (MQ=255)
                                                     atgatgGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATgcgc                                                                                               <  2:213474/150‑1 (MQ=255)
                                                        atgGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTgg                                                                                            >  2:982534/1‑150 (MQ=255)
                                                         tgGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGc                                                                                           <  1:426510/150‑1 (MQ=255)
                                                               cAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAAtt                                                                                     <  1:143848/150‑1 (MQ=255)
                                                               cAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAAtt                                                                                     <  2:1316002/150‑1 (MQ=255)
                                                                 aCCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTa                                                                                   >  1:1575141/1‑150 (MQ=255)
                                                                  ccGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTAt                                                                                  >  1:1212289/1‑150 (MQ=255)
                                                                  ccGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTAt                                                                                  <  2:1212289/150‑1 (MQ=255)
                                                                     ccATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCa                                                                               >  1:362123/1‑150 (MQ=255)
                                                                        tCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATa                                                                            >  2:301797/1‑150 (MQ=255)
                                                                        tCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATa                                                                            >  1:1255583/1‑150 (MQ=255)
                                                                        tCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATa                                                                            >  1:1151264/1‑150 (MQ=255)
                                                                         cAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATAc                                                                           >  2:1101111/1‑150 (MQ=255)
                                                                           ggATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCg                                                                         >  1:1557728/1‑150 (MQ=255)
                                                                           ggATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCg                                                                         >  2:1123180/1‑150 (MQ=255)
                                                                           ggATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCg                                                                         >  2:627082/1‑150 (MQ=255)
                                                                             aTCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCggg                                                                       >  2:488424/1‑150 (MQ=255)
                                                                              tCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGc                                                                      <  1:153177/150‑1 (MQ=255)
                                                                              tCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGc                                                                      <  1:964289/150‑1 (MQ=255)
                                                                                aCCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCg                                                                    >  2:1072908/1‑150 (MQ=255)
                                                                                 ccGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGCCCGc                                                                   <  1:1018675/150‑1 (MQ=255)
                                                                                   ggCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCaa                                                                 >  2:812062/1‑150 (MQ=255)
                                                                                    gcagcaTTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCaaa                                                                <  1:405449/150‑1 (MQ=255)
                                                                                     cagcaTTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGTCAGAATTTATTCAATACCGGGCCGCAAAg                                                               <  2:603729/150‑1 (MQ=255)
                                                                                     cagcaTTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAg                                                               <  1:1539507/150‑1 (MQ=255)
                                                                                         aTTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATa                                                           >  2:1319304/1‑150 (MQ=255)
                                                                                          ttGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGGATTTATTCAATACCGGGCCGCAAAGAATaa                                                          >  1:869392/1‑150 (MQ=255)
                                                                                            gCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAg                                                        <  2:340187/150‑1 (MQ=255)
                                                                                            gCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAg                                                        <  2:1079371/150‑1 (MQ=255)
                                                                                              cGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGt                                                      >  1:1126980/1‑150 (MQ=255)
                                                                                              cGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGt                                                      >  1:1120255/1‑150 (MQ=255)
                                                                                              cGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGt                                                      >  1:196943/1‑150 (MQ=255)
                                                                                                 gagaAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTgg                                                   <  2:196943/150‑1 (MQ=255)
                                                                                                     aaaCTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATg                                               >  1:248305/1‑150 (MQ=255)
                                                                                                       aCTGATCCCAAAGCCAGGCATTTCGCTTTCGCTCAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGcc                                             >  2:1014608/1‑150 (MQ=255)
                                                                                                        cTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGGATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCg                                            <  1:685409/150‑1 (MQ=255)
                                                                                                        cTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCg                                            <  2:352309/150‑1 (MQ=255)
                                                                                                        cTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCg                                            <  1:982534/150‑1 (MQ=255)
                                                                                                         tGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGa                                           <  1:627082/150‑1 (MQ=255)
                                                                                                          gATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGAc                                          >  2:676372/1‑150 (MQ=255)
                                                                                                          gATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGAc                                          >  2:1276546/1‑150 (MQ=255)
                                                                                                             cccAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAAt                                       >  1:315514/1‑150 (MQ=255)
                                                                                                                          aTTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATc                          <  1:1479358/150‑1 (MQ=255)
                                                                                                                          aTTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATc                          <  1:797723/150‑1 (MQ=255)
                                                                                                                           tttcgctttcgcTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCa                         <  1:208442/150‑1 (MQ=255)
                                                                                                                           tttcgctttcgcTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCa                         <  1:545326/150‑1 (MQ=255)
                                                                                                                              cgctttcgctGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCCTATATCAGTa                      <  1:1024192/150‑1 (MQ=255)
                                                                                                                               gctttcgctGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTAc                     >  1:1641729/1‑150 (MQ=255)
                                                                                                                                      ctctctCGGCGCGACCCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCt              >  1:775595/6‑150 (MQ=255)
                                                                                                                                       tctctCGGCGCGACCCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTa             >  1:1372759/5‑150 (MQ=255)
                                                                                                                                        gAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTAc            >  1:811256/1‑150 (MQ=255)
                                                                                                                                        ctctCGGCGCGACCCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTAc            >  1:30372/4‑150 (MQ=255)
                                                                                                                                         aaTCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACa           <  2:362123/150‑1 (MQ=255)
                                                                                                                                         aaTCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACa           <  2:1255583/150‑1 (MQ=255)
                                                                                                                                         aaTCGGCGGGGCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACa           <  1:1072908/150‑1 (MQ=255)
                                                                                                                                          ctCGGCGCGACCCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACAt          >  1:1320994/2‑150 (MQ=255)
                                                                                                                                          ctCGGCGCGACCCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACAt          >  1:1146381/2‑150 (MQ=255)
                                                                                                                                             ggcggGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACATCAg       >  1:656294/1‑150 (MQ=255)
                                                                                                                                               cggGTCGCGCCGAGAGAGGCTATATTCCCCTGTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGGAAAGAATAAAGGTTGGAATGCGGACAATGATAACGTATTTCAGTACCAATCCTACATCAgcg     <  2:523607/150‑1 (MQ=255)
                                                                                                                                                  gTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACATCAGCGCCa  >  1:166988/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TCACGTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACATCAGCGCCA  >  NC_000913/4521225‑4521520

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: ‑