Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,477,022 T→G 7.4% D203A (GAC→GCC)  dsdC ← dsd operon activator; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,477,0220TG 7.4% 181.8 / 3.6 76D203A (GAC→GCC) dsdCdsd operon activator; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (13/56);  new base G (4/2);  total (17/59)
Fisher's exact test for biased strand distribution p-value = 2.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATC  >  NC_000913/2476879‑2477171
                                                                                                                                               |                                                                                                                                                      
ccccAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGa                                                                                                                                                  <  1:1809051/150‑1 (MQ=255)
 cccAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAg                                                                                                                                                 <  2:426866/150‑1 (MQ=255)
     atatCATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGACAGACAGGCATGGAGCAAc                                                                                                                                             <  2:197362/150‑1 (MQ=255)
        tGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCTATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGACTGTATGTCAAGGAGCaaacaa                                                                                                                                          >  1:72695/1‑146 (MQ=255)
        tGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTa                                                                                                                                          <  2:1095677/150‑1 (MQ=255)
             cATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATTCCTGTCTGTCATGGAGCAACGTACAGTg                                                                                                                                     >  2:1082178/1‑150 (MQ=255)
               tCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGac                                                                                                                                   <  2:1718648/150‑1 (MQ=255)
                   ggCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAg                                                                                                                               <  2:1511090/150‑1 (MQ=255)
                    gCAATAACAGCTAAATCAGAAAGATCAAAGCCAATTCCAGACGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCTTACCGGAGTCGTTGCTCCATGCCGGTCTGTCATGGAGCAACGTACAGTGACTCAgg                                                                                                                              <  1:946401/150‑1 (MQ=255)
                           cAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTa                                                                                                                      <  1:557873/150‑1 (MQ=255)
                             gCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTAcc                                                                                                                    <  1:369024/150‑1 (MQ=255)
                              cTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCg                                                                                                                   <  1:260388/150‑1 (MQ=255)
                              cTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCg                                                                                                                   <  2:1309700/150‑1 (MQ=255)
                               tAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCgt                                                                                                                  >  2:278887/1‑150 (MQ=255)
                               tAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCAACATTCCAGTATGTGCTGGAGGCACGTAAAATGCCCCCAG‑TTAATTAACgt                                                                                                                  >  1:1119806/1‑150 (MQ=255)
                                aaaTCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCgtg                                                                                                                 <  1:697619/150‑1 (MQ=255)
                                   tCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTg                                                                                                              <  2:116555/150‑1 (MQ=255)
                                    cAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTgg                                                                                                             <  1:909133/150‑1 (MQ=255)
                                     aGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCAAACTATGCCATTCATCCGTACCGGAGTCGTTTATACATGCCTATAAGTCAAGGTGCAACGTAACCTGACACAGGTTTAATTACC‑TGTTGGt                                                                                                            >  1:1102725/1‑150 (MQ=255)
                                        aCGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTaa                                                                                                         <  2:755995/150‑1 (MQ=255)
                                         cGATCAAAGCCAATTCCAGAAGCTGTCGGCAAATTAACCGCATACTGTTGCGCCCAACTATGCCATTCAACCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTaaa                                                                                                        <  2:1955138/150‑1 (MQ=255)
                                          gATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAg                                                                                                       <  2:491286/150‑1 (MQ=255)
                                          gATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAg                                                                                                       <  1:1299552/150‑1 (MQ=255)
                                          gATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAg                                                                                                       <  1:607116/150‑1 (MQ=255)
                                              aaaGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATg                                                                                                   <  2:605878/150‑1 (MQ=255)
                                                      ttCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAg                                                                                           <  1:1467309/150‑1 (MQ=255)
                                                      ttCCAGAAGATGTCGGCAAAGTAACCGCATAATGTTGCGCCCAACTCTGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGGCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAg                                                                                           <  1:1282285/150‑1 (MQ=255)
                                                        ccAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCg                                                                                         <  1:162552/150‑1 (MQ=255)
                                                                tGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACATACAGTGACACAGG‑GTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCgg                                                                                 <  1:1449482/150‑1 (MQ=255)
                                                                 gtcgtcAAATTAACCGCATAATGTTGCGCCCAACTATGACATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCggg                                                                                <  2:474408/150‑1 (MQ=255)
                                                                 gTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCggg                                                                                <  1:1003466/150‑1 (MQ=255)
                                                                     gCAAATTAACCGCATAATGTTGCGCCCATCTATGCCATTCATCCGTTCCGTAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTAc                                                                            <  1:1684242/150‑1 (MQ=255)
                                                                      cAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACa                                                                           <  1:203708/150‑1 (MQ=255)
                                                                             aCCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCCTGCCGGTGGGTCATGGAGCAACGGACAGTGACACAGG‑GTAATTACCCTGTTGGGTAAAGCATGTCTTTGAGCGGATTCCGGGCGACATACTggg                                                                    >  2:1164702/1‑149 (MQ=255)
                                                                               cGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAg                                                                  <  2:243132/150‑1 (MQ=255)
                                                                                      tGTTGCTCCCAACTATGACATTCATCCGTCCCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATttct                                                           <  2:1989981/150‑1 (MQ=255)
                                                                                       tttGCGCCCAACTATGCCTTTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATttctt                                                          <  1:720037/149‑1 (MQ=255)
                                                                                          gcgcCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTcat                                                       <  1:721983/150‑1 (MQ=255)
                                                                                           cgcCCAACTATGCCATTCATCCGTACCGGAGTCGGTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTcatc                                                      <  2:1221377/150‑1 (MQ=255)
                                                                                             cccAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTcatcca                                                    <  1:1962763/150‑1 (MQ=255)
                                                                                               cAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTcatccatc                                                  <  1:1703842/150‑1 (MQ=255)
                                                                                                      gCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGt                                           <  1:390809/150‑1 (MQ=255)
                                                                                                        cATTCATCCGTACCGGAGTCGTTTCTCCATGCCGGTCTGTCATGGGGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATTTCTTTTAGCGTATTCCGGGCTACATACAGGCAGGATTTGTTCATCCAACAGAAAGtga                                         >  2:1027054/1‑150 (MQ=255)
                                                                                                          ttCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGtgatg                                       <  2:1427019/150‑1 (MQ=255)
                                                                                                          ttCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGGCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGtgatg                                       <  2:993207/150‑1 (MQ=255)
                                                                                                          ttCATCCGTAACGGAGTCGTTGCTCCATGCCTCTGTGTCATGGAGCAACGGACATTGACACAGG‑ATAATTAGCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAACGtgatg                                       <  2:324209/150‑1 (MQ=255)
                                                                                                            cATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGAtgag                                     <  2:245172/150‑1 (MQ=255)
                                                                                                              tCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGAtgagtg                                   <  1:1947153/150‑1 (MQ=255)
                                                                                                              tCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGAtgagtg                                   <  1:511324/150‑1 (MQ=255)
                                                                                                                       ccggcgctGCTCCATGACTCTCTTCCCTCGAGCAACGTACAGTGACACCTC‑ATAATAACCTTGTTGGTTACAGCATGTCTCTGAGGGTATTCCGGGCTACATACTCGCAGGAATTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCt                          <  2:345582/143‑1 (MQ=255)
                                                                                                                         gTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGa                        >  1:45561/1‑150 (MQ=255)
                                                                                                                           cGTTGCTCCGGGCCGGGCGGGCATGGTGCGGCGGAGGGTGACACAGG‑GTAATTACCGTGTTGGTTGAAGCATGTGTTTGAGCGGGTTGCGGGGTACGTACTGGCAGGATTTCTTCATCCATCAGCAAGTGATGTGTGCGTTGCGCGGACg                      >  2:2009076/1‑150 (MQ=255)
                                                                                                                           cGTTGCTCCATTCCTGTCTGTCATGGAGCAACGTACAGTGACACTGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCTGGCTACATACTGTCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACg                      <  2:199642/150‑1 (MQ=255)
                                                                                                                           cGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACg                      <  2:242661/150‑1 (MQ=255)
                                                                                                                           cGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑GTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACg                      <  2:1889282/150‑1 (MQ=255)
                                                                                                                            gTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACgg                     >  1:915921/1‑150 (MQ=255)
                                                                                                                             ttGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgc                    >  2:1374224/1‑150 (MQ=255)
                                                                                                                              tGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcg                   >  1:1318750/1‑150 (MQ=255)
                                                                                                                               gCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGGCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcgc                  <  1:278887/150‑1 (MQ=255)
                                                                                                                               gCTCCATGACAGACAGGCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGCGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcgc                  >  2:879958/1‑150 (MQ=255)
                                                                                                                               gCTCCATGACAGACAGGCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcgc                  >  2:1895662/1‑150 (MQ=255)
                                                                                                                               gCTCCATGACAGACAGGCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcgc                  >  1:369957/1‑150 (MQ=255)
                                                                                                                                cTCCATGCCTGTCTGTCATGGAGCAACGTACAGAGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGca                 <  2:743433/150‑1 (MQ=255)
                                                                                                                                 tCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGcat                <  2:1815982/150‑1 (MQ=255)
                                                                                                                                 gccATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGcat                <  1:566571/149‑1 (MQ=255)
                                                                                                                                  ccATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGcatc               <  1:504362/150‑1 (MQ=255)
                                                                                                                                   cATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGcatca              <  1:1082178/150‑1 (MQ=255)
                                                                                                                                   cATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGcatca              <  2:1858493/150‑1 (MQ=255)
                                                                                                                                   cATGCCTGTCTGTCATGGAGCAACGTAAAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGcatca              <  2:72695/150‑1 (MQ=255)
                                                                                                                                     tGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCatcatc            <  1:820097/150‑1 (MQ=255)
                                                                                                                                     tGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCatcatc            <  1:1944663/150‑1 (MQ=255)
                                                                                                                                     tGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCCGTCAGCGcaactca            <  2:1718035/150‑6 (MQ=255)
                                                                                                                                       ccTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCaa          <  2:1907144/150‑1 (MQ=255)
                                                                                                                                       ccTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCaa          <  1:663593/150‑1 (MQ=255)
                                                                                                                                       ccTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCaa          >  1:1212826/1‑150 (MQ=255)
                                                                                                                                       ccTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGACATTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTACATCCAACACAAAGTGATGAGTGAGTTGCGCTGACGGAGCATAATCaa          >  1:1219157/1‑150 (MQ=255)
                                                                                                                                           gctgtcATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGtt      <  1:346018/149‑2 (MQ=255)
                                                                                                                                            ctgtcATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCAGGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAg     <  1:527469/150‑1 (MQ=255)
                                                                                                                                             tgtcATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGa    <  1:847911/150‑1 (MQ=255)
                                                                                                                                             tgtcATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGGAAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTGCATCCATCAGACAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGa    <  1:1801478/150‑1 (MQ=255)
                                                                                                                                               tcATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGTGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATc  <  1:1526309/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                      
CCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATC  >  NC_000913/2476879‑2477171

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