Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,969,441 G→T 5.6% A137A (GCG→GCT 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,4410GT 5.6% 323.3 / 5.6 96A137A (GCG→GCTwzzEenterobacterial common antigen polysaccharide co‑polymerase
Reads supporting (aligned to +/- strand):  ref base G (44/45);  new base T (2/4);  total (47/49)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 102 total aligned reads displayed.

ATGGCTTCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTGCCAGCCAGCGTGCAGC  >  NC_000913/3969292‑3969590
                                                                                                                                                     |                                                                                                                                                     
aTGGCTTCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATgcggcg                                                                                                                                                       <  2:1473274/150‑1 (MQ=255)
    cTTCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGc                                                                                                                                                   <  1:569822/150‑1 (MQ=255)
      tCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAACGCCGCATCGGCTTTGCTg                                                                                                                                                 <  1:551590/150‑1 (MQ=255)
       cTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAACGCCGCATCGGCTTTGCTgt                                                                                                                                                <  2:1291204/150‑2 (MQ=255)
       cTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAACGCCGCATCGGCTTTGCTgt                                                                                                                                                <  1:947964/150‑2 (MQ=255)
       cTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTgg                                                                                                                                                <  1:164277/150‑1 (MQ=255)
       cTGCCGACCAACCATCGGTCATGGACGAAGCATATAAAGAGTTTGTTATGCAACTGGAATCGTGGGATACACGCAGAGAGTTCTCGCTGCAAACCGACTATTACAAACAGAGGATGGAGGCCAACAGCAAAGCCGATTCCGCGTTGCTgg                                                                                                                                                >  2:1322452/1‑150 (MQ=255)
         gCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGAt                                                                                                                                              <  1:628099/150‑1 (MQ=255)
         gCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGAt                                                                                                                                              >  2:1165512/1‑150 (MQ=255)
          ccGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGACCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATg                                                                                                                                             >  1:723330/1‑150 (MQ=255)
               caaccaTCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTgg                                                                                                                                                >  1:1207903/1‑142 (MQ=255)
               caaccaTCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTgg                                                                                                                                                <  2:1207903/142‑1 (MQ=255)
                aaccaTCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGAt                                                                                                                                      >  2:343546/1‑150 (MQ=255)
                 accaTCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGAt                                                                                                                                      <  1:464991/150‑1 (MQ=255)
                  ccaTCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGAtt                                                                                                                                     <  1:1286273/150‑1 (MQ=255)
                  ccaTCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGAtt                                                                                                                                     >  2:668550/1‑150 (MQ=255)
                  ccaTCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGAtt                                                                                                                                     <  1:1123228/150‑1 (MQ=255)
                   caTCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTa                                                                                                                                    >  1:1511711/1‑150 (MQ=255)
                     tCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTaac                                                                                                                                  <  2:35354/150‑1 (MQ=255)
                       ggTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTaacaa                                                                                                                                <  1:1441914/150‑1 (MQ=255)
                          cATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACAt                                                                                                                             >  2:1138715/1‑150 (MQ=255)
                          cATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACAt                                                                                                                             >  1:173945/1‑150 (MQ=255)
                           aTGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCGGATGCGGCGTTGCTGGATGAAATGATTAACAACATc                                                                                                                            <  1:1010455/150‑1 (MQ=255)
                           aTGGACCAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGCTACACGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACACCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATc                                                                                                                            >  2:271554/1‑150 (MQ=255)
                            tGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATcc                                                                                                                           <  1:804701/150‑1 (MQ=255)
                            tGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATcc                                                                                                                           >  1:1377617/1‑150 (MQ=255)
                             ggACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCa                                                                                                                          >  2:576803/1‑150 (MQ=255)
                             ggACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCa                                                                                                                          <  1:1237652/150‑1 (MQ=255)
                             ggACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCa                                                                                                                          >  2:1342158/1‑150 (MQ=255)
                              gACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAg                                                                                                                         >  2:1506394/1‑150 (MQ=255)
                                cGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGtt                                                                                                                       <  2:473106/150‑1 (MQ=255)
                                     ccTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATcc                                                                                                                  >  1:149368/1‑150 (MQ=255)
                                     ccTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATcc                                                                                                                  >  1:1429603/1‑150 (MQ=255)
                                     ccTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATcc                                                                                                                  <  1:1506394/150‑1 (MQ=255)
                                       tataAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATcccc                                                                                                                >  1:383948/1‑150 (MQ=255)
                                       tataAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATcccc                                                                                                                >  1:883028/1‑150 (MQ=255)
                                       tataAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATcccc                                                                                                                >  1:229314/1‑150 (MQ=255)
                                        ataAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATcc                                                                                                                           >  2:1204330/1‑138 (MQ=255)
                                        ataAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATcc                                                                                                                           <  1:1204330/138‑1 (MQ=255)
                                          aaaGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTa                                                                                                                     >  2:1042306/1‑142 (MQ=255)
                                          aaaGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTa                                                                                                                     <  1:1042306/142‑1 (MQ=255)
                                            agagTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGgaga                                                                                                           >  1:277438/1‑150 (MQ=255)
                                                tttGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCTGAGACttt                                                                                                       <  2:661352/150‑1 (MQ=255)
                                                    ttATGGAACTGGCGTCGTGGGATACCCGCAGATAGTTCTGTCTGCAAACCGAATATTACAAAGTGCGGATGGTGTTCAACAGCATAGCCGATGCGGCGTTTCTGGATGAAATGATTAACAACTTCCAGTTTATCCCAGGTTACTTTAcct                                                                                                   <  2:147859/150‑2 (MQ=255)
                                                     tATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCcg                                                                                                  <  1:1346314/150‑1 (MQ=255)
                                                     tATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCcg                                                                                                  <  1:306216/150‑1 (MQ=255)
                                                      aTGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCcgc                                                                                                 >  2:29765/1‑150 (MQ=255)
                                                        gCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCGGTTTATCCCCGGAGACTTTACCcgcgc                                                                                               <  2:733094/150‑1 (MQ=255)
                                                         cAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCcgcgcg                                                                                              <  1:504611/150‑1 (MQ=255)
                                                         cAACTGACCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCcgcgcg                                                                                              <  2:723330/150‑1 (MQ=255)
                                                               gCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCCTTGCTGGATGAAATGATTAAAAACATACACTTAATCACAGGAGACATTAAACGCACCGTCAAt                                                                                        >  1:656641/1‑150 (MQ=255)
                                                                ccTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATg                                                                                       >  2:433422/1‑150 (MQ=255)
                                                                ccTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATg                                                                                       <  2:1318856/150‑1 (MQ=255)
                                                                    gTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAg                                                                                   >  2:483813/1‑150 (MQ=255)
                                                                      gggATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCg                                                                                 <  1:1165512/150‑1 (MQ=255)
                                                                       ggATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGt                                                                                >  1:1532226/1‑150 (MQ=255)
                                                                          tACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCTGTCATTGACcgcgcggg                                                                             <  1:817637/150‑5 (MQ=255)
                                                                                  gagagTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTg                                                                     >  1:1231825/1‑150 (MQ=255)
                                                                                  gagagTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTg                                                                     >  1:1307150/1‑150 (MQ=255)
                                                                                    tagTTCTCGCTGCACACCGAATATATAAAATATCGTTGGGTTGTCAACAGCAAAGCGGATGCGTTTTTGCGGTATTAAATGATTATCATCATCCAGTTTATCCCCGGAGACTTTACCCGCTCGGTCAATGAGATCGTGAATGTTATTTcc                                                                   <  2:656641/149‑1 (MQ=255)
                                                                                    gagTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGcc                                                                   <  1:69520/150‑1 (MQ=255)
                                                                                     agTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCg                                                                  <  2:149368/150‑1 (MQ=255)
                                                                                      gTTCTGGCTGCAAACCGACTATTACAAACCGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGa                                                                 >  2:1216520/1‑150 (MQ=255)
                                                                                       ttCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGaa                                                                <  2:1300836/150‑1 (MQ=255)
                                                                                        tCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGaaa                                                               >  2:893652/1‑150 (MQ=255)
                                                                                        tCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGaaa                                                               >  2:111353/1‑150 (MQ=255)
                                                                                         cTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAAc                                                              >  2:733541/1‑150 (MQ=255)
                                                                                          tGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAAcc                                                             <  2:173945/150‑1 (MQ=255)
                                                                                             cTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACcgcg                                                          >  1:1262527/1‑150 (MQ=255)
                                                                                             cTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACcgcc                                                          >  2:80614/1‑149 (MQ=255)
                                                                                              tGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCgcgc                                                         <  2:1532226/150‑1 (MQ=255)
                                                                                               gCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGcc                                                        <  1:576803/150‑1 (MQ=255)
                                                                                                 aaaCCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTg                                                      >  1:506909/1‑150 (MQ=255)
                                                                                                      gACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCt                                                 >  2:1048234/1‑150 (MQ=255)
                                                                                                      gACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCt                                                 <  2:277438/150‑1 (MQ=255)
                                                                                                      gACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCt                                                 <  1:1259741/150‑1 (MQ=255)
                                                                                                        cTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCtaa                                               >  2:1385266/1‑150 (MQ=255)
                                                                                                          aTTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCtaata                                             <  2:229314/150‑1 (MQ=255)
                                                                                                              cAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCt                                         >  2:1388270/1‑150 (MQ=255)
                                                                                                               aaaCAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTg                                        <  2:883028/150‑1 (MQ=255)
                                                                                                               aaaCAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTg                                        <  1:1138715/150‑1 (MQ=255)
                                                                                                                aaCAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGt                                       <  1:360891/150‑1 (MQ=255)
                                                                                                                   aGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTAc                                    >  2:590399/1‑150 (MQ=255)
                                                                                                                       gATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCa                                >  1:847476/1‑150 (MQ=255)
                                                                                                                        aTGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAg                               <  2:44704/150‑1 (MQ=255)
                                                                                                                            tttcatccAGCAACGCCGCATCGGCTTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATg                           >  1:1217546/8‑150 (MQ=255)
                                                                                                                            tggGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATg                           >  2:913235/1‑150 (MQ=255)
                                                                                                                            tggGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATg                           <  2:1262527/150‑1 (MQ=255)
                                                                                                                              ggCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGtt                         <  2:202639/150‑1 (MQ=255)
                                                                                                                                cAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTgg                       >  1:286495/1‑149 (MQ=255)
                                                                                                                                cAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGc                       <  1:1216520/150‑1 (MQ=255)
                                                                                                                                  ccAGCAACGCCGCATCGGCTTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCtt                     >  2:1291782/2‑150 (MQ=255)
                                                                                                                                      cAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTgc                 <  1:483813/150‑1 (MQ=255)
                                                                                                                                       aaaGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTgcc                >  2:774688/1‑150 (MQ=255)
                                                                                                                                       aaaGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTgcc                <  2:791662/150‑1 (MQ=255)
                                                                                                                                                 cggcgTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTGCCAGCCAGCGTg      <  1:1342158/150‑1 (MQ=255)
                                                                                                                                                 cggcgTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTGCCAGCCAGCGTg      <  2:1231825/150‑1 (MQ=255)
                                                                                                                                                   gcgTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTGCCAGCCAGCGTGca    <  1:590399/150‑1 (MQ=255)
                                                                                                                                                    cgTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTGCCAGTCAGCGTGcag   <  2:1511711/150‑1 (MQ=255)
                                                                                                                                                     gTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTGCCAGCCAGCGTGcagc  <  1:668550/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATGGCTTCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTGCCAGCCAGCGTGCAGC  >  NC_000913/3969292‑3969590

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: