Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,311,830 A→G 9.6% intergenic (+22/+18) tonB → / ← yciA membrane spanning protein in TonB‑ExbB‑ExbD transport complex/acyl‑CoA esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,311,8300AG 9.6% 139.5 / 4.8 52intergenic (+22/+18)tonB/yciAmembrane spanning protein in TonB‑ExbB‑ExbD transport complex/acyl‑CoA esterase
Reads supporting (aligned to +/- strand):  ref base A (24/23);  new base G (3/2);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

CTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAAT‑ATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTAT‑AGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCATAGCAGCACACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCG  >  NC_000913/1311537‑1312124
                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                        
cTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGTCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAAT‑ATCCTGTTTAAAATTAACGGCACACCCGAAATTCAGAAAGCAGAAAGTCAAAAGCCTCCGACCGGAgg                                                                                                                                                                                                                                                                                                   >  1:458178/1‑301 (MQ=255)
         cACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAAT‑ATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACt                                                                                                                                                                                                                                                                                         >  2:640855/1‑301 (MQ=255)
                       cgcttgtccgttctgtttCCGGCACTAGCACAGGCATTGCGCTTTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACTAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAAT‑ATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGt                                                                                                                                                                                                                                                                           <  1:491696/284‑1 (MQ=255)
                           aTCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAAT‑ATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAga                                                                                                                                                                                                                                                                       >  2:446675/1‑300 (MQ=255)
                               gCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAAT‑ATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGAACGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGcgcg                                                                                                                                                                                                                                                                   >  2:625839/1‑301 (MQ=255)
                                  gCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAAT‑ATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAgg                                                                                                                                                                                                                                                                <  2:408699/301‑1 (MQ=255)
                                     gTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAAT‑ATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTagatcggaagagca                                                                                                                                                                                                                                                             >  2:492342/1‑288 (MQ=255)
                                      aaTCCGGCACGGGCACAGGGATTGGGCATTTAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAAT‑ATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGtttt                                                                                                                                                                                                                                                            <  1:639978/300‑1 (MQ=255)
                                            gcaggggaagagggatgtcggaTTGAAGTGAAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAAGCCGCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAAT‑ATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCa                                                                                                                                                                                                                                                      <  1:538262/280‑1 (MQ=255)
                                            gCACGAGCACAGGCATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAAT‑ATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCa                                                                                                                                                                                                                                                      >  2:547964/1‑301 (MQ=255)
                                                         cATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAAT‑ATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACcgcgc                                                                                                                                                                                                                                         >  1:70994/1‑300 (MQ=255)
                                                             gcgcATTGAAGGGTATGTTAAAGGTAAATTTGACGTTACGCCGGATGGGTGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAAT‑ATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATa                                                                                                                                                                                                                                     <  1:396289/301‑1 (MQ=255)
                                                             gcgcATTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAAT‑ATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATa                                                                                                                                                                                                                                     >  2:408904/1‑301 (MQ=255)
                                                                   tGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAAT‑ATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTAAGGATCAACCGCGACATACTTaaa                                                                                                                                                                                                                               >  2:615732/1‑301 (MQ=255)
                                                                      aGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAAT‑ATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCACGACATACTTAaataa                                                                                                                                                                                                                            >  2:57037/1‑301 (MQ=255)
                                                                       ggtgtggTTAAAGTTCATGTTGGTGTTACGCCGGTTGGTGGCGTGGATACCGTACAAATCCGCTCAGCCAAGCCGGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGTATTGTGGTGAAT‑ATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAaataat                                                                                                                                                                                                                           <  1:429249/296‑1 (MQ=255)
                                                                                     taaaTTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAAT‑ATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACCTACTTAAATAATGCTTCTGTCGCtat                                                                                                                                                                                                             >  2:338958/1‑299 (MQ=255)
                                                                                         tttGACGTTACGCCGGATGGTCGCGTGGATAACGTAAAAATCCTTTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAAT‑ATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTAT‑Ag                                                                                                                                                                                                        <  1:408341/301‑1 (MQ=255)
                                                                                           tGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAAT‑ATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTAT‑AGCg                                                                                                                                                                                                      <  2:71084/301‑1 (MQ=255)
                                                                                            ggctttACGCGGGGTGGTCGCGGTGATAATGTACAAATCCTTTCATCCAAGTCTGGGAACATGTTTGCGCGTGAGGTGAATAATGCGATGCGCAGATGGCGTTAGGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAAT‑ATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGTTTCTGTCGCTTTAT‑AGCGt                                                                                                                                                                                                     <  1:614874/297‑1 (MQ=255)
                                                                                                 tACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGTCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAAT‑ATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTAT‑AGCGTTGccc                                                                                                                                                                                                <  2:458742/301‑1 (MQ=255)
                                                                                                               gcgTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAAT‑ATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTAT‑AGCGTTGCCCAATTGGTTCAGACg                                                                                                                                                                                  >  2:560103/1‑301 (MQ=255)
                                                                                                                cgTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAAT‑ATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑CCTATTACTCAACAGGTAAGGCACGAGGTTTTCATTCCGGATCAACCCAGCCATACTTAACTAATGCTTCTGTCGCTTTAT‑AGCATTTCCCAATTGCATCAGAAGc                                                                                                                                                                                 >  1:397730/1‑301 (MQ=255)
                                                                                                                                  ccggtctgCCAAGCCTTCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAAT‑ATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTAT‑AGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCacac                                                                                                                                                               <  1:628095/294‑1 (MQ=255)
                                                                                                                                    tctgagccaatgtgtctaaCAGGTTTGAGCGTGGGGTGAAAAATGCGATGCGGAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATGGTGGTGAAT‑ATCCGGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTAT‑AGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACtt                                                                                                                                                             <  1:329849/284‑1 (MQ=255)
                                                                                                                                       cAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAAT‑ATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTCAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTAT‑AGCGTTGCCCAATTGGTTCAGACGCTAATTTTTTAACCCAAACTCCaa                                                                                                                                                          >  1:35750/1‑301 (MQ=255)
                                                                                                                                         gCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAAT‑ATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTAT‑AGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCacacagatcgg                                                                                                                                                        >  2:628972/1‑294 (MQ=255)
                                                                                                                                                      aCATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAAT‑ATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTAT‑AGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGAcca                                                                                                                                           >  2:500428/1‑300 (MQ=255)
                                                                                                                                                         tGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGCTGGTGG‑GAATAACCCTGTTTAAAATAACCGCCACCCCCGAAATTCAGTAAGCAGAAAGTAAAAAGCC‑CCGACCCGAGGCCTTTGAACTATTATCCAAAAGTAAAGGCCCAAGGTTTCCCTCAGGGATCAACCCCGCCTTATTAAAATAAT‑CATCCTTCCCTTTTTAAGCTTTGCCAAATTGGTCCAGACACACCTTTTTCCAcaccacctccaaaatcaacccacccccaat                                                                                                                                        >  1:161324/1‑272 (MQ=255)
                                                                                                                                                          gTTTGGGCGTGATGTGAAAAATGGGGTGCGCAGATGGCGTTATGAGCCGGGTAAGACAGGCAGTGGGATTGTGGTGAAT‑ATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTAT‑AGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATgt                                                                                                                                       <  1:111492/301‑1 (MQ=255)
                                                                                                                                                                gCGTGAGGTGAAAAATGCGATGCGCATATGGGGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAAT‑ATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTAT‑AGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTccc                                                                                                                                 <  1:547291/301‑1 (MQ=255)
                                                                                                                                                                    gAGGTGAAAAATGCGATGCGCAGATGTCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGGAT‑ATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTAT‑AGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTc                                                                                                                             <  1:305986/301‑1 (MQ=255)
                                                                                                                                                                    aaGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAAT‑ATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTAT‑AGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTc                                                                                                                             <  1:559389/300‑1 (MQ=255)
                                                                                                                                                                                         aGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAAT‑ATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTAT‑AGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCATagcag                                                                                                        <  2:398244/301‑1 (MQ=255)
                                                                                                                                                                                                      aGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAAT‑ATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTAT‑AGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCATAGCAGCACACCACATCGc                                                                                           <  2:95909/301‑1 (MQ=255)
                                                                                                                                                                                                            ttatgCCAGGCAGTGGTATTTTGTTGAAT‑ATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGTGGCTTTTG‑ACTATGAATCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCACGACATACTTAAATAATGCTTCTGTCGCTTTAT‑AGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCATAGCAGCACACCACATCGCCGACCg                                                                                     <  1:56969/297‑1 (MQ=255)
                                                                                                                                                                                                                gtctggaaGTGGGAGTGTGGTTCAT‑ATGCTGTTTAAAATTAACGGGACCACCGAAATTCAGTAAGCAGAAAGTCGAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTAT‑AGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCATAGCAGCACACCACATCGCCGACCGCAAc                                                                                 <  1:181497/294‑1 (MQ=255)
                                                                                                                                                                                                                      cAGTGGGATTGTGGTGAAT‑ATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTAT‑AGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTATGGACACAGCGTGAATAGCAGCACACAACATCGCCGAACGAAACCggacc                                                                           >  2:334536/1‑298 (MQ=255)
                                                                                                                                                                                                                                gtggtgAAT‑ATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTAT‑AGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACCCACCATGCAAACAAGCACACCACATCACCGACCCCAACCGACCGAAACAAAgcc                                                                 >  1:469188/1‑299 (MQ=255)
                                                                                                                                                                                                                                                          cgatctccaccGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTGTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTAT‑AGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCATAGCAGCACACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTa                                        <  1:135730/295‑1 (MQ=255)
                                                                                                                                                                                                                                                               ctcgcgccttacctgttgagtaataGTCAAAAGCCTCCGGTCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTAT‑AGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCATAGCAGCACACCACATCGCCGACCGCAACCGCCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTAcgcgc                                   >  2:356447/25‑300 (MQ=255)
                                                                                                                                                                                                                                                                ccaccGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTAT‑AGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCATAGCAGCACACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTaagatcg                                  >  2:135855/1‑295 (MQ=255)
                                                                                                                                                                                                                                                                               tgagtaataGTCAAAAGCCTCCGGTCGGTGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACGTAAATAATGCTTCTGTCGCTTTAT‑AGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCATAGCAGCACACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCtt                   <  1:355963/293‑1 (MQ=255)
                                                                                                                                                                                                                                                                                aGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTAT‑AGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCATAGCAGCACACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCttt                  <  2:35801/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                aGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTAT‑AGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCATAGCAGCACACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTACGCAACCGTAGGCAATTTCttt                  >  2:383627/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                 tcttccgatctaaaGCCTCCGGTCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTAT‑AGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCATAGCAGCACACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCtttt                 <  2:164938/290‑1 (MQ=255)
                                                                                                                                                                                                                                                                                   aGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTAT‑AGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCATAGCAGCACACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGc               <  2:469789/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                          cAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGCTCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTAT‑AGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCCTAGCAGCACACCACATAGCAGACCGCAACCGGCCGACAGCAAGTCATTCCTTCAACCAGCACAGTCACTACGCCCACGTGGGCAATTTCTTTTCCCAGaaca        >  1:500018/1‑299 (MQ=255)
                                                                                                                                                                                                                                                                                          cAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTAT‑AGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAAGATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCATAGCAACACACCACACAGCCGACCGCAACCGACCGTCAGACACTCAATCCTTCAACCCACCCAGTCACAACGCCACCCGGCACAATTTCTTTTCCCataata        >  1:499997/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                            aaaGCCTCCGGTCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTAT‑AGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCATAGCAGCACACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAga      >  2:640656/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                            aaaGCCTCCGGTCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTAT‑AGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCATAGCAGCACACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGCCACTAACCGACCGTGGGCAATTTCttttagatcggaaga      >  1:164786/1‑290 (MQ=255)
                                                                                                                                                                                                                                                                                                ccTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTAT‑AGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCATAGCAGCACACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGcgccc  >  2:236299/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                        
CTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAAT‑ATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTG‑ACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTAT‑AGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCATAGCAGCACACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCG  >  NC_000913/1311537‑1312124

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑