Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,971,929 A→C 13.4% S255A (TCA→GCA)  tar ← methyl‑accepting chemotaxis protein Tar

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,971,9290AC13.4% 221.6 / 30.2 90S255A (TCA→GCA) tarmethyl‑accepting chemotaxis protein Tar
Reads supporting (aligned to +/- strand):  ref base A (48/29);  new base C (6/6);  total (54/36)
Fisher's exact test for biased strand distribution p-value = 5.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.05e-01

GCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGC  >  NC_000913/1971781‑1972077
                                                                                                                                                    |                                                                                                                                                    
gcAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGaa                                                                                                                                                     <  1:220806/150‑1 (MQ=255)
     ttcttcCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGct                                                                                                                                                >  2:864878/1‑150 (MQ=255)
      tcttcCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTTCACATATGCAACGctc                                                                                                                                               <  1:140850/150‑1 (MQ=255)
       cttcCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGctct                                                                                                                                              <  1:606635/150‑1 (MQ=255)
         tcCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTTCACATATGCAACGctcttt                                                                                                                                            <  1:121983/150‑3 (MQ=255)
         tcCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTTCACATATGCAACGctcttt                                                                                                                                            <  2:1388057/150‑3 (MQ=255)
         tcCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTTCACATATGCAACGctcttt                                                                                                                                            <  2:1339698/150‑3 (MQ=255)
         tcCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTgc                                                                                                                                            <  1:1392823/150‑1 (MQ=255)
          ccAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTTCACATATGCAACGctctttg                                                                                                                                           <  1:1488115/150‑4 (MQ=255)
          ccAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTgcg                                                                                                                                           <  1:1300037/150‑1 (MQ=255)
              cgcgGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAg                                                                                                                                       <  1:1153989/150‑1 (MQ=255)
               gcgGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAgg                                                                                                                                      >  2:450962/1‑150 (MQ=255)
               gcgGATGCCTGCTGTGCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAgg                                                                                                                                      <  2:1182008/150‑1 (MQ=255)
                      ccTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCa                                                                                                                               >  2:1330145/1‑150 (MQ=255)
                        tgctgTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCAtt                                                                                                                             >  1:1208553/1‑150 (MQ=255)
                        tgctgTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCAtt                                                                                                                             >  1:401826/1‑150 (MQ=255)
                          ctgTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTc                                                                                                                           <  2:372877/150‑1 (MQ=255)
                          ctgTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTc                                                                                                                           <  1:570515/150‑1 (MQ=255)
                            gTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCAc                                                                                                                         >  2:860671/1‑150 (MQ=255)
                            gTTCAGTACGGGAGGAAAGATCGGTGTTACCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCAc                                                                                                                         >  2:812873/1‑150 (MQ=255)
                                 gTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTgcgc                                                                                                                    >  1:690364/1‑150 (MQ=255)
                                   aCGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGccc                                                                                                                  >  1:1612767/1‑150 (MQ=255)
                                      ggaggaAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTc                                                                                                               >  2:187778/1‑150 (MQ=255)
                                      ggaggaAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTc                                                                                                               <  1:880130/150‑1 (MQ=255)
                                         ggaAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAAt                                                                                                            <  1:217129/150‑1 (MQ=255)
                                         ggaAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATATGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAAt                                                                                                            <  1:497574/150‑1 (MQ=255)
                                              gATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCa                                                                                                       >  2:950271/1‑150 (MQ=255)
                                               aTCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGGGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAg                                                                                                      <  1:1536468/150‑1 (MQ=255)
                                                  ggTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGt                                                                                                   <  1:864878/150‑1 (MQ=255)
                                                   gtgtTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTa                                                                                                  >  1:777642/1‑150 (MQ=255)
                                                    tgtTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTAt                                                                                                 >  1:639608/1‑150 (MQ=255)
                                                    tgtTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTAt                                                                                                 >  2:1046739/1‑150 (MQ=255)
                                                        gCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGc                                                                                             >  1:753072/1‑150 (MQ=255)
                                                         cccGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGcc                                                                                            >  2:776126/1‑150 (MQ=255)
                                                          ccgccgCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCa                                                                                           >  2:1191180/1‑150 (MQ=255)
                                                            gccgcAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGAGCCCGTCAATGGTCAGGGTATTCGCCAgg                                                                                         >  2:28443/1‑150 (MQ=255)
                                                               gcAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTa                                                                                      >  2:1180201/1‑150 (MQ=255)
                                                                 aaTTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTacc                                                                                    >  2:1369485/1‑150 (MQ=255)
                                                                     tCACGGGTACCGGCATAGATGGAATCTGAACCTTCGCGGACCTGAGTGACGGTGTCAGTCAAAGAGCGTGGCATATGTGAAACGCTCTGCGACAGGTCGCCCCTTTCCCTGCGCCCGGCAATGGTCAGGGTATTCGCCAGGTTACAACCg                                                                                >  2:41902/1‑150 (MQ=255)
                                                                      cACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCgg                                                                               >  2:734149/1‑150 (MQ=255)
                                                                            tACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGAttt                                                                         >  1:461897/1‑150 (MQ=255)
                                                                             aCCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTc                                                                        <  1:228037/150‑1 (MQ=255)
                                                                              ccGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTcg                                                                       >  2:587692/1‑150 (MQ=255)
                                                                               cGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTcgc                                                                      <  2:690364/150‑1 (MQ=255)
                                                                                ggCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTcgcg                                                                     <  2:401826/150‑1 (MQ=255)
                                                                                 gCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTc                                                                                                                           <  2:441380/95‑1 (MQ=255)
                                                                                 gCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTc                                                                                                                           >  1:441380/1‑95 (MQ=255)
                                                                                   aTAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAAt                                                                  >  1:1114773/1‑150 (MQ=255)
                                                                                      gATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAAtgtg                                                               <  1:1304721/150‑1 (MQ=255)
                                                                                        tGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAg                                                             >  2:430861/1‑150 (MQ=255)
                                                                                        tGGCATCTGAACCTTCGCGGACATGAGTGACGGTGGCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAg                                                             <  1:635921/150‑1 (MQ=255)
                                                                                           cATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCaa                                                          >  1:155994/1‑150 (MQ=255)
                                                                                            aTCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTTAATATGTGAAACGCTCTGCGCCAGGTCGCACATTTCACTGCGCCCGTCAATGTTAAGGGTATTCGCCAGTTTACACCAGGCTATTTCGCGAATGTGAGCaat                                                         >  2:1583904/1‑150 (MQ=255)
                                                                                            aTCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCaat                                                         >  2:636338/1‑150 (MQ=255)
                                                                                                 aaCCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAAttt                                                    <  1:860671/150‑1 (MQ=255)
                                                                                                  aCCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAAtttt                                                   >  1:1551723/1‑150 (MQ=255)
                                                                                                  aCCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAAtttt                                                   <  1:1330145/150‑1 (MQ=255)
                                                                                                        gcgGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAg                                             >  2:1346146/1‑150 (MQ=255)
                                                                                                             cATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAg                                        <  2:1551723/150‑1 (MQ=255)
                                                                                                                gAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTaa                                     <  1:1377824/150‑1 (MQ=255)
                                                                                                                 aGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAg                                    >  2:796929/1‑150 (MQ=255)
                                                                                                                  gTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGc                                   >  2:1346219/1‑150 (MQ=255)
                                                                                                                    gACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCaa                                 >  2:1025459/1‑150 (MQ=255)
                                                                                                                        gtgtCAGTTAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATa                             >  1:1425070/1‑150 (MQ=255)
                                                                                                                        gtgtCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATa                             >  2:683833/1‑150 (MQ=255)
                                                                                                                         tgtCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATAc                            >  2:496761/1‑150 (MQ=255)
                                                                                                                            cagtcaAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGc                         >  2:1508868/1‑150 (MQ=255)
                                                                                                                            cagtcaAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGc                         >  1:1360608/1‑150 (MQ=255)
                                                                                                                            cagtcaAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGc                         <  1:430861/150‑1 (MQ=255)
                                                                                                                             agtcaAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCg                        >  2:1105777/1‑150 (MQ=255)
                                                                                                                              gtcaAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGa                       >  1:1180645/1‑150 (MQ=255)
                                                                                                                              gtcaAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGa                       <  1:450962/150‑1 (MQ=255)
                                                                                                                               gcgcagagCGTTTCACATATGCAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGaa                      >  1:786149/5‑150 (MQ=255)
                                                                                                                                cgcagagCGTTTCACATATGCAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTg                                                 <  2:882056/119‑1 (MQ=255)
                                                                                                                                cgcagagCGTTTCACATATGCAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTg                                                 >  1:882056/4‑122 (MQ=255)
                                                                                                                                cgcagagCGTTTCACATATGCAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAAt                     >  2:406673/4‑150 (MQ=255)
                                                                                                                                cgcagagCGTTTCACATATGCAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAAt                     >  1:1450538/4‑150 (MQ=255)
                                                                                                                                 aaaGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATg                    >  2:1512940/1‑150 (MQ=255)
                                                                                                                                  cagagCGTTTCACATATGCAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGc                   >  1:1398377/2‑150 (MQ=255)
                                                                                                                                  cagagCGTTTCACATATGCAACGCTCAGCGCCAGGTCGCCCATTTCACTGCGCCCGACAATGGTCAGGGTATTCGCCAGGTTACCCCCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGc                   >  2:1471615/2‑150 (MQ=255)
                                                                                                                                  aaGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGc                   >  1:405735/1‑150 (MQ=255)
                                                                                                                                     agCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGt                <  1:1191180/150‑1 (MQ=255)
                                                                                                                                     agCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGt                <  1:1180201/150‑1 (MQ=255)
                                                                                                                                           gCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACg          <  1:776126/150‑1 (MQ=255)
                                                                                                                                           gCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAAAGCCGTACCACg          >  2:742742/1‑150 (MQ=255)
                                                                                                                                                tgtgTAACGCTCTGCGCCAGGTCGCCCATTTCAGTGCGCCGGTCATTGGTAATTGTATTCGCCAGGTTACCACCTGCTATTTCGCTAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACGACGccacc     <  1:777090/150‑1 (MQ=255)
                                                                                                                                                tgtgAAACGCTCTGCGCCATGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGccacc     <  2:907472/150‑1 (MQ=255)
                                                                                                                                                tgtgAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGccacc     <  1:1294812/150‑1 (MQ=255)
                                                                                                                                                 gtgAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGccacca    >  1:1223654/1‑150 (MQ=255)
                                                                                                                                                 gtgAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGccacca    >  1:774788/1‑150 (MQ=255)
                                                                                                                                                   gAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACcagc  >  2:1419319/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGC  >  NC_000913/1971781‑1972077

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