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

CCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATC  >  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)
                           cAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTa                                                                                                                      <  1:557873/150‑1 (MQ=255)
                             gCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTAcc                                                                                                                    <  1:369024/150‑1 (MQ=255)
                              cTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCg                                                                                                                   <  1:260388/150‑1 (MQ=255)
                              cTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCg                                                                                                                   <  2:1309700/150‑1 (MQ=255)
                               tAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCgt                                                                                                                  >  2:278887/1‑150 (MQ=255)
                               tAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCAACATTCCAGTATGTGCTGGAGGCACGTAAAATGCCCCCAGTTAATTAACgt                                                                                                                  >  1:1119806/1‑150 (MQ=255)
                                aaaTCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCgtg                                                                                                                 <  1:697619/150‑1 (MQ=255)
                                   tCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTg                                                                                                              <  2:116555/150‑1 (MQ=255)
                                    cAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTgg                                                                                                             <  1:909133/150‑1 (MQ=255)
                                     aGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCAAACTATGCCATTCATCCGTACCGGAGTCGTTTATACATGCCTATAAGTCAAGGTGCAACGTAACCTGACACAGGTTTAATTACCTGTTGGt                                                                                                            >  1:1102725/1‑150 (MQ=255)
                                        aCGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTaa                                                                                                         <  2:755995/150‑1 (MQ=255)
                                         cGATCAAAGCCAATTCCAGAAGCTGTCGGCAAATTAACCGCATACTGTTGCGCCCAACTATGCCATTCAACCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTaaa                                                                                                        <  2:1955138/150‑1 (MQ=255)
                                          gATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAg                                                                                                       <  2:491286/150‑1 (MQ=255)
                                          gATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAg                                                                                                       <  1:1299552/150‑1 (MQ=255)
                                          gATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAg                                                                                                       <  1:607116/150‑1 (MQ=255)
                                              aaaGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATg                                                                                                   <  2:605878/150‑1 (MQ=255)
                                                      ttCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAg                                                                                           <  1:1467309/150‑1 (MQ=255)
                                                      ttCCAGAAGATGTCGGCAAAGTAACCGCATAATGTTGCGCCCAACTCTGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGGCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAg                                                                                           <  1:1282285/150‑1 (MQ=255)
                                                        ccAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCg                                                                                         <  1:162552/150‑1 (MQ=255)
                                                                tGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACATACAGTGACACAGGGTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCgg                                                                                 <  1:1449482/150‑1 (MQ=255)
                                                                 gtcgtcAAATTAACCGCATAATGTTGCGCCCAACTATGACATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCggg                                                                                <  2:474408/150‑1 (MQ=255)
                                                                 gTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCggg                                                                                <  1:1003466/150‑1 (MQ=255)
                                                                     gCAAATTAACCGCATAATGTTGCGCCCATCTATGCCATTCATCCGTTCCGTAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTAc                                                                            <  1:1684242/150‑1 (MQ=255)
                                                                      cAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACa                                                                           <  1:203708/150‑1 (MQ=255)
                                                                             aCCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCCTGCCGGTGGGTCATGGAGCAACGGACAGTGACACAGGGTAATTACCCTGTTGGGTAAAGCATGTCTTTGAGCGGATTCCGGGCGACATACTggg                                                                    >  2:1164702/1‑149 (MQ=255)
                                                                               cGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAg                                                                  <  2:243132/150‑1 (MQ=255)
                                                                                      tGTTGCTCCCAACTATGACATTCATCCGTCCCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATttct                                                           <  2:1989981/150‑1 (MQ=255)
                                                                                       tttGCGCCCAACTATGCCTTTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATttctt                                                          <  1:720037/149‑1 (MQ=255)
                                                                                          gcgcCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTcat                                                       <  1:721983/150‑1 (MQ=255)
                                                                                           cgcCCAACTATGCCATTCATCCGTACCGGAGTCGGTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTcatc                                                      <  2:1221377/150‑1 (MQ=255)
                                                                                             cccAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTcatcca                                                    <  1:1962763/150‑1 (MQ=255)
                                                                                               cAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTcatccatc                                                  <  1:1703842/150‑1 (MQ=255)
                                                                                                      gCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGt                                           <  1:390809/150‑1 (MQ=255)
                                                                                                        cATTCATCCGTACCGGAGTCGTTTCTCCATGCCGGTCTGTCATGGGGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATTTCTTTTAGCGTATTCCGGGCTACATACAGGCAGGATTTGTTCATCCAACAGAAAGtga                                         >  2:1027054/1‑150 (MQ=255)
                                                                                                          ttCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGtgatg                                       <  2:1427019/150‑1 (MQ=255)
                                                                                                          ttCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGGCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGtgatg                                       <  2:993207/150‑1 (MQ=255)
                                                                                                          ttCATCCGTAACGGAGTCGTTGCTCCATGCCTCTGTGTCATGGAGCAACGGACATTGACACAGGATAATTAGCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAACGtgatg                                       <  2:324209/150‑1 (MQ=255)
                                                                                                            cATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGAtgag                                     <  2:245172/150‑1 (MQ=255)
                                                                                                              tCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGAtgagtg                                   <  1:1947153/150‑1 (MQ=255)
                                                                                                              tCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGAtgagtg                                   <  1:511324/150‑1 (MQ=255)
                                                                                                                       ccggcgctGCTCCATGACTCTCTTCCCTCGAGCAACGTACAGTGACACCTCATAATAACCTTGTTGGTTACAGCATGTCTCTGAGGGTATTCCGGGCTACATACTCGCAGGAATTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCt                          <  2:345582/143‑1 (MQ=255)
                                                                                                                         gTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGa                        >  1:45561/1‑150 (MQ=255)
                                                                                                                           cGTTGCTCCGGGCCGGGCGGGCATGGTGCGGCGGAGGGTGACACAGGGTAATTACCGTGTTGGTTGAAGCATGTGTTTGAGCGGGTTGCGGGGTACGTACTGGCAGGATTTCTTCATCCATCAGCAAGTGATGTGTGCGTTGCGCGGACg                      >  2:2009076/1‑150 (MQ=255)
                                                                                                                           cGTTGCTCCATTCCTGTCTGTCATGGAGCAACGTACAGTGACACTGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCTGGCTACATACTGTCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACg                      <  2:199642/150‑1 (MQ=255)
                                                                                                                           cGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACg                      <  2:242661/150‑1 (MQ=255)
                                                                                                                           cGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGGTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACg                      <  2:1889282/150‑1 (MQ=255)
                                                                                                                            gTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACgg                     >  1:915921/1‑150 (MQ=255)
                                                                                                                             ttGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgc                    >  2:1374224/1‑150 (MQ=255)
                                                                                                                              tGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcg                   >  1:1318750/1‑150 (MQ=255)
                                                                                                                               gCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGGCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcgc                  <  1:278887/150‑1 (MQ=255)
                                                                                                                               gCTCCATGACAGACAGGCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGCGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcgc                  >  2:879958/1‑150 (MQ=255)
                                                                                                                               gCTCCATGACAGACAGGCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcgc                  >  2:1895662/1‑150 (MQ=255)
                                                                                                                               gCTCCATGACAGACAGGCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcgc                  >  1:369957/1‑150 (MQ=255)
                                                                                                                                cTCCATGCCTGTCTGTCATGGAGCAACGTACAGAGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGca                 <  2:743433/150‑1 (MQ=255)
                                                                                                                                 tCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGcat                <  2:1815982/150‑1 (MQ=255)
                                                                                                                                 gccATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGcat                <  1:566571/149‑1 (MQ=255)
                                                                                                                                  ccATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGcatc               <  1:504362/150‑1 (MQ=255)
                                                                                                                                   cATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGcatca              <  1:1082178/150‑1 (MQ=255)
                                                                                                                                   cATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGcatca              <  2:1858493/150‑1 (MQ=255)
                                                                                                                                   cATGCCTGTCTGTCATGGAGCAACGTAAAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGcatca              <  2:72695/150‑1 (MQ=255)
                                                                                                                                     tGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCatcatc            <  1:820097/150‑1 (MQ=255)
                                                                                                                                     tGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCatcatc            <  1:1944663/150‑1 (MQ=255)
                                                                                                                                     tGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCCGTCAGCGcaactca            <  2:1718035/150‑6 (MQ=255)
                                                                                                                                       ccTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCaa          <  2:1907144/150‑1 (MQ=255)
                                                                                                                                       ccTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCaa          <  1:663593/150‑1 (MQ=255)
                                                                                                                                       ccTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCaa          >  1:1212826/1‑150 (MQ=255)
                                                                                                                                       ccTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGACATTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTACATCCAACACAAAGTGATGAGTGAGTTGCGCTGACGGAGCATAATCaa          >  1:1219157/1‑150 (MQ=255)
                                                                                                                                           gctgtcATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGtt      <  1:346018/149‑2 (MQ=255)
                                                                                                                                            ctgtcATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCAGGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAg     <  1:527469/150‑1 (MQ=255)
                                                                                                                                             tgtcATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGa    <  1:847911/150‑1 (MQ=255)
                                                                                                                                             tgtcATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGGAAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTGCATCCATCAGACAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGa    <  1:1801478/150‑1 (MQ=255)
                                                                                                                                               tcATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGTGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATc  <  1:1526309/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                      
CCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATC  >  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: