Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,976,691 T→C 13.3% E150G (GAG→GGG)  motA ← proton conductor component of flagella motor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,976,6910TC13.3% 170.9 / 3.0 75E150G (GAG→GGG) motAproton conductor component of flagella motor
Reads supporting (aligned to +/- strand):  ref base T (26/36);  new base C (4/7);  total (30/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGATGCGTGGGTAGCTGGCGAAGATCTCGCTCTCACGGGGATTTTCAATATCACGTTCCAGCGAAAACATCCCCATCTGCCGCGATTTCGCCATCAACCGATAAAGCAGAGCCAGCAGATCCATATACATTGCT  >  NC_000913/1976453‑1976924
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                         
cTGACTAAAAACGGGCCCTAATGGGGAAATAAACCCGTAAGCCCACAAAATGCCGAAGAAAGCCCCCACCCACGCATGTGCGATAAGCGCCCCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGCCCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAAtct                                                                                                                                                                                                                                           <  2:183351/239‑1 (MQ=255)
cTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGGGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAAtct                                                                                                                                                                                                                                           <  2:173457/239‑1 (MQ=255)
     gTAAAACAGTTGGTAATTGGGAAATAAATCCGTAAGCAAATAAAAGGCCGAAGAAAGTCCCCACCATCGCAAGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCGAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTca                                                                                                                                                                                                                                      <  1:822005/239‑1 (MQ=255)
          aCAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGCCCCCACCATCGCATGGGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTcatccat                                                                                                                                                                                                                                 <  1:270800/239‑1 (MQ=255)
              tCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATCCCAAACGCCGGACGTGAGTCCCCGACCACGGCCAGACCGTTCGCCGGGACTTCTGCCCCCCTTTCGGGCGGCCCAATCCCTTctcccccccaa                                                                                                                                                                                                                             >  2:1064254/1‑229 (MQ=255)
              tCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCTTCaga                                                                                                                                                                                                                             >  1:791784/1‑239 (MQ=255)
              tCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCCCTTCATCCATCata                                                                                                                                                                                                                             >  1:869467/1‑237 (MQ=255)
              tCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGGGAGTCCCCGACCAGCGCCACACCGTTCGCCGGGCCTCCTGCCTCGCTTTCGTCCGCCCCCACCCCTTcatcccccaca                                                                                                                                                                                                                             >  2:162886/1‑233 (MQ=255)
               cGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTCCAGCCGCAACAATACCAAACGCCGGAAGTGGGGCCCCGACCAGCCCCGACCTGTTCCCCGGGACTTCTCCCCCCCTTTCCTGTGTCCCAAACCCTTctcccataaaca                                                                                                                                                                                                                            >  2:256974/1‑234 (MQ=255)
                           aataataCCGTAAGCCTAAAAAAGGCCGAGGAAAGTCCCCCCCATCGCAAGTGGCATAAGCGGACCCAGCCCGGCGGCAGGACGATCGGCTGAACCTAACGCGGGAACGCCCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGa                                                                                                                                                                                                                <  2:183392/239‑1 (MQ=255)
                           aataaatCCGTAAGCCAATAAAAGGGCGAGGAAAGGCCCCCCCAGCGTCTGGGCGATAAGCGCACCCAGGCCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGa                                                                                                                                                                                                                <  2:133118/239‑1 (MQ=255)
                            ataaatCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGACTCGCTTTCGTGCGGCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGaa                                                                                                                                                                                                               >  1:1050412/1‑239 (MQ=255)
                             taaatCCGTAAGCCAAAAAAATGCCGAAGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAg                                                                                                                                                                                                              <  1:656530/239‑1 (MQ=255)
                             taaaccctgaaGCCAATAAAATGCCGAGGAAAGTCCCCCCCATCGCATGGGCGAAAAGCCCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGCCCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAg                                                                                                                                                                                                              <  2:791784/230‑1 (MQ=255)
                             taaacCCGTAAGCAAAAAAAAATCGGAGGAAAGTCCCCCCCATCGCAGGTGGGATAAGCGCACCCAGCTCGGGCGCAGGACGACCGGCAGAACCTACCGCGTGAACGCCCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAg                                                                                                                                                                                                              <  2:793684/239‑1 (MQ=255)
                                atCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACCGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCCCAATCTCTTCATCCTTCAGCGCTTCGATTTCGAAGgtg                                                                                                                                                                                                           >  1:1001015/1‑239 (MQ=255)
                                acccGTAAGCCAAAAAAAGGCCGAGGAAACTCCCCCCCATCGCATGGGCGATAAGGGCCCCCAGCTCGGCGGCAGGACGAACGGCTGAACCTAACGCGTGAACGCCCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGgtg                                                                                                                                                                                                           <  2:1050412/237‑1 (MQ=255)
                                          cAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGAc                                                                                                                                                                                                 <  1:149888/239‑1 (MQ=255)
                                          cAAAAAAATGACGAGGAAAGTCCCCCCCCGCGCAGGGGCGAAAAGCGCCCCCACCCCGGCGGCAGGACGATCGGCGGAACCCACCGCGTGAACGACCCCCCGTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGGTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGAc                                                                                                                                                                                                 <  2:565518/239‑1 (MQ=255)
                                          cAAAAAAAAGCCGGGGAAAGTCCCCACCATCTGATGGTGCAAAAGCGCCCCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGGGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGCCTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGAc                                                                                                                                                                                                 <  2:726264/239‑1 (MQ=255)
                                            aTAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGGTTTCATTTGAcc                                                                                                                                                                                                >  1:576600/1‑238 (MQ=255)
                                                agcgCCGGGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCCCGGCGGCAGGACGACCGGATGAACCTAGCGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCGGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGa                                                                                                                                                                                           <  1:1064254/236‑1 (MQ=255)
                                                       aGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGCCTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCCCTGATAATCAg                                                                                                                                                                                    >  1:951344/1‑239 (MQ=255)
                                                        ggAAAGTCCCCCGCATCGCAAGTGCGAAAAGCGGACCCAGCTGGGAGGCAGGACGATCGGCTGAACCTAGCGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAg                                                                                                                                                                                    <  1:115762/238‑1 (MQ=255)
                                                         gaagggccccACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAgg                                                                                                                                                                                   <  1:917429/232‑1 (MQ=255)
                                                                   accaTCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCACTCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCCCCGATAATc                                                                                                                                                                        >  1:758915/1‑239 (MQ=255)
                                                                   accaTCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGGGTTCATGTGACCGCTGCTAATCAGGCGCAGATAATc                                                                                                                                                                        >  1:45936/1‑239 (MQ=255)
                                                                        cGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCTCTTTCGAAGGTGTTCATGTGACCGCTGGTAATCAGGCGCAGATAATcgacgc                                                                                                                                                                   >  1:1049661/1‑238 (MQ=255)
                                                                        cGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGCTTTCGAAGGGGTTCATGGGCCCGCTGATAACCAGGCGCAGATAATcgccga                                                                                                                                                                   >  2:444792/1‑239 (MQ=255)
                                                                          gcagtggCGATAAGCGCCCCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGGGTGAACGACCCCCATTACAGCCGCAACAATCCCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGAt                                                                                                                                                                  <  1:444792/233‑1 (MQ=255)
                                                                          ctcttgtccgaaaGCCCCCCCAGCTCGGCGCCAGGCCGATCGGTGTAACCCAACCCGTGACCGACCCCGATTACAGCCGCAACAATCCCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGAt                                                                                                                                                                  <  2:134289/228‑1 (MQ=255)
                                                                                gCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCCACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGCTTTCGAAGGGGTTCATGTGACCGCTGATAATAAAGCGCAGATAATCGCGGAAAAAAta                                                                                                                                                            >  2:99381/1‑238 (MQ=255)
                                                                                ccatcaaGCCCCCCCAGCCCGGCGGCAGGACGACCGGCTGAACCTAACGCGTGACCGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATc                                                                                                                                                            <  2:919837/234‑1 (MQ=255)
                                                                                      aGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGCTAATCGACGATAAAATCAAGCAt                                                                                                                                                      >  1:607148/1‑239 (MQ=255)
                                                                                          gaCCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATGCCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACg                                                                                                                                                  <  1:403048/238‑1 (MQ=255)
                                                                                            cccACCCCGGGGGCCGGAGCACCGGCTGCACCTAAGGCGTGAACCACCCCCATTCCAGCCGCAAAAATGCCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCCCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCt                                                                                                                                                <  2:869467/239‑1 (MQ=255)
                                                                                                    ggcggcAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGa                                                                                                                                        <  1:162886/239‑1 (MQ=255)
                                                                                                      cggcAGGACGATCGGCTAAACCGAACGGGTGAACGACCCCCATTACATCCGCAACAATACCAACCGCCGGAAGTGGGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAg                                                                                                                                       <  1:1081966/238‑1 (MQ=255)
                                                                                                        gcAGGACGATCGGCTAAACCTAACGCTGGAACGACCCACATTACAGCCGCAACAATACCAAACGGCGGAAGTGAGTCCCCGGCCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGa                                                                                                                                     <  1:675317/238‑1 (MQ=255)
                                                                                                                gtCGGCTTACCCCACCCCGGGAACCACCCCCATTCCCGCCGAAACAATCCCAAACGCCGGAAGGGAGTCCCCGACCAGCGCCAGGCTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCGCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGATGCGTGGGt                                                                                                                            <  2:870335/238‑1 (MQ=255)
                                                                                                                                  gTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGGGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCAGAACGTTATCCGCAGGGATGCGTGGGTAGCTGGCGAAGCTATCg                                                                                                           >  1:991703/1‑238 (MQ=255)
                                                                                                                                   ggAACGCCCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGGGGCCCCGACCAGCGCCATGCTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGATGCGTGGGTAGCTGGCGAAGATCTCGct                                                                                                         <  2:758915/238‑1 (MQ=255)
                                                                                                                                    gAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCGTCCATCAGAGCTTCGATTTCGAAGGCGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGATGCGTGGGCAGCTGGCGAAGATCTCGct                                                                                                         >  1:160571/1‑238 (MQ=255)
                                                                                                                                     aaCGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGCGACCGCTGATAATCAGGCGCAGATAATCGACGTAAAAATCAAGCATGACGCTATCCGCGAGGATGCGTGGGTAGCTGGCGAAGATCTCGctc                                                                                                        >  1:1056808/1‑238 (MQ=255)
                                                                                                                                       ccccccccATTCCACCCCCTACTATTCCAGACGCCGCGTGTGGTGCCCCGACCCGTGCCAGACTGTTCGCCGTGCCTTCTTCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCAGCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGATGCGTGGGTAGCTGGCGAAGATCTCGCTCTCa                                                                                                     <  2:633304/236‑1 (MQ=255)
                                                                                                                                                 tACCGCCCCAACAATCCCAAACGCCGGAAGTGAGTCCCCGACCAGGGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGATGCGTGGGTAGCTGGCGAAGATCTCGCTCTCACGGGGAtttt                                                                                           <  2:991703/239‑1 (MQ=255)
                                                                                                                                                    aGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCCCGAGGATGCGTGGGTAGCGGGCGAAGATTCCCCTCCCCCGGGGATTTTCaa                                                                                        >  2:768743/1‑239 (MQ=255)
                                                                                                                                                      ccGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGATGCGTGGGTAGCTGGGGAAGATCTCGCCCTCACGGGGATTTTCaaa                                                                                       >  2:537379/1‑237 (MQ=255)
                                                                                                                                                        ccacaaaTCCCAAACGCCGGAAGGGGGCCCCCGACCAGCGCCAGGCTGTGCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCCCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGATGCGTGGGTAGCTGGCGAAGATCTCGCTCTCACGGGGATTTTCAATATc                                                                                    <  2:756528/234‑1 (MQ=255)
                                                                                                                                                         caacaaTACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGATGCGTGGGTAGCTGGCGAAGATCTCGCTCTCACGGGGATTTTCAATATc                                                                                    >  1:793033/1‑238 (MQ=255)
                                                                                                                                                          aacaaTCCCAAAGCCCGGAAGTGCGTCCCCGCCCAGGGCCAGACTGGTCGCCGGGACTTCTGCCCCGCTTTCGTGCGTCTCAATCTCTGCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGATGCGTGGGTAGCTGGCGAAGATCTCGCTCTCACGGGGATTTTCAATATc                                                                                    <  2:98529/237‑1 (MQ=255)
                                                                                                                                                                    aaCGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGATGCGTGGGTAGCTGGCGAAGATCTCGCTCTCACGGGGATTTTCAATATCACGTTCCAGc                                                                          >  1:98529/1‑237 (MQ=255)
                                                                                                                                                                    aaCGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGATGCGTGGGTAGCTGGCGAAGATCTCGCTCTCACGGGGATTTTCAATATCACGTTCCAGc                                                                          >  1:756528/1‑237 (MQ=255)
                                                                                                                                                                    aaCGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGATCCGTGGGTAGCTGGCGAAGATCTCGCTCTCACGGGGATTTTCAATATCACGTTCCAGCGa                                                                        >  1:1027016/1‑239 (MQ=255)
                                                                                                                                                                       gCCGGAAGTGAGTCCCCGACCGGCGCCAGACTGTTCGCCGGGACTTCGGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCAGCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGATGCGTGGGTAGCTGGCGAAGATCTCGCTCTCACGGGGATTTTCAATATCACGTTCCAGCGaaaa                                                                     <  1:424272/239‑1 (MQ=255)
                                                                                                                                                                           gAAGTGGGTCCCCGGCCAGCGCCAGCCTTTTCGCCGGGACTTCTGCCCCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGATGCGTGGGTAGCTGGCGAAGATCTCGCTCTCACGGGGATTTTCAATATCACGTTCCAGCGAAAACATc                                                                 <  2:953905/239‑1 (MQ=255)
                                                                                                                                                                                             gcgcCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGATGCGTGGGTAGCTGGCGAAGATCTCGCCCCACCGGGGATTTTCATTATCACGTTCCAGCGAAAAAATCCCCATCTACCGCGATTTc                                               >  2:474789/1‑239 (MQ=255)
                                                                                                                                                                                              cgcCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGAGGCGTGGGTAGCTGGGGAAGATCTCGCTCTAACGGGGATTTTCAAAATAAGGTTCCAGCGAAAAAATACCAACCTGCCGCGATTTc                                               >  2:497836/1‑238 (MQ=255)
                                                                                                                                                                                                    gcctttCGCCGGTTCTTTCTCCTCCCTTTCGTGCGTCTCAATCTCTTCATCCATCAGGGTTTCGAGTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAAGCAAGCATGACGCTATCCGCGAGGATGCGTGGGTAGCTGGCGAAGATTTCGCTCTCACGGGGATTTTCAATATCACGTTCCAGCGAAAACATCCCCATCTGCCGCGATTTCGCCATCa                                        <  2:425305/235‑1 (MQ=255)
                                                                                                                                                                                                     cttttCGCCGGGACTTCTGCCTCGCTTTCGGGCGTCCCAATCCCTTCATCCCTCAGAGCTTGGATTTCGAAGGTGTTCATGGGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGATGCGTGGGTAGCTGGCGAAGATCTCGCTCTCACGGGGATTTTCAATATCACGTTCCAGCGAAAACATCCCCATCTGCCGCGATTTCGCCATCa                                        <  2:793033/235‑1 (MQ=255)
                                                                                                                                                                                                     cttttCGCCGGGACTTCTGCCTCGCTTGTGGGCCTCTCAATCCCTTAATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGGGCCCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGATGCGTGGGTAGCTGGCGAAGATCTCGCTCTCACGGGGATTTTCAATATCACGTTCCAGCGAAAACATCCCCATCTGCCGCGATTTCGCCATCa                                        <  1:474789/235‑1 (MQ=255)
                                                                                                                                                                                                        ttCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGATGCGTGGGTAGCTGGCGAAGATCTCGCTCTCCCGGGGATTTTCAATATCCCGTTCCAGCGAAAACATCCCCATCTTCCCCGTTTTTCCCACCAACCg                                    >  2:882177/1‑239 (MQ=255)
                                                                                                                                                                                                        ttCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGAGGCGTGGGTAGCTGGCGAAGATCTCGCTCTCACGGGGATTTTCAATATCACGTTCCAGCGAAAACATCCCCATCTGCCGCGATTTCGCCATCAACCg                                    >  1:90843/1‑239 (MQ=255)
                                                                                                                                                                                                          cGCCGGGCCTTTTGCTGTGCTTTCGTTGGTCTCAATCCCTTAATCCGACAGAGCTTAGTGTTCGAAGGTGGTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGATGCGTGGGTAGCTGGCGAAGATCTCGCTCTCACGGGGATTTTCAATATCACGTTCCAGCGAAAACATCCCCATCTGCCGCGATTTCGCCATCAACCGAt                                  <  2:160571/239‑1 (MQ=255)
                                                                                                                                                                                                                         cctgttttctggCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGGGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGATGCGTGGGTAGCTGGCGAAGATCTCGCTCTCACGGGGATTTTCAATATCACGTTCCAGCGAAAACATCCCCATCTGCCGCGATTTCGCCATCAACCGATAAAGCAGAGCcagca                   <  2:1049661/228‑1 (MQ=255)
                                                                                                                                                                                                                         cccggcttCGTGGGTCACAATCCCTTCATCCCACAGGGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGATGCGTGGGTAGCTGGCGAAGATCTCGCTCTCACGGGGATTTTCAATATCACGTTCCAGCGAAAACATCCCCATCTGCCGCGATTTCGCCATCAACCGATAAAGCAGAGCcagca                   <  1:918220/233‑1 (MQ=255)
                                                                                                                                                                                                                         cacGCTTTCGGGCGTCTCAATCGCTTCAACCAACAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGATGCGTGGGTAGCTGGCGAAGATCTCGCTCTCACGGGGATTTTCAATATCACGTTCCAGCGAAAACATCCCCATCTGCCGCGATTTCGCCATCAACCGATAAAGCAGAGCcagca                   <  1:882177/237‑1 (MQ=255)
                                                                                                                                                                                                                         cTCGGTTTCGTGCCGCTCAATCCCTTCATCCAACAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGATGCGTGGGTAGCTGGCGAAGATCTCGCTCTCACGGGGATTTTCAATATCACGTTCCAGCGAAAACATCCCCATCTGCCGCGATTTCGCCATCAACCGATAAAGCAGAGCcagca                   <  1:486136/239‑1 (MQ=255)
                                                                                                                                                                                                                         cTCGCTTTCGTGCGTCTCAATCTCATCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGATGCGTGGGTAGCTGGCGAAGATCTCGCTCTCACGGGGATTTTCAATATCACGTTCCAGCGAAAACATCCCCATCTGCCGCGATTTCGCCATCAACCGATAAAGCAGAGCcagca                   <  1:768743/239‑1 (MQ=255)
                                                                                                                                                                                                                         cTCGCTTTCGGTCGTCTCAATCTCTTCATCCATCAGGGCTTTGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGATGCGTGGGTAGCTGGCGAAGATCTCGCTCTCACGGGGATTTTCAATATCACGTTCCAGCGAAAACATCCCCATCTGCCGCGATTTCGCCATCAACCGATAAAGCAGAGCcagca                   <  2:1001015/239‑1 (MQ=255)
                                                                                                                                                                                                                         cTCGCATTCGGTCGTCTCAATCTCTTCATCCATCGATGCTTCGATGTCGAAGGTGTTCATGTGCGCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGATGCGTGGGTAGCTGGCGAAGATCTCGCTCTCACGGGGATTTTCAATATCACGTTCCAGCGAAAACATCCCCATCTGCCGCGATTTCGCCATCAACCGATAAAGCAGAGCcagca                   <  1:584983/239‑1 (MQ=255)
                                                                                                                                                                                                                          tctgtttCGGGCGTCCCAATCCCTTCATCCAACAGGGCTTCGATTTCGAAGGTGTTCATGGGCCCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGATGCGTGGGTAGCTGGCGAAGATCTCGCTCTCACGGGGATTTTCAATATCACGGTCCAGCGAAAACATCCCCATCTGCCGCGATTTCGCCATCAACCGATAAAGCAGAGCcagca                   <  1:19466/234‑1 (MQ=255)
                                                                                                                                                                                                                           cGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCTATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGATGCGTGGGTAGCTGGCGAAGATCTCGCTCTCACGGGGATTTTCAATATCACGTTCCAGCGAAAACATCCCCATCTGCCGCGATTTCGCCATCAACCGATAAAGCAGAGCCAGCAGa                 <  1:497836/239‑1 (MQ=255)
                                                                                                                                                                                                                                     cGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGATGCGTGGGTAGCTGGCGAAGATCTCGCTCTCACGGGGATTTTCAATATCACGTTCCAGCGAAAACATCCCCATCTGCCGCGATTTCGCCATCAACCGATAAAGCAGAGccacca                   >  1:700812/1‑227 (MQ=255)
                                                                                                                                                                                                                                      gtttccaTCTCTTCATCCATCAGAGCTTGGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGATGCGTGGGTAGCTGGCGAAGATCTCGCTCTCACGGGGATTTTCAATATCACGTTCCAGCGAAAACATCCCCATCTGCCGCGATTTCGCCATCAACCGATAAAGCAGAGCCAGCAGATCCATATACa       <  2:45936/232‑1 (MQ=255)
                                                                                                                                                                                                                                       tctcAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGATGCGTGGGTAGCTGGCGAAGATCTCGCTCTCACGGGGATTTTCAATATCACGTTCCAGCGAAAAAATCCCCACCTGCCGCGATTTCTCCATCAACCGATAAAGCAACGCCAGCCCATCCCTATACAt      >  1:481162/1‑238 (MQ=255)
                                                                                                                                                                                                                                          cAATCTCTTCATCCATACGAGCTTCGATTTCGAAGGTTTTCATGTGACCGCTGATAAACAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGATGCGTGGGTAGCTGGCGAAGATCTCGCTCTCACGGGGATTTTCAATATCACGTTCCAGCGAAAACATCCCCATCTGCCGCGATTTCGCCATCAACCGATAAAGCAGAGCCAGCAGATCCATATACATTGCt  <  1:468733/239‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                         
CTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGATGCGTGGGTAGCTGGCGAAGATCTCGCTCTCACGGGGATTTTCAATATCACGTTCCAGCGAAAACATCCCCATCTGCCGCGATTTCGCCATCAACCGATAAAGCAGAGCCAGCAGATCCATATACATTGCT  >  NC_000913/1976453‑1976924

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: