Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,521,374 G→C 7.1% 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 7.1% 264.1 / 8.1 80A216G (GCG→GGG) yjhFKpLE2 phage‑like element; putative transporter YjhF
Reads supporting (aligned to +/- strand):  ref base G (30/44);  new base C (2/4);  total (32/48)
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

CACGTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACATCAGCGCCATA  >  NC_000913/4521226‑4521522
                                                                                                                                                    |                                                                                                                                                    
cACGTTGTAAAATAAGCCTGAATCCGCCGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAAATTATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGCTcgc                                                                                                                                                     >  2:1285043/1‑150 (MQ=255)
 aCGTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTcgcg                                                                                                                                                    <  1:1110091/150‑1 (MQ=255)
   gTTGTAAAATAAGCCTGACTCCGCGGCACTTTTGGCATGAGTCATCTGGAGGATGGTGACAACCACCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGCGACCCGcc                                                                                                                                                  <  1:1959478/150‑1 (MQ=255)
   gTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGcc                                                                                                                                                  >  2:183449/1‑150 (MQ=255)
    ttGTAAAATAAGCCTGAATCGGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTTAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCg                                                                                                                                                 <  1:835481/150‑1 (MQ=255)
     tGTAAAATAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCga                                                                                                                                                >  2:1984675/1‑150 (MQ=255)
        aaaaTAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCgagag                                                                                                                                             >  2:323987/1‑150 (MQ=255)
         aaaTAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCgagaga                                                                                                                                            >  2:797189/1‑150 (MQ=255)
         aaaTAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGCGACCCGCCgattca                                                                                                                                            <  2:1252661/150‑5 (MQ=255)
          aaTAAGCCTGAATCCGCGTCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGCGACCCGCCgattcag                                                                                                                                           <  1:878482/150‑6 (MQ=255)
           aTAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAgg                                                                                                                                          <  1:381278/150‑1 (MQ=255)
           aTAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGCGACCCGCCgattcagc                                                                                                                                          <  2:762203/150‑7 (MQ=255)
              aGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCta                                                                                                                                       <  1:1903160/150‑1 (MQ=255)
                     aTCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTccc                                                                                                                                >  2:1516331/1‑150 (MQ=255)
                        cgcgGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCctc                                                                                                                             >  1:181508/1‑150 (MQ=255)
                        cgcgGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCctc                                                                                                                             >  1:739536/1‑150 (MQ=255)
                         gcgGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCctct                                                                                                                            <  1:1869510/150‑1 (MQ=255)
                             cACTTTTGGCATGTGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGAGGGTCGAGCCGAGAGAGGCTATATTCCCCTCttttt                                                                                                                        <  2:1048060/150‑1 (MQ=255)
                             cACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCttttt                                                                                                                        <  1:797189/150‑1 (MQ=255)
                               cTTTTGGCATGAGTCATCTGGATGATGGTGACAATCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTa                                                                                                                      <  2:973580/150‑1 (MQ=255)
                                 tttGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTa                                                                                                                      <  2:1432842/148‑1 (MQ=255)
                                 tttGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTa                                                                                                                      >  1:1432842/1‑148 (MQ=255)
                                  ttGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAg                                                                                                                   <  1:1318321/150‑1 (MQ=255)
                                  ttGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAg                                                                                                                   <  2:904008/150‑1 (MQ=255)
                                     gCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCg                                                                                                                <  2:987601/150‑1 (MQ=255)
                                     gCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCg                                                                                                                <  1:1803500/150‑1 (MQ=255)
                                      cATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGt                                                                                                               >  1:1459036/1‑150 (MQ=255)
                                      cATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGt                                                                                                               >  2:1580976/1‑150 (MQ=255)
                                        tGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTcc                                                                                                             <  1:1107583/150‑1 (MQ=255)
                                             cATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGaat                                                                                                        <  2:1404801/150‑1 (MQ=255)
                                             cATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGaat                                                                                                        >  1:91033/1‑150 (MQ=255)
                                             cATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGaat                                                                                                        <  1:827922/150‑1 (MQ=255)
                                                 tGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATg                                                                                                    >  1:37084/1‑150 (MQ=255)
                                                       atgGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTgg                                                                                              >  1:216773/1‑150 (MQ=255)
                                                        tgGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGc                                                                                             <  2:1683420/150‑1 (MQ=255)
                                                        tgGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGc                                                                                             <  2:181508/150‑1 (MQ=255)
                                                         ggTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCa                                                                                            >  1:1257173/1‑150 (MQ=255)
                                                             aCAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAAt                                                                                        <  1:1322058/150‑1 (MQ=255)
                                                              cAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAAtt                                                                                       >  2:1287467/1‑150 (MQ=255)
                                                              cAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAAtt                                                                                       >  2:259292/1‑150 (MQ=255)
                                                                 ccGCCCTCAGCATCACCGGCAGAATTGCAGTGAGAAAACTGATCCCAAAGCATGGCATTTAGCTTACGCTGAATCGGCGGGTCCCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGAGTCCGCAAAAAAGCGCTGGCAGAATTTaa                                                                                    >  2:113774/1‑149 (MQ=255)
                                                                  cGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTAtt                                                                                   <  1:1516331/150‑1 (MQ=255)
                                                                   gcaaTCAGGATCCCCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGGGGGTCGCGCCGAGAGAGGCTATATTCGCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTc                                                                                  <  1:1634114/147‑1 (MQ=255)
                                                                   gCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTc                                                                                  <  2:218767/150‑1 (MQ=255)
                                                                   gCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTc                                                                                  <  1:977732/150‑1 (MQ=255)
                                                                   gCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGGCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTc                                                                                  <  1:431638/150‑1 (MQ=255)
                                                                     cATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCaa                                                                                >  2:1284703/1‑150 (MQ=255)
                                                                            aTCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCggg                                                                         >  2:1887841/1‑150 (MQ=255)
                                                                                 cGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCa                                                                    <  1:511773/150‑1 (MQ=255)
                                                                                    caatgctGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAg                                                                 <  1:961914/144‑1 (MQ=255)
                                                                                       caTTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAAt                                                              >  2:1155935/1‑150 (MQ=255)
                                                                                            ccGTGAGAAATCTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAgg                                                         >  1:71965/1‑150 (MQ=255)
                                                                                            ccGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAgg                                                         <  2:739536/150‑1 (MQ=255)
                                                                                              gTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGtt                                                       <  1:1564158/150‑1 (MQ=255)
                                                                                                gagaAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTgg                                                     <  2:216773/150‑1 (MQ=255)
                                                                                                 agaAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGa                                                    >  2:962954/1‑150 (MQ=255)
                                                                                                  gaAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGaa                                                   <  1:1400546/150‑1 (MQ=255)
                                                                                                   aaaaCTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTGTATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAAt                                                  >  2:1276240/1‑150 (MQ=255)
                                                                                                   aaaaCTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAAt                                                  <  2:385910/150‑1 (MQ=255)
                                                                                                   aaaaCTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAAt                                                  >  1:740437/1‑150 (MQ=255)
                                                                                                     aaCTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCATAATTTATTCAATACCGTGCCGCAAAGAATATAGGTTGGAATGc                                                <  2:37084/150‑1 (MQ=255)
                                                                                                       cTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCg                                              >  1:463329/1‑150 (MQ=255)
                                                                                                       cTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCg                                              >  2:1315627/1‑150 (MQ=255)
                                                                                                            cccTAAGCGTGGCATTTGGGTTTCGCTGATTGGGCGGGTCGGGGAGAGTGAGGCTATATTCACGTCTTTTTTAAAGGCGTGCGGAATAATGCGCTGGCTGAATTTATTCAATACCGGGAGGCAAAGATTAAAGGTTGGAATGCCGACAAt                                         <  2:1255511/150‑1 (MQ=255)
                                                                                                            cccAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAAt                                         <  2:740437/150‑1 (MQ=255)
                                                                                                            cccAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAAt                                         >  2:412805/1‑150 (MQ=255)
                                                                                                                aaGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATa                                     >  1:1440174/1‑150 (MQ=255)
                                                                                                                   ccTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCg                                  <  2:1857533/150‑1 (MQ=255)
                                                                                                                   ccTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCg                                  <  2:1440174/150‑1 (MQ=255)
                                                                                                                    cTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGt                                 >  1:419394/1‑150 (MQ=255)
                                                                                                                    cTGGCATTTCGCTTTCGCTGAATCGGAGGGTCGCGCCGAGAGAGGCTATATTCCCATCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGt                                 >  1:115578/1‑150 (MQ=255)
                                                                                                                     tGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGta                                <  1:147217/150‑1 (MQ=255)
                                                                                                                       gCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGtata                              >  2:862952/1‑150 (MQ=255)
                                                                                                                         aTTTAGCTTTAGCTGAATCTGAGGTTCGCGACGAGAGATGCTATATTCCCATCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTAAATACCGGGAAGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATc                            <  2:1257173/150‑1 (MQ=255)
                                                                                                                           ttcgctttcgctGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAg                          <  2:1911868/150‑1 (MQ=255)
                                                                                                                               ctttcgctGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTAcc                      <  1:1984675/150‑1 (MQ=255)
                                                                                                                               ctttcgctGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTAcc                      <  1:252343/150‑1 (MQ=255)
                                                                                                                                  tcgctGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAAt                   <  1:1580976/150‑1 (MQ=255)
                                                                                                                                        tctCGGCGCGACCCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACa             >  1:1306334/3‑150 (MQ=255)
                                                                                                                                        aaTCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAATAATAAAGGTTGTAATGCCGACAATGATAGCGTATATCAGTAGCAATCCTACa             <  2:91033/150‑1 (MQ=255)
                                                                                                                                        aaTCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACa             <  1:962954/150‑1 (MQ=255)
                                                                                                                                        aaTCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACa             <  2:71965/150‑1 (MQ=255)
                                                                                                                                              cggcccgcgCCTCGAGTGCCTATATACCCCTCTTATTTAAACGCTTCCGCAATAATTCGCTGTCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACATCAgcg       <  1:861080/145‑1 (MQ=255)
                                                                                                                                                   cgcgCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACATCAGCGCCata  >  1:1764540/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CACGTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACATCAGCGCCATA  >  NC_000913/4521226‑4521522

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