Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 429,862 G→A 5.2% E120K (GAA→AAA)  secF → SecYEG protein translocase auxillary subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913429,8620GA 5.2% 539.9 / 11.7 153E120K (GAA→AAA) secFSecYEG protein translocase auxillary subunit
Reads supporting (aligned to +/- strand):  ref base G (60/85);  new base A (4/4);  total (64/89)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 153 total aligned reads displayed.

AATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTCGTGTACGTAGGTTTCCGCTTTGAGTGGCGAC  >  NC_000913/429714‑430009
                                                                                                                                                    |                                                                                                                                                   
aaTGCGTGATGCATTGCAAACAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGa                                                                                                                                                    <  2:1471269/150‑1 (MQ=255)
 aTGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGaa                                                                                                                                                   >  2:164141/1‑150 (MQ=255)
 aTGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGaa                                                                                                                                                   <  2:450742/150‑1 (MQ=255)
 aTGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGaa                                                                                                                                                   >  2:568982/1‑150 (MQ=255)
  tGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAAt                                                                                                                                                  <  1:1504817/150‑1 (MQ=255)
        aTGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATccacca                                                                                                                                            >  1:754014/1‑150 (MQ=255)
          gCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCTTGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAAt                                                                                                                                          <  1:1557061/150‑1 (MQ=255)
          gCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAAt                                                                                                                                          <  2:768262/150‑1 (MQ=255)
           cATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATc                                                                                                                                         <  1:2036958/150‑1 (MQ=255)
                 aaaaaGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATg                                                                                                                                   >  1:1404703/1‑150 (MQ=255)
                 aaaaaGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATg                                                                                                                                   >  1:1704291/1‑150 (MQ=255)
                        cGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGt                                                                                                                            <  2:1961611/150‑1 (MQ=255)
                        cGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGt                                                                                                                            <  2:1033410/150‑1 (MQ=255)
                        cGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGt                                                                                                                            <  2:2577663/150‑1 (MQ=255)
                        cGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGt                                                                                                                            <  2:614324/150‑1 (MQ=255)
                             ttGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGc                                                                                                                       <  1:853094/150‑1 (MQ=255)
                                aaGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTa                                                                                                                    <  2:304214/150‑1 (MQ=255)
                                 agagCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTAt                                                                                                                   >  2:337887/1‑150 (MQ=255)
                                 agagCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTAt                                                                                                                   >  2:337905/1‑150 (MQ=255)
                                  gagCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTAtt                                                                                                                  <  1:1736478/150‑1 (MQ=255)
                                  gagCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTAtt                                                                                                                  <  1:1736680/150‑1 (MQ=255)
                                  gagCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTAtt                                                                                                                  >  2:1186036/1‑150 (MQ=255)
                                  gagCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTAtt                                                                                                                  >  1:1900071/1‑150 (MQ=255)
                                  gagCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTAtt                                                                                                                  <  1:2697118/150‑1 (MQ=255)
                                      cGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGt                                                                                                              >  1:2582114/1‑150 (MQ=255)
                                      cGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGt                                                                                                              >  1:1117714/1‑150 (MQ=255)
                                      cGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGt                                                                                                              <  1:2680564/150‑1 (MQ=255)
                                      cGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGt                                                                                                              <  1:1463527/150‑1 (MQ=255)
                                          gctgcAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTtcgt                                                                                                          <  2:2446854/150‑1 (MQ=255)
                                          gctgcAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTtcgt                                                                                                          <  2:729297/150‑1 (MQ=255)
                                          gctgcAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTtcgt                                                                                                          <  2:411466/150‑1 (MQ=255)
                                          gctgcAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTtcgt                                                                                                          <  2:411436/150‑1 (MQ=255)
                                                  aCTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGa                                                                                                  <  2:1002091/150‑1 (MQ=255)
                                                   cTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAg                                                                                                 <  1:1938657/150‑1 (MQ=255)
                                                   cTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAg                                                                                                 <  1:1485515/150‑1 (MQ=255)
                                                          agcagcCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTgggg                                                                                          >  2:1233847/1‑150 (MQ=255)
                                                          agcagcCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACAAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTgggg                                                                                          >  1:2361084/1‑150 (MQ=255)
                                                          agcagcCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACAAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTgggg                                                                                          >  1:1132709/1‑150 (MQ=255)
                                                          agcagcCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACAAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTgggg                                                                                          >  1:1160422/1‑150 (MQ=255)
                                                          agcagcCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACAAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTgggg                                                                                          >  1:2748352/1‑150 (MQ=255)
                                                              gcCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAg                                                                                      <  1:1011920/150‑1 (MQ=255)
                                                               ccATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGa                                                                                     >  1:459379/1‑150 (MQ=255)
                                                                cATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGAc                                                                                    <  2:2223433/150‑1 (MQ=255)
                                                                cATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGAc                                                                                    <  2:2218463/150‑1 (MQ=255)
                                                                   gACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCtt                                                                                 >  2:1123615/1‑150 (MQ=255)
                                                                   gACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCtt                                                                                 >  2:1123616/1‑150 (MQ=255)
                                                                    acctccgggTCCGGAGGCCGCCTGCTGCCGGCGAACCCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGCAGGTGATGAACGAAGCCACCAATCAGAATGCAGCAGTGCAGCGTATTGAGTTCGTCGGTCCGCGCGTGGGGGCAGACCTTg                                                                                <  2:2460950/143‑1 (MQ=255)
                                                                    aCAGCAGGGGCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTg                                                                                <  2:1092996/150‑1 (MQ=255)
                                                                     catcatGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTgc                                                                               <  2:1608684/150‑1 (MQ=255)
                                                                     catcatGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTgc                                                                               <  1:805657/150‑1 (MQ=255)
                                                                     catcatGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTgc                                                                               <  1:1469186/150‑1 (MQ=255)
                                                                          tGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCaaa                                                                          >  2:1554896/1‑150 (MQ=255)
                                                                          tGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCaaa                                                                          >  1:1862822/1‑150 (MQ=255)
                                                                           ggTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAAc                                                                         <  2:754014/150‑1 (MQ=255)
                                                                             tCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCg                                                                       >  1:1000355/1‑150 (MQ=255)
                                                                                   tGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGa                                                                 <  2:1677072/150‑1 (MQ=255)
                                                                                    gccgccTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGAt                                                                <  2:2186781/150‑1 (MQ=255)
                                                                                    gccgccTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGAt                                                                <  1:203327/150‑1 (MQ=255)
                                                                                      cgccTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATgg                                                              <  2:1878520/150‑1 (MQ=255)
                                                                                              gAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATg                                                      <  1:220699/150‑1 (MQ=255)
                                                                                                  gCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAg                                                  >  1:332964/1‑150 (MQ=255)
                                                                                                    gAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAgcg                                                >  2:1231873/1‑150 (MQ=255)
                                                                                                    gAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAgcg                                                <  1:954053/150‑1 (MQ=255)
                                                                                                    gAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAgcg                                                >  1:114926/1‑150 (MQ=255)
                                                                                                    gAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGGTGATGGCAgcg                                                >  1:2607732/1‑150 (MQ=255)
                                                                                                      aaCCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCgct                                              >  2:1370793/1‑150 (MQ=255)
                                                                                                      aaCCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCgct                                              >  2:438402/1‑150 (MQ=255)
                                                                                                         cggcggTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCgctgct                                           <  2:2613974/150‑1 (MQ=255)
                                                                                                          ggcggTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGctgctg                                          >  1:2520302/1‑150 (MQ=255)
                                                                                                                   gtgtTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCACACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTc                                 <  1:214679/150‑1 (MQ=255)
                                                                                                                   gtgtTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTc                                 <  1:2707423/150‑1 (MQ=255)
                                                                                                                   gtgtTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTc                                 <  1:2636702/150‑1 (MQ=255)
                                                                                                                   gtgtTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTc                                 <  2:1611934/150‑1 (MQ=255)
                                                                                                                   gtgtTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTc                                 <  1:652061/150‑1 (MQ=255)
                                                                                                                   gtgtTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTc                                 <  1:1226081/150‑1 (MQ=255)
                                                                                                                   gtgtTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTc                                 <  2:101039/150‑1 (MQ=255)
                                                                                                                               cAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTCGTGTACGTAGGt                     >  1:927736/1‑150 (MQ=255)
                                                                                                                               cAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTCGTGTACGTAGGt                     >  1:1794694/1‑150 (MQ=255)
                                                                                                                               cAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTCGTGTACGTAGGt                     >  1:378696/1‑150 (MQ=255)
                                                                                                                               cAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTCGTGTACGTAGGt                     >  1:295060/1‑150 (MQ=255)
                                                                                                                                aGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTCGTGTACGTAGGtt                    >  1:1525475/1‑150 (MQ=255)
                                                                                                                                aGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTCGTGTACGTAGGtt                    >  1:1844206/1‑150 (MQ=255)
                                                                                                                                aGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTCGTGTACGTAGGtt                    >  1:1862748/1‑150 (MQ=255)
                                                                                                                                   ttCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTCGTGTACGTAGGTTTcc                 <  1:1233847/150‑1 (MQ=255)
                                                                                                                                   ttCTGAAGGTGATTAACAAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTCGTGTACGTAGGTTTcc                 <  2:2361084/150‑1 (MQ=255)
                                                                                                                                   ttCTGAAGGTGATTAACAAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTCGTGTACGTAGGTTTcc                 <  2:1160422/150‑1 (MQ=255)
                                                                                                                                   ttCTGAAGGTGATTAACAAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTCGTGTACGTAGGTTTcc                 <  2:2748352/150‑1 (MQ=255)
                                                                                                                                     cTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTCGTGTACGTAGGTTTCCGc               >  1:2367556/1‑150 (MQ=255)
                                                                                                                                       gAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTCGTGTACGTAGGTTTCCGCtt             >  1:1429358/1‑150 (MQ=255)
                                                                                                                                          ggTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTCGTGTACGTAGGTTTCCGCTTTGa          <  2:1404703/150‑1 (MQ=255)
                                                                                                                                          ggTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTCGTGTACGTAGGTTTCCGCTTTGa          >  1:938394/1‑150 (MQ=255)
                                                                                                                                           gTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTCGTGTACGTAGGTTTCCGCTTTGAg         >  1:387056/1‑150 (MQ=255)
                                                                                                                                            tGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTCGTGTACGTAGGTTTCCGCTTTGAGt        <  2:864075/150‑1 (MQ=255)
                                                                                                                                             gATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTCGTGTACGTAGGTTTCCGCTTTGAGTg       >  1:714302/1‑150 (MQ=255)
                                                                                                                                              aTTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTCGTGTACGTAGGTTTCCGCTTTGAGTgg      <  1:346233/150‑1 (MQ=255)
                                                                                                                                              aTTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTCGTGTACGTAGGTTTCCGCTTTGAGTgg      <  2:1339568/150‑1 (MQ=255)
                                                                                                                                              aTTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTCGTGTACGTAGGTTTCCGCTTTGAGTgg      <  2:828654/150‑1 (MQ=255)
                                                                                                                                              aTTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTCGTGTACGTAGGTTTCCGCTTTGAGTgg      >  2:1550793/1‑150 (MQ=255)
                                                                                                                                                 aaCGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTCGTGTACGTAGGTTTCCGCTTTGAGTGGCGa   >  2:2522791/1‑150 (MQ=255)
                                                                                                                                                 aaCGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTAGTCTGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTCGTGTACGTAGGTTTCCGCTTTGAGTGGCGa   >  2:1450057/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTCGTGTACGTAGGTTTCCGCTTTGAGTGGCGAC  >  NC_000913/429714‑430009

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