Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,969,429 A→C 6.5% K133N (AAA→AAC 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,4290AC 6.5% 189.0 / 4.3 77K133N (AAA→AACwzzEEntobacterial Common Antigen (ECA) polysaccharide chain length modulation protein
Reads supporting (aligned to +/- strand):  ref base A (41/31);  new base C (3/2);  total (44/33)
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.97e-01

ACATGGCTTCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTGCCA  >  NC_000913/3969290‑3969577
                                                                                                                                           |                                                                                                                                                    
aCATGGCGTCTGACTTCAAACCCTCGGTCATGGACGAAGCCTATAAAGAGTTTGTAATGCAACTGGCCTCGTGGGATCCCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATgcg                                                                                                                                             <  2:504708/149‑1 (MQ=255)
aCATGGCATCTGCCGACCAACCATCGGTCATGGAAGAAGCCTATAAAGAGTTTGTTAAGCAACTGGCCTCCTGGGATACCCGCATAGAGTTCTGGCTGCAATCCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATgcg                                                                                                                                             >  2:923526/1‑149 (MQ=255)
          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:541351/149‑1 (MQ=255)
                                                                                               cTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACcgc                                              <  1:446641/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)
                                                                                                                                           |                                                                                                                                                    
ACATGGCTTCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTGCCA  >  NC_000913/3969290‑3969577

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: