Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,674,177 G→T 5.5% E89* (GAG→TAG) ‡ tqsA → autoinducer 2 exporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,674,1770GT 5.5% 305.6 / 5.8 91E89* (GAG→TAG) ‡tqsAautoinducer 2 exporter
Reads supporting (aligned to +/- strand):  ref base G (41/45);  new base T (3/2);  total (44/47)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.09e-01

TCGCATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCTCA  >  NC_000913/1674028‑1674324
                                                                                                                                                     |                                                                                                                                                    
tcGCATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACg                                                                                                                                                      <  1:519294/150‑1 (MQ=255)
tcGCATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACg                                                                                                                                                      <  2:79040/150‑1 (MQ=255)
tcGCATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACg                                                                                                                                                      >  1:958067/1‑150 (MQ=255)
  gCATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAg                                                                                                                                                    <  1:840063/150‑1 (MQ=255)
   cATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGt                                                                                                                                                   >  2:886435/1‑150 (MQ=255)
   cATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGt                                                                                                                                                   <  1:404352/150‑1 (MQ=255)
           attGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCgg                                                                                                                                          <  2:459399/150‑1 (MQ=255)
            ttGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGTCAACTCGTTGAGCGCGGa                                                                                                                                         <  1:548386/150‑1 (MQ=255)
                tGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGa                                                                                                                                         <  2:1392698/146‑1 (MQ=255)
                tGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGa                                                                                                                                         >  1:1392698/1‑146 (MQ=255)
                tGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGtt                                                                                                                                     <  1:1379606/150‑1 (MQ=255)
                    aTTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCg                                                                                                                                 >  1:164773/1‑150 (MQ=255)
                      tCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCa                                                                                                                               >  1:665416/1‑150 (MQ=255)
                       cTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCaa                                                                                                                              >  2:690971/1‑150 (MQ=255)
                            cccGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGTCAACTCGTTGagcgcg                                                                                                                                           <  1:1195762/132‑6 (MQ=255)
                            cccGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGTCAACTCGTTGagcgcg                                                                                                                                           >  2:1195762/1‑127 (MQ=255)
                             ccGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGc                                                                                                                        <  1:721134/150‑1 (MQ=255)
                             ccGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGc                                                                                                                        <  1:1180727/150‑1 (MQ=255)
                              ggCTGGTGCAACAGATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCa                                                                                                                       <  2:361534/149‑1 (MQ=255)
                                  ggTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAActc                                                                                                                   <  1:346642/150‑1 (MQ=255)
                                    tGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCta                                                                                                                 >  2:1121472/1‑150 (MQ=255)
                                     gCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCtat                                                                                                                >  2:511504/1‑150 (MQ=255)
                                       aaCACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCtatta                                                                                                              <  2:253728/150‑1 (MQ=255)
                                         cacaTGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATg                                                                                                            <  1:1078939/150‑1 (MQ=255)
                                         cacaTGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATg                                                                                                            >  1:1209164/1‑150 (MQ=255)
                                           caTGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGAc                                                                                                          <  1:952120/150‑1 (MQ=255)
                                           caTGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGAc                                                                                                          >  1:1438091/1‑150 (MQ=255)
                                            aTGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGTGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGAcg                                                                                                         <  2:209787/150‑1 (MQ=255)
                                             tGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGATATTAGCTTATATGGGTTCCGCGCTCAAAGAGTTCACGCGGACGTTACCGCAATATAGCAACTCTATTATGActc                                                                                                        >  2:122743/1‑148 (MQ=255)
                                              ggTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCGTCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGAcgcc                                                                                                       <  2:220759/150‑1 (MQ=255)
                                              ggTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGAcgcc                                                                                                       >  2:576417/1‑150 (MQ=255)
                                                tCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGAcgccgc                                                                                                     <  1:886435/150‑1 (MQ=255)
                                                 ccGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCt                                                                                                    >  1:1268207/1‑150 (MQ=255)
                                                   gCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGc                                                                                                  >  2:338909/1‑150 (MQ=255)
                                                      ggCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAg                                                                                               >  1:642976/1‑150 (MQ=255)
                                                        cGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAgcg                                                                                             <  1:584024/150‑1 (MQ=255)
                                                          ggtgCCGCGGGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGGGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGCGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCt                                                                                           <  2:1087003/149‑1 (MQ=255)
                                                            tgCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTg                                                                                         <  2:1438091/150‑1 (MQ=255)
                                                             gCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGa                                                                                        >  1:58340/1‑150 (MQ=255)
                                                             gCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGa                                                                                        >  2:727509/1‑150 (MQ=255)
                                                                  gtgtACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCgt                                                                                   <  1:578087/150‑1 (MQ=255)
                                                                   tgtACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCgtt                                                                                  >  1:738612/1‑150 (MQ=255)
                                                                    gtACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCgttg                                                                                 >  1:333552/1‑150 (MQ=255)
                                                                      aCTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCgttgtt                                                                               >  2:1022042/1‑150 (MQ=255)
                                                                      aCTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCgttgtt                                                                               <  1:191330/150‑1 (MQ=255)
                                                                      aCTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATATGAGTTCAGCGATCAACGAGTTGAAGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCgttgtt                                                                               >  2:1394850/1‑150 (MQ=255)
                                                                       cTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGttgttg                                                                              <  1:206222/150‑1 (MQ=255)
                                                                         ggCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCa                                                                            >  1:28398/1‑150 (MQ=255)
                                                                         ggCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCa                                                                            >  2:328810/1‑150 (MQ=255)
                                                                         ggCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCa                                                                            >  2:1110120/1‑150 (MQ=255)
                                                                           cGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCTATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAAc                                                                          <  2:164773/150‑1 (MQ=255)
                                                                             gtgtAGATTTTGATGACCATCATCGTGATGGCGATGGTGTTCCTATTAGCTTATATGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTCAACCGTTGTTGCAAcgc                                                                        >  1:1122824/1‑150 (MQ=255)
                                                                                 cGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAg                                                                    <  2:654462/150‑1 (MQ=255)
                                                                                 cGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAg                                                                    <  2:1244375/150‑1 (MQ=255)
                                                                                  gATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAgg                                                                   <  1:724688/150‑1 (MQ=255)
                                                                                   aTTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAggg                                                                  <  2:1136977/150‑1 (MQ=255)
                                                                                       tgatgaCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTg                                                              <  2:1202975/150‑1 (MQ=255)
                                                                                       tgatgaCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTg                                                              <  2:721331/150‑1 (MQ=255)
                                                                                       tgatgaCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTg                                                              <  2:31222/150‑1 (MQ=255)
                                                                                         atgaCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGAc                                                            >  1:1430769/1‑150 (MQ=255)
                                                                                          tgaCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACg                                                           <  1:511504/150‑1 (MQ=255)
                                                                                             ccATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGtct                                                        >  1:502411/1‑150 (MQ=255)
                                                                                              catcatAGTTATGGCGATGGTGTTGCTATTAGGTTATCTGGGTTCCTCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTGTATTATGACGGCGCTTCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGtctc                                                       <  2:642976/150‑1 (MQ=255)
                                                                                                 catcGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGt                                                    <  1:63902/150‑1 (MQ=255)
                                                                                                  atcGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGtt                                                   <  1:834923/150‑1 (MQ=255)
                                                                                                         tGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGc                                            <  1:1022042/150‑1 (MQ=255)
                                                                                                                ggTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCa                                     <  1:727509/150‑1 (MQ=255)
                                                                                                                  tgtTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCAtt                                   <  2:1431466/150‑1 (MQ=255)
                                                                                                                   gtTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATta                                  >  1:387930/1‑150 (MQ=255)
                                                                                                                   gtTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATCACGCCGCTGCAAGCGCTTCAACCGTTCTTGCCACGCGTAGGGATTGACGTCTCAGTTGACCCGCTGGCGCAAta                                  >  2:1532133/1‑150 (MQ=255)
                                                                                                                   gtTGCTATTAGCTTACCTGGGTTCCGCCCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATta                                  <  1:1479921/150‑1 (MQ=255)
                                                                                                                   gtTGCTATTAGCTTACCTGGGTTCCGCCCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATta                                  <  1:1485125/150‑1 (MQ=255)
                                                                                                                        tATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTg                             >  1:200851/1‑150 (MQ=255)
                                                                                                                         aTTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGa                            >  1:1126796/1‑150 (MQ=255)
                                                                                                                           tAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATc                          >  1:419860/1‑150 (MQ=255)
                                                                                                                             gCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCg                        <  2:1429000/150‑1 (MQ=255)
                                                                                                                             gCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGATGTTGCAACGCGTAGGGATTGACGTCTCAGTTGCCCAGCTGGCGCATTATATTGATCCg                        >  1:613844/1‑150 (MQ=255)
                                                                                                                                tATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAAc                     >  1:1005145/1‑150 (MQ=255)
                                                                                                                                   cTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAAcgcg                  >  1:1177507/1‑150 (MQ=255)
                                                                                                                                    tGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACgcgg                 >  1:821112/1‑150 (MQ=255)
                                                                                                                                     gggTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACgcggc                >  2:659626/1‑150 (MQ=255)
                                                                                                                                     gggTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACgcggc                >  2:943406/1‑150 (MQ=255)
                                                                                                                                     gggTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACgcggc                <  1:1223470/150‑1 (MQ=255)
                                                                                                                                      ggTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACgcggcg               <  2:419860/150‑1 (MQ=255)
                                                                                                                                         tCCGCGTCAACTCGTTGAGCGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATg            >  1:909316/1‑150 (MQ=255)
                                                                                                                                          ccGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGa           >  1:70221/1‑150 (MQ=255)
                                                                                                                                           cgcgCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGAc          >  2:927659/1‑150 (MQ=255)
                                                                                                                                           ccgcgtCAACTCGTTGAGCGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTAATCCGAACGCGGCGATGa           >  2:484062/6‑150 (MQ=255)
                                                                                                                                             cgcTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGt        <  1:659626/150‑1 (MQ=255)
                                                                                                                                              gcTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGtt       <  1:582440/150‑1 (MQ=255)
                                                                                                                                                tCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGGTGc     <  2:665416/150‑1 (MQ=255)
                                                                                                                                                 cAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCt    <  2:1209164/150‑1 (MQ=255)
                                                                                                                                                 cAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCt    >  2:198588/1‑150 (MQ=255)
                                                                                                                                                   aCGAGTTGACGCGGACGTAACCGAAATAACGAAAAGATATTATGACCCCGATGCAAGCGATTGACCCGTTTTTGCAACGCGTAGGGATTGACGGAGAAGTTGACCAGATGGCGCATTTTATTGAACCGAACGCGGCAATGCCGTTGCtaa  >  2:98173/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCGCATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCTCA  >  NC_000913/1674028‑1674324

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