Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,311,285 A→C 7.1% E66A (GAG→GCG) ‡ tonB → membrane spanning protein in TonB‑ExbB‑ExbD transport complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,311,2850AC 7.1% 243.4 / 9.3 85E66A (GAG→GCG) ‡tonBmembrane spanning protein in TonB‑ExbB‑ExbD transport complex
Reads supporting (aligned to +/- strand):  ref base A (36/43);  new base C (3/3);  total (39/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.57e-01

GGTCTGCATTCATGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAAAACGCGATGTCAAAC  >  NC_000913/1311136‑1311434
                                                                                                                                                     |                                                                                                                                                     
ggTCTGCATTCATGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGa                                                                                                                                                       >  2:1591047/1‑150 (MQ=255)
   cTGCATTCATGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCGGCGCAGCCGATTTCTGTCACGATGGTTACGCCTACTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGcc                                                                                                                                                    <  1:324034/150‑1 (MQ=255)
    tGCATTCATGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCCCAAGCCGTTCAGCCGCCCCCGGAGCCg                                                                                                                                                   >  2:288464/1‑150 (MQ=255)
     gCATTCATGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCgg                                                                                                                                                  >  2:876458/1‑150 (MQ=255)
      cATTCATGGTGCTTTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCggt                                                                                                                                                 <  2:1309200/150‑1 (MQ=255)
      cATTCATGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACCAGCCGTTCAGCCGCCACCGGAGCCggt                                                                                                                                                 >  2:138328/1‑150 (MQ=255)
      cATTCATGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCggt                                                                                                                                                 <  1:862536/150‑1 (MQ=255)
      cATTCATGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCggt                                                                                                                                                 <  2:103008/150‑1 (MQ=255)
         tCATGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGCTCCggtggc                                                                                                                                              <  2:854704/150‑2 (MQ=255)
           aTGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTag                                                                                                                                            <  1:967694/150‑1 (MQ=255)
           aTGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTag                                                                                                                                            <  1:66657/150‑1 (MQ=255)
            tGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGCTCCggtggcggc                                                                                                                                           <  1:424013/150‑5 (MQ=255)
                gCTGTTGTCGCGGGGCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTGCAGCCGCCACCGGAGCCGGTGGTAGAGCCa                                                                                                                                       <  1:288464/150‑1 (MQ=255)
                 cTGTTGTCGCGGGTCTGCTCTATACCTCCGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGAGGTAGAGCCAg                                                                                                                                      >  1:273210/1‑150 (MQ=255)
                   gttgtGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGaa                                                                                                                                    <  1:495322/150‑1 (MQ=255)
                    gtgtgGCGGGTCTGCGCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAAc                                                                                                                                   <  1:136013/149‑1 (MQ=255)
                     tgtgGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAAcc                                                                                                                                  <  1:400008/150‑1 (MQ=255)
                     tgtgGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAAcc                                                                                                                                  <  2:404560/150‑1 (MQ=255)
                      gtgGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACCAGCCGTTCAGCCGCCCCCGGCGCCGGTGGTAGCGCCAGAACCg                                                                                                                                 >  2:1241930/1‑150 (MQ=255)
                      gtgGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCg                                                                                                                                 >  2:1470037/1‑150 (MQ=255)
                         gCGGGTCTGCGCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCACCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGaa                                                                                                                              <  2:434996/150‑1 (MQ=255)
                             gTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTg                                                                                                                          >  1:104729/1‑150 (MQ=255)
                               cTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAg                                                                                                                        <  1:1591047/150‑1 (MQ=255)
                               cTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAg                                                                                                                        >  1:1424461/1‑150 (MQ=255)
                                tGCTCTATACCTCGGTACATCAGGGTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGc                                                                                                                       <  1:1214211/150‑1 (MQ=255)
                                      ataCCTCGTTACTTCAGTTTATTGAACTTCCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCAGTTCAGCGGCCTCCGGAGCCGGTGGTAGATCCAGAACCGGAACCTGAGCCGATcc                                                                                                                 <  2:665352/150‑1 (MQ=255)
                                      ataCCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCAGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATcc                                                                                                                 <  2:81790/150‑1 (MQ=255)
                                         ccTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACACGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGCACCGGAACCTGAGCCGATCCCCg                                                                                                              >  2:779288/1‑150 (MQ=255)
                                         ccTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGCTCCggtggcggct                                                                                                                                          >  2:1182342/1‑117 (MQ=255)
                                         ccTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGCTCCggtggcggct                                                                                                                                          <  1:1182342/122‑6 (MQ=255)
                                         ccTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCg                                                                                                              >  2:808120/1‑150 (MQ=255)
                                            cGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGACCCGGAACCTGAGCCGATCCCCGAAc                                                                                                           <  1:481601/150‑1 (MQ=255)
                                                     aGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGCGCCTGAGCCGATCCCCGCACCGCCAAaag                                                                                                  >  1:1169403/1‑150 (MQ=255)
                                                     aGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAaag                                                                                                  >  1:1433409/1‑150 (MQ=255)
                                                     aGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAAACGCCACAAGCCGTTCAGCCGCCACAGGAGCCGGTGGTAGAGCCAGAACCTGAACCTGAGCCGATCCCCGCACCGCCAAaag                                                                                                  >  1:1219268/1‑150 (MQ=255)
                                                        ttATTGAACTACCTGCGACTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAaagaag                                                                                               <  2:484551/150‑1 (MQ=255)
                                                         tATTGAACTACCTGCGCCTGCGCAGCCGAGTTCTGTCACGATGGTGACGCCTGCTGATCTCGCACCGCCCCACGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGCGCCAGCACCGGCACCTGAGCCGCTCCCCGCACCGCCCATAGAAGc                                                                                              >  1:136545/1‑150 (MQ=255)
                                                            tGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACACGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCAcc                                                                                           >  1:453673/1‑150 (MQ=255)
                                                                  aCCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCTCGAACCGCCAAAAGAAGCACCggtggt                                                                                     >  2:1198738/1‑150 (MQ=255)
                                                                   cctgcgcctgcgCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTc                                                                                    >  1:1356582/1‑150 (MQ=255)
                                                                    ctgcgcctgcgcAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCa                                                                                   <  1:530136/150‑1 (MQ=255)
                                                                       cgcctgcgcAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGCACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGCACCGCCAAAAGCAGCACCGGTGGTCATTg                                                                                >  2:1164937/1‑150 (MQ=255)
                                                                        ccctgcgcAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGa                                                                               <  1:769271/149‑1 (MQ=255)
                                                                         cctgcgcAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGaa                                                                              >  1:1536929/1‑150 (MQ=255)
                                                                            gcgcAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAaag                                                                           >  2:1187486/1‑150 (MQ=255)
                                                                              gcAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTTGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAaagcc                                                                         <  1:438147/150‑1 (MQ=255)
                                                                              gcAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAaagcc                                                                         >  2:1438519/1‑150 (MQ=255)
                                                                              gcAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAaagcc                                                                         <  2:1433409/150‑1 (MQ=255)
                                                                               cAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAaagccg                                                                        >  2:1528135/1‑150 (MQ=255)
                                                                                 gCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAaagccgaa                                                                      >  1:674969/1‑150 (MQ=255)
                                                                                  ccGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAaagccgaag                                                                     <  1:1427540/150‑1 (MQ=255)
                                                                                   cGATTTCTGTCACGATGGTTACGCCTGCTGAACTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAaagccgaagc                                                                    <  2:945292/150‑1 (MQ=255)
                                                                                   cGATTTCTGGCACGAGGGTTACGCCTGCTTATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAaagccgaagc                                                                    <  1:1582569/150‑1 (MQ=255)
                                                                                         cTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAAc                                                              >  1:1092598/1‑150 (MQ=255)
                                                                                            tCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGACGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTa                                                           <  1:553941/150‑1 (MQ=255)
                                                                                              ccGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTATAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTatg                                                         <  2:1639252/149‑3 (MQ=255)
                                                                                               cGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGc                                                        <  1:122269/150‑1 (MQ=255)
                                                                                                   ggTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCaaa                                                    <  1:1187486/150‑1 (MQ=255)
                                                                                                      tACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAAcc                                                 <  1:1262762/150‑1 (MQ=255)
                                                                                                       aCGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCg                                                >  1:739114/1‑150 (MQ=255)
                                                                                                        cGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGCACCGGACCCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCAGTGAAAAGCCGAAGCCGAAGCCGCAGCCAAACCCGa                                               >  2:94706/1‑150 (MQ=255)
                                                                                                            tgctgATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAGCCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTACGCCAAAACCGTAGcc                                           >  2:93670/1‑150 (MQ=255)
                                                                                                            tgctgATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGcc                                           <  2:1092598/150‑1 (MQ=255)
                                                                                                             gctgATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCg                                          >  2:1516314/1‑150 (MQ=255)
                                                                                                                gATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTg                                       <  2:674969/150‑1 (MQ=255)
                                                                                                                 aTCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGa                                      <  1:1470037/150‑1 (MQ=255)
                                                                                                                     cGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGaaaaa                                  <  2:1698852/150‑1 (MQ=255)
                                                                                                                     cGAACCGCCACAAGCCGTGCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGACCCTGAGCCGAGCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGaaaaa                                  <  2:1019299/150‑1 (MQ=255)
                                                                                                                      gAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAACCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAg                                 >  1:363775/1‑150 (MQ=255)
                                                                                                                        aCCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGt                               <  1:418475/150‑1 (MQ=255)
                                                                                                                           gCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACa                            <  2:273210/150‑1 (MQ=255)
                                                                                                                            ccACAAGCCGTACAGCCGCCACCGGAGCCGGTGGTAGCGCCAGAACCCGAACCTGAGCCGATCCCCGACCCGCCAAAAGAAGCCCCGGTGGTCATTGACAAGCCGCAGCCGCAACCTCAGCCACAACCGCAGCCGGTGAAACAGGTACAg                           >  1:133275/1‑150 (MQ=255)
                                                                                                                              acaAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGaaaaa                                  >  1:1504017/1‑141 (MQ=255)
                                                                                                                              acaAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGaaaaa                                  <  2:1504017/141‑1 (MQ=255)
                                                                                                                                   ccGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGagcagc                    <  1:525426/150‑1 (MQ=255)
                                                                                                                                    cGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGcc                   >  2:1681590/1‑150 (MQ=255)
                                                                                                                                    cGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGcc                   >  1:1387120/1‑150 (MQ=255)
                                                                                                                                        gctctaccACCGGCTCCGGTGGTAGAGCCAGAACCGGCACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCAGTGAAAAGCCGAAGCCGAAGCCTAAGCCAAAACCGAAGCCGGAGAAAAAGGTACAGGAGCAGCCataa               >  1:36912/7‑150 (MQ=255)
                                                                                                                                         aGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGCACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAAAAc              >  1:1009064/1‑150 (MQ=255)
                                                                                                                                         aGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAAAAc              >  2:739018/1‑150 (MQ=255)
                                                                                                                                          gccgccACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGCAGCCGGAGAAAAAGGTACAGGAGCAGCCAAACcg             >  2:1302973/1‑150 (MQ=255)
                                                                                                                                            cgccACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAAAAcgcg           <  2:1536929/150‑1 (MQ=255)
                                                                                                                                              ccACCGGAGCCGGTGGTAGAGCCAGAACCGGCACCTTAGCCGATCCCCGAAACGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCAGAATCCGAAACCTAAGCCAAAACCGACGCCGGTGAATAAGGTACAGGAGCCGCCAAACCGCGag         >  2:220190/1‑149 (MQ=255)
                                                                                                                                              caaCCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAAAACGCGAt         <  1:753226/148‑1 (MQ=255)
                                                                                                                                               cACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAAAACGCGATg        <  2:1568290/150‑1 (MQ=255)
                                                                                                                                                 ccGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAAAACGCGATGTc      <  2:453673/150‑1 (MQ=255)
                                                                                                                                                     aGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAAAACGCGATGTCAAAc  <  1:1528135/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGTCTGCATTCATGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAAAACGCGATGTCAAAC  >  NC_000913/1311136‑1311434

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