Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,405,933 C→G 6.1% A237P (GCA→CCA)  lrhA ← transcriptional repressor of flagellar, motility and chemotaxis genes

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,405,9330CG 6.1% 323.2 / 9.3 98A237P (GCA→CCA) lrhAtranscriptional repressor of flagellar, motility and chemotaxis genes
Reads supporting (aligned to +/- strand):  ref base C (48/43);  new base G (4/2);  total (53/45)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

GCTCTCCATCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGA  >  NC_000913/2405784‑2406078
                                                                                                                                                     |                                                                                                                                                 
gCTCTCCATCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGc                                                                                                                                                   >  1:577897/1‑150 (MQ=255)
 ctctCCATCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGcc                                                                                                                                                  >  2:397831/1‑150 (MQ=255)
  tctcCATCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCg                                                                                                                                                 <  1:556010/150‑1 (MQ=255)
       aTCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACg                                                                                                                                            <  2:1349077/150‑1 (MQ=255)
          gCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCa                                                                                                                                         <  2:1654887/150‑1 (MQ=255)
          gCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCa                                                                                                                                         >  2:729282/1‑150 (MQ=255)
            cTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAg                                                                                                                                       <  2:704850/150‑1 (MQ=255)
              gATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGTTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGAc                                                                                                                                     <  1:1080027/150‑1 (MQ=255)
              gATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGAc                                                                                                                                     >  1:544730/1‑150 (MQ=255)
              gATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGAc                                                                                                                                     >  1:551355/1‑150 (MQ=255)
               aTAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGAcc                                                                                                                                    <  1:537353/150‑1 (MQ=255)
                   aTCACCTGGGCCAGCTCATTATTACTGGCCGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCt                                                                                                                                <  1:1312108/150‑1 (MQ=255)
                      aCCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGTCTTGGCGTGACGGCAAGGCCggttgag                                                                                                                             <  2:979135/150‑6 (MQ=255)
                      aCCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGTCTTGGCGTGACGGCAAGGCCggttgag                                                                                                                             <  2:494810/150‑6 (MQ=255)
                      aCCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTc                                                                                                                             >  1:1658505/1‑150 (MQ=255)
                      aCCTGGGCCAGCTCAATATTACTCGACGGAACATACCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCCCTCAGCACACGCAGGTCCGTGCGCATGATCTCAACCGTCCTTGCCGTCACGCCAAGACCGGCTTTc                                                                                                                             <  2:1046655/150‑1 (MQ=255)
                       ccTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCa                                                                                                                            >  2:1623483/1‑150 (MQ=255)
                       ccTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATATCAACCTGCCTTGCCGTCACGCCAAGACCGGCTTTCa                                                                                                                            <  2:132017/150‑1 (MQ=255)
                         tGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACt                                                                                                                          >  1:1190641/1‑150 (MQ=255)
                          gggCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTc                                                                                                                             >  2:193727/1‑146 (MQ=255)
                          gggCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTc                                                                                                                             <  1:193727/146‑1 (MQ=255)
                          gggCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTg                                                                                                                         <  2:482096/150‑1 (MQ=255)
                           ggCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTCTCACTGc                                                                                                                        <  2:1205913/150‑1 (MQ=255)
                            gCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGcc                                                                                                                       >  2:1501840/1‑150 (MQ=255)
                             ccAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCg                                                                                                                      >  1:1234823/1‑150 (MQ=255)
                             ccAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCg                                                                                                                      <  1:878079/150‑1 (MQ=255)
                             ccAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCg                                                                                                                      >  2:485128/1‑150 (MQ=255)
                               aGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCa                                                                                                                    >  2:30452/1‑150 (MQ=255)
                                   cATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGac                                                                                                                >  1:1438155/1‑149 (MQ=255)
                                     ttattaCTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACg                                                                                                              <  2:398704/150‑1 (MQ=255)
                                      tattaCTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCCGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACgg                                                                                                             <  2:1233328/150‑1 (MQ=255)
                                      tattaCTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACgg                                                                                                             >  2:1069968/1‑150 (MQ=255)
                                       attaCTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCTGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGc                                                                                                            <  2:1186947/150‑1 (MQ=255)
                                        ttaCTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGcc                                                                                                           >  1:904650/1‑150 (MQ=255)
                                        ttaCTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGcc                                                                                                           >  1:221152/1‑150 (MQ=255)
                                          aCTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCgg                                                                                                         >  1:1492857/1‑150 (MQ=255)
                                              gACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGc                                                                                                     >  2:766007/1‑150 (MQ=255)
                                                cGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGt                                                                                                   <  1:655665/150‑1 (MQ=255)
                                                cGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCGTTCACTGCCGCACGACCGGCCGGAcgcgg                                                                                                   >  1:129765/1‑149 (MQ=255)
                                                     cATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAgg                                                                                              <  1:932192/150‑1 (MQ=255)
                                                      aTAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGc                                                                                             >  1:1666353/1‑150 (MQ=255)
                                                       tAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCg                                                                                            >  1:1529244/1‑150 (MQ=255)
                                                         gCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGAc                                                                                          >  1:1192085/1‑150 (MQ=255)
                                                          cacaGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGCGGCGTCa                                                                                         >  2:1649310/1‑150 (MQ=255)
                                                          cacaGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACa                                                                                         <  2:1716222/150‑1 (MQ=255)
                                                           acaGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCTGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTAGAGGCGACAt                                                                                        >  2:194230/1‑150 (MQ=255)
                                                              gCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGACCGGCCGGAAGCGTCGAGGCGCCATCAg                                                                                     >  1:1409197/1‑150 (MQ=255)
                                                              gCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAg                                                                                     >  1:8995/1‑150 (MQ=255)
                                                                  atatTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAg                                                                                 <  2:484994/150‑1 (MQ=255)
                                                                    atTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTc                                                                               >  2:1428518/1‑150 (MQ=255)
                                                                    atTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTc                                                                               <  1:428240/150‑1 (MQ=255)
                                                                        ggTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCa                                                                           >  1:622693/1‑150 (MQ=255)
                                                                        ggTATCAGGTAGCGGCGGTAAGCCAGCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCa                                                                           <  2:738198/150‑1 (MQ=255)
                                                                           aTCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATgg                                                                        <  1:10074/150‑1 (MQ=255)
                                                                           aTCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATgg                                                                        <  2:1658505/150‑1 (MQ=255)
                                                                                   gcggcgGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCATCGAGGCGACATAAGCCAGTCGCCATGGAATATCTg                                                                <  2:1190641/150‑1 (MQ=255)
                                                                                    cggaggTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGACGTCACGCCAAGACCGGCTTTCACTGCCGCACCAACGGCCGGAAGCGTCGAGGCGACATAAGCAAGTCGCCATGGAATATCTGc                                                               >  2:1487512/1‑150 (MQ=255)
                                                                                     ggcggTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCt                                                              <  2:577897/150‑1 (MQ=255)
                                                                                       cggTAAGCGATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCATTTCCGTCACGCCAAGAGCGGCTTTCACTGCCGCTCGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCttt                                                            <  2:1148926/150‑1 (MQ=255)
                                                                                           aaGCCAGCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTttatta                                                        <  2:324356/150‑1 (MQ=255)
                                                                                            aGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTaa                                                       <  1:1623483/150‑1 (MQ=255)
                                                                                            aGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTaa                                                       <  2:8995/150‑1 (MQ=255)
                                                                                             gCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAg                                                      <  1:773753/150‑1 (MQ=255)
                                                                                             gCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAg                                                      >  2:602960/1‑150 (MQ=255)
                                                                                              ccATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAga                                                                                    <  1:313937/119‑2 (MQ=255)
                                                                                              ccATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGc                                                                                    >  2:313937/1‑119 (MQ=255)
                                                                                                 tCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTc                                                  <  1:1069968/150‑1 (MQ=255)
                                                                                                    aCTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGACCGGCCGGCAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGcc                                               >  2:1480577/1‑150 (MQ=255)
                                                                                                    aCTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGCCATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGcc                                               >  1:511577/1‑150 (MQ=255)
                                                                                                     cTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATTAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGccc                                              >  1:624667/1‑149 (MQ=255)
                                                                                                     cTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCa                                              >  1:291648/1‑150 (MQ=255)
                                                                                                       ccACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGCAGCGTCGAGGCGACATCAGCCAGTCGCCATGGAATATCTGCGTTTTTATGCGTCGCCCGc                                            >  1:148102/1‑150 (MQ=255)
                                                                                                        cACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCa                                           >  1:777810/1‑150 (MQ=255)
                                                                                                          cTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCAcc                                         >  2:1441882/1‑150 (MQ=255)
                                                                                                           tCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCa                                        >  1:851008/1‑150 (MQ=255)
                                                                                                           gcAGCACTCGCAGGGCCGGGCTCATCATCTCAACCGGCCGTGCCGGCACGCCAAGTCCGGCTGTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGGCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGGACCa                                        <  1:946944/149‑1 (MQ=255)
                                                                                                            cAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCat                                       >  2:1504658/1‑150 (MQ=255)
                                                                                                               cacaCGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCAATGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATc                                    >  2:259473/1‑150 (MQ=255)
                                                                                                                 cacGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATcgc                                  <  1:117440/150‑1 (MQ=255)
                                                                                                                  acGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATcgcg                                 <  2:1234823/150‑1 (MQ=255)
                                                                                                                   cGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCCCGACCGGCCGGAAGCGTCGAGGCGACATAAGCCCGTCGCCATGGACTATCTGCTTTATTAAGCGTCGCCCGCCCCATATCGCGa                                >  1:218063/1‑150 (MQ=255)
                                                                                                                    gCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCTGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGaa                               <  1:397831/150‑1 (MQ=255)
                                                                                                                       ggTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACg                            <  2:1114774/150‑1 (MQ=255)
                                                                                                                       ggTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACg                            >  2:673693/1‑150 (MQ=255)
                                                                                                                         tCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCCCGAACGGCCGGAAGCGTCGAGGCGCCATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCCCCCTATCGCGAAACggg                          >  1:1618810/1‑150 (MQ=255)
                                                                                                                          ccGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGc                         <  2:1325708/150‑1 (MQ=255)
                                                                                                                              gCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCAAATCGCGAAACGGGCTAgg                     >  1:1186005/1‑150 (MQ=255)
                                                                                                                                 catcatCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCCCCATATCGCGAAACGGGCTAGGctc                  >  1:5440/1‑150 (MQ=255)
                                                                                                                                     atcTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATcc              <  1:1148294/150‑1 (MQ=255)
                                                                                                                                     atcTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATcc              >  2:294427/1‑150 (MQ=255)
                                                                                                                                           gccGGTCTTGGCGTGACGGCAAGGCCGGCTTTCACTGCCGCACGAACGGCCGGACGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTACGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCaaa        >  1:268054/2‑149 (MQ=255)
                                                                                                                                           aCCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAAt        <  1:324130/150‑1 (MQ=255)
                                                                                                                                            ccGGTCTTGGCGTGACGGCAAGGCCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATa       >  1:1587359/1‑150 (MQ=255)
                                                                                                                                            ccGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAATCGGGCTAGGATCATCCAGCAATa       >  1:1590635/1‑150 (MQ=255)
                                                                                                                                             cGGTCTTGGCGTGACGGCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCCGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATAc      >  1:742452/1‑150 (MQ=255)
                                                                                                                                              ggtcTTGGCGTGACGGCAAGGCCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACa     >  2:2610/4‑150 (MQ=255)
                                                                                                                                              ggCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACa     <  2:1529244/150‑1 (MQ=255)
                                                                                                                                              ggCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACa     >  1:1302476/1‑150 (MQ=255)
                                                                                                                                               gCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACaa    <  2:221152/150‑1 (MQ=255)
                                                                                                                                                 cTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGACGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACCagc  >  2:1295864/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GCTCTCCATCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGA  >  NC_000913/2405784‑2406078

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