Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,969,441 G→T 6.4% A137A (GCG→GCT wzzE → Entobacterial Common Antigen (ECA) polysaccharide chain length modulation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,969,4410GT 6.4% 203.9 / 3.2 78A137A (GCG→GCTwzzEEntobacterial Common Antigen (ECA) polysaccharide chain length modulation protein
Reads supporting (aligned to +/- strand):  ref base G (41/32);  new base T (3/2);  total (44/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-01

TGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTGCCAGCCAGCGTGC  >  NC_000913/3969300‑3969587
                                                                                                                                             |                                                                                                                                                  
tGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTgg                                                                                                                                             >  2:843902/1‑149 (MQ=255)
tGCCGACCAACCAGCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAACGCCGCATCGGCTTTGCTgt                                                                                                                                             <  2:294086/149‑2 (MQ=255)
 gCCGACCAACCATCGCTCATGGACGAAGCCTATAAAGAGGTTGTTATGCAACTGGCCTCGAGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAACCAGCGGATGGTGGGCAACAGCAAAGCCAAGGCGGCGTTGCTGGa                                                                                                                                            >  2:1042441/1‑149 (MQ=255)
  ccGCCCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCTACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTAGAAACAGCGGAGGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTggct                                                                                                                                           >  2:85607/1‑147 (MQ=255)
   cggccacccaTCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCACAGCCGATGCGGCGTTGCTGGATg                                                                                                                                          >  1:907579/4‑149 (MQ=255)
   cGACCAACCTTCGGTCCTGGACGAAGCCTATAAAGAGTTTGTTATGAAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATg                                                                                                                                          <  2:703913/149‑1 (MQ=255)
      ccaaccaTCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGaaa                                                                                                                                       <  1:406198/149‑1 (MQ=255)
         accaTCGGTCATGGACGAAGCCTATAAAGAGTTTGTGATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCTGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGa                                                                                                                                    <  2:447362/149‑1 (MQ=255)
          ccaTCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGAt                                                                                                                                   <  2:833896/149‑1 (MQ=255)
          ccaTCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGAt                                                                                                                                   >  1:793053/1‑149 (MQ=255)
             tCGGTCATGGACGAAGCCTATAAAGAGTTTGATATGCAACTGGCCTCGTCGGATACCCGCAGAGAGTTCTGGCACCAAAACCACTATTACAAAAAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGATGAAGGATGAAATGATTaa                                                                                                                                >  2:890738/1‑149 (MQ=255)
              cGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTaac                                                                                                                               <  1:751443/149‑1 (MQ=255)
                gTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGCCTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTtacac                                                                                                                             >  2:173691/1‑148 (MQ=255)
                gTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTaacaa                                                                                                                             >  2:1165081/1‑149 (MQ=255)
                 tCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTaacaac                                                                                                                            <  2:166726/149‑1 (MQ=255)
                 tCATGGACGAAGACTATAAAGATTTTGTTATGCAACTGTCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTaacaac                                                                                                                            <  1:251991/149‑1 (MQ=255)
                     ggACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATcc                                                                                                                        >  1:79358/1‑149 (MQ=255)
                     ggACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATcc                                                                                                                        >  2:913356/1‑149 (MQ=255)
                       aCGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGCTTAACAACATCCAg                                                                                                                      >  2:167655/1‑149 (MQ=255)
                       aCGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAg                                                                                                                      <  2:907579/149‑1 (MQ=255)
                         gAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGtt                                                                                                                    >  2:1142943/1‑149 (MQ=255)
                          aaGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGttt                                                                                                                   <  1:373437/149‑1 (MQ=255)
                                ataAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATcccc                                                                                                             >  2:340417/1‑149 (MQ=255)
                                ataAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATcccc                                                                                                             >  2:541351/1‑149 (MQ=255)
                                 aaaaGAGTTGGTTATGCAACTGCCCTCGGGGGCTACCCACAGAGAGTTCTGGCTGCAAACCGACAAATACAAACAGCGGAAGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCg                                                                                                            <  2:1059047/148‑1 (MQ=255)
                                     gagTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGgaga                                                                                                        >  1:972868/1‑149 (MQ=255)
                                        tttGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACtt                                                                                                     <  1:483538/149‑1 (MQ=255)
                                              aTGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCcg                                                                                               >  2:53126/1‑149 (MQ=255)
                                                gCAACTGGCCTCGTGGGCTACCCGCAGAGAGTTCTGGCTGCTAACCGACTATTAGAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACATCATCCAGGTTATCCCCGGAGACTTTACCcgcg                                                                                             >  2:603393/1‑149 (MQ=255)
                                                  tactgtgcTCGTGTGAGACCCTCAGAGAGTTCTTGCTCCAAACCCCCTATTACAAACAGCGGCTGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGACATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCcgcgcg                                                                                           <  2:771738/142‑1 (MQ=255)
                                                  aaCTGGCCTCGTGGGATACCCTCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACCGCGGAGGGGGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTCTCCCCGGAGACTTTACCcgcgcg                                                                                           <  1:565571/149‑1 (MQ=255)
                                                  aaCTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCcgcgcg                                                                                           >  1:358882/1‑149 (MQ=255)
                                                  aaCTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCcgcgcg                                                                                           >  1:757310/1‑149 (MQ=255)
                                                     tGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTc                                                                                        >  1:7409/1‑149 (MQ=255)
                                                       gCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCaa                                                                                      >  1:925984/1‑149 (MQ=255)
                                                         cTCGTGGGAGACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATg                                                                                    >  1:157813/1‑149 (MQ=255)
                                                          tCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCCGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGCAGACTTTACCCGCGCGGTCAATGa                                                                                   >  2:1047051/1‑149 (MQ=255)
                                                           cgggggATACCCGCAGAGAGATCTGGCTCCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGAc                                                                                  <  1:1144589/146‑1 (MQ=255)
                                                                gATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGt                                                                             <  1:843902/149‑1 (MQ=255)
                                                                     ccGCAGAGTGTTCTGGCTGCAAACCTTCTATTACAAACAGCGGATCGTGGGCTACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCTGTTTATCCCCGGAGACTTAACCCGCGCGGTCAATGACAGCGTGAAGc                                                                        <  1:890738/149‑1 (MQ=255)
                                                                     ccGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACCTCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGc                                                                        >  2:1157214/1‑149 (MQ=255)
                                                                           agagTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTg                                                                  >  1:397621/1‑149 (MQ=255)
                                                                            gagTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGc                                                                 >  2:91471/1‑149 (MQ=255)
                                                                            gagTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGATGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGc                                                                 >  2:1115392/1‑149 (MQ=255)
                                                                             agTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGACGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGCCTTTACCCGCGCGGTCAATGACAGCGTGAAGCGTATTGcc                                                                >  1:529009/1‑149 (MQ=255)
                                                                               ttCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGAGGCGGCGTTGCTGGATGAAATGATTAACAACATCCTGTTTATCACAGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGa                                                              >  2:1136657/1‑149 (MQ=255)
                                                                               ttCTGGCTGCAAACCGACTATTACAAACAGCGGATGGGGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGa                                                              <  1:793963/149‑1 (MQ=255)
                                                                                    gCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACcg                                                         >  1:862251/1‑149 (MQ=255)
                                                                                     cTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACcgc                                                        <  1:446641/149‑1 (MQ=255)
                                                                                     cTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACcgc                                                        <  1:541351/149‑1 (MQ=255)
                                                                                     cTGCAAACCGACTATTACAAACACCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACcgc                                                        <  1:305923/149‑1 (MQ=255)
                                                                                       gCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCgcgc                                                      >  1:530572/1‑149 (MQ=255)
                                                                                          aaCCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAACTGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTg                                                   <  2:793053/149‑1 (MQ=255)
                                                                                            acGACTATTACATACAGCTGATGGTGGGAAACATCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGAc                                                 <  2:241935/148‑1 (MQ=255)
                                                                                              gACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGAGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGc                                               <  1:923526/149‑1 (MQ=255)
                                                                                                cTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCta                                             <  2:972868/149‑1 (MQ=255)
                                                                                                 tATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCACCGGAGACTTTACCCGCGCGGTCATTGACAGCGGGAAGCTTATTGCCGAAACCGCGCCTGACGCtaa                                            >  1:89251/1‑149 (MQ=255)
                                                                                                    tACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCtaataa                                         <  1:53126/149‑1 (MQ=255)
                                                                                                         aCAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGt                                    >  1:116177/1‑149 (MQ=255)
                                                                                                           aGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTa                                  <  1:657637/149‑1 (MQ=255)
                                                                                                            gCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAAGAACATCCAGTTTATCCCCGGAGACTTTACACGCGCGGTCAATGACAGCGTGACGCTTATTGCCGAATCCGCGCCTGACGCTAATAACCTGTTAc                                 >  1:535926/1‑149 (MQ=255)
                                                                                                             cGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGCCGCTAATAACCTGTTACg                                >  2:95380/1‑149 (MQ=255)
                                                                                                                     gggCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATg                        <  2:79358/149‑1 (MQ=255)
                                                                                                                     gggCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTGATCCCCGGAGACGTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCCAAACCGCGCCTGACGCTACTAACCTGTTACGTCAGTATg                        >  1:792374/1‑149 (MQ=255)
                                                                                                                      ggCAACAGCAAAGACGATGCGGCGTTGCTGGATGAAATGATTAAAAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGATTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGACAGTATGt                       >  2:910079/1‑149 (MQ=255)
                                                                                                                        cAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTg                     >  1:30335/1‑149 (MQ=255)
                                                                                                                        atccAGCAACGCCGCATCGGCTTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGaaa                                                            >  1:343627/4‑110 (MQ=38)
                                                                                                                        atccAGCAACGCCGCATCGGCTTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGaaa                                                            <  2:343627/107‑1 (MQ=38)
                                                                                                                         tccAGCAACGCCGCATCGGCTTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGCCGCTAATAACCTGTTACGTCAGTATGTTGc                    >  1:1066716/3‑149 (MQ=255)
                                                                                                                         aaCAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGc                    <  1:1165081/149‑1 (MQ=255)
                                                                                                                           cAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCtt                  <  1:120630/149‑1 (MQ=255)
                                                                                                                            aGCAACGCCGCATCGGCTTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCttt                 >  2:271230/1‑149 (MQ=255)
                                                                                                                            aGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCttt                 >  1:991629/1‑149 (MQ=255)
                                                                                                                             gCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCtttt                <  1:481330/149‑1 (MQ=255)
                                                                                                                              cAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTg               >  2:1011588/1‑149 (MQ=255)
                                                                                                                               aaaGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTgc              <  1:913356/149‑1 (MQ=255)
                                                                                                                                aaGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTgcc             >  1:433321/1‑149 (MQ=255)
                                                                                                                                 aGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTgcca            <  2:733685/149‑1 (MQ=255)
                                                                                                                                          ggcgTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTGCCAGCCAGCGTg   >  2:257170/1‑149 (MQ=255)
                                                                                                                                           gcgTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTGCCAGCCAGCGTGc  >  1:800028/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
TGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTGCCAGCCAGCGTGC  >  NC_000913/3969300‑3969587

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