Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,405,941 C→G 5.9% G234A (GGC→GCC)  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,9410CG 5.9% 339.7 / 9.4 102G234A (GGC→GCC) lrhAtranscriptional repressor of flagellar, motility and chemotaxis genes
Reads supporting (aligned to +/- strand):  ref base C (52/44);  new base G (4/2);  total (56/46)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.63e-01

Only 100 of 105 total aligned reads displayed.

GCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGG‑GATCGGCTCG  >  NC_000913/2405794‑2406090
                                                                                                                                                   |                                                                                                                                                      
gCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCa                                                                                                                                                      >  2:729282/1‑150 (MQ=255)
gCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCa                                                                                                                                                      <  2:1654887/150‑1 (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: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                                                                                                                                          <  1:193727/146‑1 (MQ=255)
                gggCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTc                                                                                                                                          >  2:193727/1‑146 (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:878079/150‑1 (MQ=255)
                   ccAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCg                                                                                                                                   >  2:485128/1‑150 (MQ=255)
                   ccAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCg                                                                                                                                   >  1:1234823/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)
                                                              ggTATCAGGTAGCGGCGGTAAGCCAGCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCa                                                                                        <  2:738198/150‑1 (MQ=255)
                                                                 aTCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATgg                                                                                     <  2:1658505/150‑1 (MQ=255)
                                                                 aTCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATgg                                                                                     <  1:10074/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                                                                   >  2:602960/1‑150 (MQ=255)
                                                                                   gCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAg                                                                   <  1:773753/150‑1 (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)
                                                                                                                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                  >  1:1302476/1‑150 (MQ=255)
                                                                                                                                    ggCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACa                  <  2:1529244/150‑1 (MQ=255)
                                                                                                                                     gCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACaa                 <  2:221152/150‑1 (MQ=255)
                                                                                                                                       cTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGACGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACCagc               >  2:1295864/1‑149 (MQ=255)
                                                                                                                                             gTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGG‑GATc        <  1:485128/150‑1 (MQ=255)
                                                                                                                                                aCGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCCGCAATACAAGAGG‑GATCGGc     >  2:389772/1‑150 (MQ=255)
                                                                                                                                                aCGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGG‑GATCGGc     <  2:904650/150‑1 (MQ=255)
                                                                                                                                                 cGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGG‑GATCGGCt    >  1:1019386/1‑150 (MQ=255)
                                                                                                                                                   ccAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGG‑GATCGGCTCg  <  1:30452/150‑1 (MQ=255)
                                                                                                                                                   ccAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCCCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGG‑GATCGGCTCg  >  2:663743/1‑150 (MQ=255)
                                                                                                                                                   ccAAGACCGGCTATCACTGCCGCACGAGCGGCCGGAAGCGTCGAGGCGCCATAAGCCAGTCGCCATGGAATCTCTGCTATATTAAGCGTCGCCAGCACCAAAACGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGGAGATCggaag   >  1:249494/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
GCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGG‑GATCGGCTCG  >  NC_000913/2405794‑2406090

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