Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,969,429 A→C 6.3% K133N (AAA→AAC wzzE → enterobacterial common antigen polysaccharide co‑polymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,969,4290AC 6.3% 889.7 / 39.7 286K133N (AAA→AACwzzEenterobacterial common antigen polysaccharide co‑polymerase
Reads supporting (aligned to +/- strand):  ref base A (145/123);  new base C (8/10);  total (153/133)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.06e-01

Only 100 of 296 total aligned reads displayed.

GTCCGTTCAAACATGGCTTCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTGCCAG  >  NC_000913/3969280‑3969578
                                                                                                                                                     |                                                                                                                                                     
 tCCGTTCAAACATGGCTTCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAg                                                                                                                                                      <  2:196552/150‑1 (MQ=255)
  ccGTTCAAACATGGCTTCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGc                                                                                                                                                     >  1:1368337/1‑150 (MQ=255)
       cAAACATGGCTTCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATg                                                                                                                                                >  1:1707082/1‑150 (MQ=255)
         aaCATGGCTTCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATgcg                                                                                                                                              >  1:905182/1‑150 (MQ=255)
          aCATGGCTTCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATgcgg                                                                                                                                             <  1:1389909/150‑1 (MQ=255)
                 ttCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAACGCCGCATCGGCTTTGCt                                                                                                                                      <  2:1546236/150‑1 (MQ=255)
                   cTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAACGCCGCATCGGCTTTGCTgt                                                                                                                                    <  1:1460278/150‑2 (MQ=255)
                    tGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGa                                                                                                                                   >  1:1647051/1‑150 (MQ=255)
                     gCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAACGCCGCATCGGCTTTGCTgttg                                                                                                                                  <  1:46642/150‑4 (MQ=255)
                     gCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAACGCCGCATCGGCTTTGCTgttg                                                                                                                                  <  1:664496/150‑4 (MQ=255)
                      ccGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTTTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCt                                                                                                                                      <  2:1574551/145‑1 (MQ=255)
                      ccGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTTGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATg                                                                                                                                 <  2:1415234/150‑1 (MQ=255)
                          ccaaccaTCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAAt                                                                                                                             >  1:643695/1‑150 (MQ=255)
                           caaccaTCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTTCAATCCGACTATTACAAACAGCGTATGTTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATg                                                                                                                            <  2:1041206/150‑1 (MQ=255)
                              ccaTCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGAtt                                                                                                                         <  1:858736/150‑1 (MQ=255)
                               caTCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTa                                                                                                                        <  1:943075/150‑1 (MQ=255)
                                  cGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTaaca                                                                                                                     <  1:1805848/150‑1 (MQ=255)
                                   ggTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTaacaa                                                                                                                    <  2:1444645/150‑1 (MQ=255)
                                   ggTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTaacaa                                                                                                                    <  2:1451627/150‑1 (MQ=255)
                                   ggTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCCACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTaacaa                                                                                                                    <  2:1269443/150‑1 (MQ=255)
                                        tGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATcc                                                                                                               <  1:1679789/150‑1 (MQ=255)
                                        tGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATcc                                                                                                               <  1:1409944/150‑1 (MQ=255)
                                        tGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATcc                                                                                                               <  1:18432/150‑1 (MQ=255)
                                           aCGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGt                                                                                                            <  1:1488663/150‑1 (MQ=255)
                                            cGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGtt                                                                                                           >  1:1925053/1‑150 (MQ=255)
                                            cGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGtt                                                                                                           <  1:265736/150‑1 (MQ=255)
                                             gAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGttt                                                                                                          <  2:1925053/150‑1 (MQ=255)
                                              aaGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTa                                                                                                         <  2:730818/150‑1 (MQ=255)
                                               aGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGt                                                                                                                                          >  1:238879/1‑116 (MQ=255)
                                                gCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATc                                                                                                       <  1:1578580/150‑1 (MQ=255)
                                                  cTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATccc                                                                                                     >  2:1660794/1‑150 (MQ=255)
                                                  cTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATccc                                                                                                     >  1:1332665/1‑150 (MQ=255)
                                                  cTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATccc                                                                                                     >  1:1812455/1‑150 (MQ=255)
                                                    ataAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCg                                                                                                   >  1:1111077/1‑150 (MQ=255)
                                                     taAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCgg                                                                                                  >  2:1172891/1‑150 (MQ=255)
                                                      aaaGAGTTTGTTATGCAACTGGCCTCGTGCGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGTGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGga                                                                                                 <  1:1888578/150‑1 (MQ=255)
                                                       aaGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGgag                                                                                                >  2:1233272/1‑150 (MQ=255)
                                                        agagTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGgaga                                                                                               >  1:30615/1‑150 (MQ=255)
                                                             ttGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGAAACAGCAAAGACGATGCGGCGTAGCTGGATGAAATGATTAACAAAAACAAGTGTATCACGCGAGACTGTa                                                                                          >  2:873527/1‑150 (MQ=255)
                                                              tGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTAc                                                                                         >  2:245104/1‑150 (MQ=255)
                                                               gTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCAGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTAcc                                                                                        >  2:1506554/1‑150 (MQ=255)
                                                                        cTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTc                                                                               <  1:125629/150‑1 (MQ=255)
                                                                          ggCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCaa                                                                             >  1:803714/1‑150 (MQ=255)
                                                                          ggCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCaa                                                                             <  2:30615/150‑1 (MQ=255)
                                                                          ggCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCaa                                                                             >  2:652884/1‑150 (MQ=255)
                                                                           gCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAAt                                                                            <  1:1506554/150‑1 (MQ=255)
                                                                            ccTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAACGCCGCATCGGCTTTGCTgt                                                                                                                                    >  2:983644/1‑92 (MQ=38)
                                                                            ccTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAACGCCGCATCGGCTTTGCTgt                                                                                                                                    <  1:983644/93‑2 (MQ=38)
                                                                             cTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGa                                                                          >  2:267140/1‑150 (MQ=255)
                                                                                 tGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGc                                                                      <  1:585125/150‑1 (MQ=255)
                                                                                        cccGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGc                                                               >  2:1343674/1‑150 (MQ=255)
                                                                                         ccGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCt                                                              >  1:1047452/1‑150 (MQ=255)
                                                                                         ccGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCt                                                              <  1:1750928/150‑1 (MQ=255)
                                                                                              gagagTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTg                                                         <  2:1380078/150‑1 (MQ=255)
                                                                                               agagTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGc                                                        >  2:512441/1‑150 (MQ=255)
                                                                                               agagTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGc                                                        <  1:547602/150‑1 (MQ=255)
                                                                                                 agTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCg                                                      <  2:1812455/150‑1 (MQ=255)
                                                                                                  gTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCcgcgcg                                                                                  <  2:1001894/121‑1 (MQ=255)
                                                                                                  gTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCcgcgcg                                                                                  >  1:1001894/1‑121 (MQ=255)
                                                                                                    tCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGaaa                                                   <  2:521308/150‑1 (MQ=255)
                                                                                                     cTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAAc                                                  >  2:60117/1‑150 (MQ=255)
                                                                                                      tGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAAcc                                                 <  2:1077290/150‑1 (MQ=255)
                                                                                                      tGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAAcc                                                 >  2:610055/1‑150 (MQ=255)
                                                                                                         cTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACcgcg                                              >  1:194431/1‑150 (MQ=255)
                                                                                                          tGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCgcgc                                             >  1:568142/1‑150 (MQ=255)
                                                                                                          tGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCgcgc                                             >  2:190550/1‑150 (MQ=255)
                                                                                                           gCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGcc                                            <  2:1643464/150‑1 (MQ=255)
                                                                                                            cAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGcc                                                       <  2:119265/138‑1 (MQ=255)
                                                                                                            cAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGcc                                                       >  1:119265/1‑138 (MQ=255)
                                                                                                            cAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCt                                           >  2:85521/1‑150 (MQ=255)
                                                                                                            cAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCt                                           >  2:670319/1‑150 (MQ=255)
                                                                                                              aaCCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGa                                         <  2:1687090/150‑1 (MQ=255)
                                                                                                                ccGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATcccc                                                                                                    <  2:473628/89‑1 (MQ=255)
                                                                                                                ccGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACg                                       >  2:1294547/1‑150 (MQ=255)
                                                                                                                      aTTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCtaata                                 >  2:15866/1‑150 (MQ=255)
                                                                                                                      aTTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCtaata                                 >  2:287372/1‑150 (MQ=255)
                                                                                                                         aCAAACTGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAAcc                              <  2:1398342/150‑1 (MQ=255)
                                                                                                                         aCAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAAcc                              >  2:1724909/1‑150 (MQ=255)
                                                                                                                          cAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCgcgc                                             <  2:918908/134‑1 (MQ=255)
                                                                                                                          cAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCt                             >  2:1363723/1‑150 (MQ=255)
                                                                                                                           aaaCAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTg                            <  1:1913916/150‑1 (MQ=255)
                                                                                                                             aCAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGtt                          >  2:238218/1‑150 (MQ=255)
                                                                                                                              cAGCTGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTa                         <  1:425985/150‑1 (MQ=255)
                                                                                                                                 cGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTg                                                                   >  1:850869/1‑105 (MQ=255)
                                                                                                                                     tggtggGCAACAGCAAAGCCGATGCGGCTTTGCTGGATGAAATGATTAACCACATCCATGTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAACCCGCGCCTGACGCTAATAACCTGTTACGTCAGt                  <  2:1767303/150‑1 (MQ=255)
                                                                                                                                      ggtggGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTa                 >  1:1075730/1‑150 (MQ=255)
                                                                                                                                       gtggGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTAt                >  2:761753/1‑150 (MQ=255)
                                                                                                                                        tggGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATg               >  2:456922/1‑150 (MQ=255)
                                                                                                                                        tggGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATg               >  2:701434/1‑150 (MQ=255)
                                                                                                                                           gCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTg            >  2:1939930/1‑150 (MQ=255)
                                                                                                                                           gCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTg            >  2:1839258/1‑150 (MQ=255)
                                                                                                                                           catccAGCAACGCCGCATCGGCTTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTg            >  2:3757/5‑150 (MQ=255)
                                                                                                                                            cAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGc           >  1:1568482/1‑150 (MQ=255)
                                                                                                                                               cAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGGCTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCttt        <  1:116880/150‑1 (MQ=255)
                                                                                                                                                aGCAACGCCGCATCGGCTTTGCTGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTg      <  1:225938/150‑1 (MQ=255)
                                                                                                                                                 gCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTg      >  2:113289/1‑150 (MQ=255)
                                                                                                                                                  cAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTgc     >  1:462193/1‑150 (MQ=255)
                                                                                                                                                   aaaGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTgct    <  2:8198/150‑2 (MQ=255)
                                                                                                                                                   aaaGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTgcc    <  1:73842/150‑1 (MQ=255)
                                                                                                                                                   aaaGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTgcc    >  2:341665/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTCCGTTCAAACATGGCTTCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTGCCAG  >  NC_000913/3969280‑3969578

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: