Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,969,436 G→T 5.8% A136S (GCG→TCG)  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,4360GT 5.8% 328.7 / 4.1 96A136S (GCG→TCG) wzzEenterobacterial common antigen polysaccharide co‑polymerase
Reads supporting (aligned to +/- strand):  ref base G (44/46);  new base T (3/3);  total (47/49)
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

AACATGGCTTCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTGCC  >  NC_000913/3969289‑3969576
                                                                                                                                                   |                                                                                                                                            
aagaTGGCTTCTGCCGACCAACCATCGGTGATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATgcg                                                                                                                                            <  1:1232387/147‑1 (MQ=255)
aaCATGGCTTCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATgcg                                                                                                                                            >  2:1123228/1‑150 (MQ=255)
   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:1123228/150‑1 (MQ=255)
                     ccaTCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGAtt                                                                                                                       <  1:1286273/150‑1 (MQ=255)
                     ccaTCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGAtt                                                                                                                       >  2:668550/1‑150 (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:1342158/1‑150 (MQ=255)
                                ggACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCa                                                                                                            >  2:576803/1‑150 (MQ=255)
                                ggACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCa                                                                                                            <  1:1237652/150‑1 (MQ=255)
                                 gACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAg                                                                                                           >  2:1506394/1‑150 (MQ=255)
                                   cGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGtt                                                                                                         <  2:473106/150‑1 (MQ=255)
                                        ccTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATcc                                                                                                    <  1:1506394/150‑1 (MQ=255)
                                        ccTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATcc                                                                                                    >  1:149368/1‑150 (MQ=255)
                                        ccTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATcc                                                                                                    >  1:1429603/1‑150 (MQ=255)
                                          tataAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATcccc                                                                                                  >  1:229314/1‑150 (MQ=255)
                                          tataAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATcccc                                                                                                  >  1:883028/1‑150 (MQ=255)
                                          tataAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATcccc                                                                                                  >  1:383948/1‑150 (MQ=255)
                                           ataAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATcc                                                                                                             <  1:1204330/138‑1 (MQ=255)
                                           ataAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATcc                                                                                                             >  2:1204330/1‑138 (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:1318856/150‑1 (MQ=255)
                                                                   ccTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATg                                                                         >  2:433422/1‑150 (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:111353/1‑150 (MQ=255)
                                                                                           tCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGaaa                                                 >  2:893652/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)
                                                                                                     aaCCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGa                                       >  2:993253/1‑150 (MQ=255)
                                                                                                         gACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCt                                   <  1:1259741/150‑1 (MQ=255)
                                                                                                         gACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCt                                   <  2:277438/150‑1 (MQ=255)
                                                                                                         gACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCt                                   >  2:1048234/1‑150 (MQ=255)
                                                                                                           cTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCtaa                                 >  2:1385266/1‑150 (MQ=255)
                                                                                                             aTTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCtaata                               <  2:229314/150‑1 (MQ=255)
                                                                                                                 cAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCt                           >  2:1388270/1‑150 (MQ=255)
                                                                                                                  aaaCAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTg                          <  1:1138715/150‑1 (MQ=255)
                                                                                                                  aaaCAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTg                          <  2:883028/150‑1 (MQ=255)
                                                                                                                   aaCAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGt                         <  1:360891/150‑1 (MQ=255)
                                                                                                                   aaCAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGt                         <  1:993253/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:1262527/150‑1 (MQ=255)
                                                                                                                               tggGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATg             >  2:913235/1‑150 (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:791662/150‑1 (MQ=255)
                                                                                                                                          aaaGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTgcc  >  2:774688/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
AACATGGCTTCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTGCC  >  NC_000913/3969289‑3969576

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: