Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,000,474 A→G 8.4% D420G (GAT→GGT)  elfC → putative outer membrane fimbrial subunit export usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,000,4740AG 8.4% 207.6 / 7.9 71D420G (GAT→GGT) elfCputative outer membrane fimbrial subunit export usher protein
Reads supporting (aligned to +/- strand):  ref base A (34/31);  new base G (3/3);  total (37/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 = 1.00e+00

AGAACCAATAGTAATGAACAGCAGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATAAAGGTCAGTCATATCGTTTTCTGTATGCCAAAACGCTCAACCACTTGGGCACTAACTTTCAATTGATGGGCTATCGCTATTC  >  NC_000913/1000326‑1000621
                                                                                                                                                    |                                                                                                                                                   
aGAACCAATAGTAATGAACAGCAGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGa                                                                                                                                                     >  7:285515/1‑149 (MQ=255)
   aCCAATAGTAATGAACAGCAGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGGtt                                                                                                                                                   >  2:398738/1‑148 (MQ=255)
   aCCAATAGTAATGAACAGCAGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGAttt                                                                                                                                                  <  1:178585/149‑1 (MQ=255)
     cAATAGTAATGAACAGCAGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGGTTTCg                                                                                                                                                <  1:398738/149‑1 (MQ=255)
             aTGAACAGCAGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACTACCTGGTATGGTGGTGGACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAAtt                                                                                                                                        <  5:173071/149‑1 (MQ=255)
             aTGAACAGCAGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAAtt                                                                                                                                        >  6:6169/1‑149 (MQ=255)
                aaCAGCAGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATttc                                                                                                                                      <  4:180110/148‑1 (MQ=255)
                  cagcagGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATttcgtt                                                                                                                                   <  8:3072/149‑1 (MQ=255)
                      agGAGAGCAAATTTGCCCAGGCCACGTGGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGAt                                                                                                                               <  2:180871/149‑1 (MQ=255)
                           agCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGAc                                                                                                                          <  8:285515/149‑1 (MQ=255)
                             cAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATttcgtt                                                                                                                                   >  8:67494/1‑138 (MQ=255)
                             cAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATttcgtt                                                                                                                                   <  7:67494/138‑1 (MQ=255)
                              aaaTTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCa                                                                                                                       <  5:6169/149‑1 (MQ=255)
                               aaTTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTACAGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAg                                                                                                                      >  4:263836/1‑149 (MQ=255)
                                     cccAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCg                                                                                                                   >  1:363892/1‑146 (MQ=255)
                                     cccAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGaa                                                                                                                 <  3:160802/148‑1 (MQ=255)
                                        aGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAagag                                                                                                              >  5:385552/1‑148 (MQ=255)
                                               gTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAg                                                                                                      <  7:393244/149‑1 (MQ=255)
                                               gTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAg                                                                                                      <  7:254893/149‑1 (MQ=255)
                                               gTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAg                                                                                                      <  5:288367/149‑1 (MQ=255)
                                                      tGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGGTTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATc                                                                                                >  2:195067/1‑148 (MQ=255)
                                                        gggCGGACCGGGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCaaa                                                                                             <  3:211678/149‑1 (MQ=255)
                                                          gCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGGTTTCGGAGCAATTTCGTTCGATGCGACCCa                                                                                                                       >  3:293416/1‑121 (MQ=255)
                                                          gCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGGTTTCGGAGCAATTTCGTTCGATGCGACCCa                                                                                                                       <  4:293416/121‑1 (MQ=255)
                                                             gACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGaa                                                                                        >  5:173381/1‑149 (MQ=255)
                                                                  tGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATaa                                                                                   <  2:363892/149‑1 (MQ=255)
                                                                  tGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATa                                                                                    >  5:187886/1‑148 (MQ=255)
                                                                   ggggAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATaaa                                                                                  <  3:208820/149‑1 (MQ=255)
                                                                      gAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGAt                                                                                                 >  5:337563/1‑131 (MQ=255)
                                                                      gAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGAt                                                                                                 <  6:337563/131‑1 (MQ=255)
                                                                            ccTGGTATGGTGGTGGACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATAAAGgtcagtca                                                                         >  6:172600/1‑149 (MQ=255)
                                                                             cTGGTATGGTGGTGGACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATAAAGgtcagtca                                                                         >  1:130395/1‑148 (MQ=255)
                                                                              tGGTATGGTGGTGGACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATAAAGGTCAGTCata                                                                       >  3:263721/1‑149 (MQ=255)
                                                                                    ggtggtggACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATAAAGGTCAGTCATATCGttt                                                                 >  4:85986/1‑149 (MQ=255)
                                                                                     gtggtggACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGTCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATAAAGGTCAGTCATATCGttt                                                                 <  2:130395/148‑1 (MQ=255)
                                                                                     gtggtggACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATAAAGGTCAGTCATATCGtttt                                                                >  2:383586/1‑149 (MQ=255)
                                                                                     gtggtggACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATAAAGGTCAGTCATATCGttt                                                                 <  6:173381/148‑1 (MQ=255)
                                                                                        gtggACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATAAAGGTCAGTCATATCGTTTTCTg                                                             >  6:77180/1‑149 (MQ=255)
                                                                                                 atGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGaa                                                                                                                 <  7:90849/88‑1 (MQ=255)
                                                                                                 atGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGaa                                                                                                                 >  8:90849/1‑88 (MQ=255)
                                                                                                 atGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATAAAGGTCAGTCATATCGTTTTCTGTATGCCaa                                                     >  6:6624/1‑148 (MQ=255)
                                                                                                  tGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATAAAGGTCAGTCATATCGTTTTCTGTATGCCaaaa                                                   >  3:276972/1‑149 (MQ=255)
                                                                                                   gCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATAAAg                                                                                 <  3:94432/118‑1 (MQ=255)
                                                                                                   gCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATAAAg                                                                                 >  4:94432/1‑118 (MQ=255)
                                                                                                    cTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATAAAGgtcagtc                                                                          >  3:138245/1‑124 (MQ=255)
                                                                                                    cTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATAAAGgtcagtc                                                                          <  4:138245/124‑1 (MQ=255)
                                                                                                    cTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATAAAGGTCAGTCATATCGTTTTCTGTATGCCAAAACg                                                 >  8:382144/1‑149 (MQ=255)
                                                                                                     tGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTACCCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATAAAGGTCAGTCATATCGTTTTCTGTATGCCAAAACGc                                                >  5:34111/1‑149 (MQ=255)
                                                                                                     tGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATAAAGGTCAGTCATATCGTTTTCTGTATGCCAAAACGc                                                >  7:376106/1‑149 (MQ=255)
                                                                                                       aaTATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATAAAGGTCAGTCATATCGTTTTCTGTATGCCAAAACGCt                                               >  4:221796/1‑148 (MQ=255)
                                                                                                              ccGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATc                                                                                                <  7:33615/92‑1 (MQ=255)
                                                                                                              ccGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATc                                                                                                >  8:33615/1‑92 (MQ=255)
                                                                                                              ccGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATAAAGGTCAGTCATATCGTTTTCTGTATGCCAAAACGCTCAACCACt                                       <  5:172600/149‑1 (MQ=255)
                                                                                                               cGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATAAAGGTCAGTCATATCGTTTTCTGTATGCCAAAACGCTCAACCACtt                                      <  6:385552/149‑1 (MQ=255)
                                                                                                                   ccgccATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATAAAGGTCAGTCATATCGTTTTCTGTATGCCAAAACGCTCAACCACTTggg                                   <  5:77180/148‑1 (MQ=255)
                                                                                                                      ccATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATa                                                                                    <  6:375081/96‑1 (MQ=255)
                                                                                                                      ccATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATa                                                                                    >  5:375081/1‑96 (MQ=255)
                                                                                                                      ccATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATAAAGGTCAGTCATATCGTTTTCTGTATGCCAAAACGCTCAACCACTTGGGCAc                                >  4:181022/1‑148 (MQ=255)
                                                                                                                        aTGTTTGGTCTGGGATTTAACCTCGGCGGTTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATAAAGGTCAGTCATATCGTTTTCTGTATGCCAAAACGCTCAACCACTTGGGCACTaa                             <  1:195067/149‑1 (MQ=255)
                                                                                                                             tGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATAAAGGTCAGTCATATCGTTTTCTGTATGCCAAAACGCTCCACCACTTGGGCACTAACttt                         >  8:391588/1‑148 (MQ=255)
                                                                                                                                  tGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGACCATAAAGGTCAGTCATATCGTTTTCTGTATGCCAAAACGCt                                               <  5:333514/121‑1 (MQ=255)
                                                                                                                                  tGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATAAAGGTCAGTCATATCGTTTTCTGTATGCCAAAACGCt                                               >  6:333514/1‑121 (MQ=255)
                                                                                                                                  tGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATAAAGGTCAGTCATATCGTTTTCTGTATGCCAAAACGCTCAACCACTTGGGCACTAACTTTCAAtt                    >  7:68035/1‑148 (MQ=255)
                                                                                                                                  tGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATAAAGGTCAGTCATATCGTTTTCTGTATGCCAAAACGCTAAACCACTTGGGCACTAACTTTCAAtt                    >  1:190841/1‑148 (MQ=255)
                                                                                                                                      aTTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATAAAGGTCAGTCATATCGTTTTCTGTATGCCAAAACGCTCAACCACTTGGGCACTAACTTTCAATTGATgg               <  1:173353/149‑1 (MQ=255)
                                                                                                                                       tttAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATAAAGGTCAGTCATATCGTTTTCTGTATGCCAAAACGCTCAACCACTTGGGCACTAACTTTCAATTGATggg              >  6:1380/1‑149 (MQ=255)
                                                                                                                                            ccTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATAAAGGTCAGTCATATCGTTTTCTGTATGCCAAAACGCTc                                              <  3:285220/112‑1 (MQ=255)
                                                                                                                                            ccTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATAAAGGTCAGTCATATCGTTTTCTGTATGCCAAAACGCTc                                              >  4:285220/1‑112 (MQ=255)
                                                                                                                                            ccTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATAAAGGTCAGTCATATCGTTTTCTGTATGCCAAAACGCTCACCCACTTGGGCACTAACTTTCAATTGATGGGCTATc         >  4:129583/1‑149 (MQ=255)
                                                                                                                                              tCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATAAAGGTCAGTCATATCGTTTTCTGTATGCCAAAACGCTCAACCACTTGGGCACTAACTTTCaa                      <  4:65338/134‑1 (MQ=255)
                                                                                                                                              tCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATAAAGGTCAGTCATATCGTTTTCTGTATGCCAAAACGCTCAACCAATTGGGCACTAACTTTCaa                      >  3:65338/1‑134 (MQ=255)
                                                                                                                                                   gATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATAAAGGTCAGTCATATCGTTTTCTGTATGCCAAAACGCTCAACCACTTGGGCACTAACTTTCAATTGATGGGCTATCGCTATTc  <  7:382144/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGAACCAATAGTAATGAACAGCAGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTACCGTGCCGCCATGTTTGGTCTGGGATTTAACCTCGGCGATTTCGGAGCAATTTCGTTCGATGCGACCCAGGCGAAGAGTACGCTAGCAGATCAAAGCGAACATAAAGGTCAGTCATATCGTTTTCTGTATGCCAAAACGCTCAACCACTTGGGCACTAACTTTCAATTGATGGGCTATCGCTATTC  >  NC_000913/1000326‑1000621

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