Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,232,484 A→G 5.9% intergenic (+30/+16) umuC → / ← dsbB translesion error‑prone DNA polymerase V subunit; DNA polymerase activity/oxidoreductase that catalyzes reoxidation of DsbA protein disulfide isomerase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,232,4840AG 5.9% 252.5 / 5.0 84intergenic (+30/+16)umuC/dsbBtranslesion error‑prone DNA polymerase V subunit; DNA polymerase activity/oxidoreductase that catalyzes reoxidation of DsbA protein disulfide isomerase I
Reads supporting (aligned to +/- strand):  ref base A (42/37);  new base G (2/3);  total (44/40)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.75e-01

GAGGGCAGAGGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCC  >  NC_000913/1232335‑1232627
                                                                                                                                                     |                                                                                                                                               
gagGGCAGAGGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCa                                                                                                                                                 <  2:1447233/150‑1 (MQ=255)
  gggCAGAGGAACACTCTATTTTGCCGGGCAGGGGATCCCGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAgg                                                                                                                                               >  1:243005/1‑150 (MQ=255)
   ggCAGAGGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCCCGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGa                                                                                                                                              >  1:677170/1‑150 (MQ=255)
    gCAGAGGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAGCGCTCCCGCAGGAgc                                                                                                                                            >  1:946438/1‑150 (MQ=255)
        agGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgct                                                                                                                                         <  2:1462992/150‑2 (MQ=255)
           aaCACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTGTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGaa                                                                                                                                      >  2:957887/1‑150 (MQ=255)
           aaCACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACCTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgctttt                                                                                                                                      <  1:1245375/150‑5 (MQ=255)
            acacTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGGCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAAt                                                                                                                                     <  2:950748/150‑1 (MQ=255)
             cacTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATg                                                                                                                                    <  2:584211/150‑1 (MQ=255)
              acTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATgg                                                                                                                                   >  1:475662/1‑150 (MQ=255)
                tctATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATTAATATAGCTGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGAt                                                                                                                                 <  1:1303559/150‑1 (MQ=255)
                 ctATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGAtt                                                                                                                                >  2:75964/1‑150 (MQ=255)
                   aTTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAg                                                                                                                              <  2:1420962/150‑1 (MQ=255)
                   aTTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAg                                                                                                                              <  2:14138/150‑1 (MQ=255)
                    ttttGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGc                                                                                                                             <  1:1177493/150‑1 (MQ=255)
                        gCCGGGCAGGGGATCCATCAACAATGGCAGATGAAGCGAGACATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGATGGTCTAATAAATATAGCGGCAGGAAAAAAGCGCGCCCGCAGGAGCGCCGAATGGATAAGCGAcc                                                                                                                         <  2:540293/150‑1 (MQ=255)
                         ccGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCg                                                                                                                        <  2:496809/150‑1 (MQ=255)
                         ccGGGCAGGGGATCCAGCAACAATGGAAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCg                                                                                                                        >  1:457920/1‑150 (MQ=255)
                           gggCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGTTTAGCGACCGaa                                                                                                                      >  2:357266/1‑150 (MQ=255)
                           gggCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGaa                                                                                                                      >  1:1789956/1‑150 (MQ=255)
                            ggCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAAc                                                                                                                     <  2:653563/150‑1 (MQ=255)
                               aGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCAGGTTATACAACGCGATGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGTGCGCCGAATGGATTAGCGACCGAACAGa                                                                                                                  <  2:1171171/150‑1 (MQ=255)
                                ggggATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGAt                                                                                                                 >  2:482101/1‑150 (MQ=255)
                                ggggATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGAt                                                                                                                 >  2:386793/1‑150 (MQ=255)
                                  ggATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAAAGGATTAGCGACCGAACAGATCa                                                                                                               >  2:1796941/1‑150 (MQ=255)
                                   gATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCAc                                                                                                              <  2:1596618/150‑1 (MQ=255)
                                       cAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAGAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGttt                                                                                                          >  1:1657257/1‑150 (MQ=255)
                                          caacaaTGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGCTCACGtttttt                                                                                                       >  1:455418/1‑150 (MQ=255)
                                             aaaTGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGc                                                                                                    <  1:1494401/149‑1 (MQ=255)
                                                 ggCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTa                                                                                                <  1:1212239/150‑1 (MQ=255)
                                                  gCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTATCGATTTaa                                                                                               >  1:1113841/1‑150 (MQ=255)
                                                      aTGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTAGCTTTAAACgg                                                                                           >  2:1452296/1‑150 (MQ=255)
                                                       tGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGc                                                                                          >  1:833073/1‑150 (MQ=255)
                                                         aaGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTg                                                                                        <  1:1340589/150‑1 (MQ=255)
                                                         aaGCGAGCAATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTTCTGAGTGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGATCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTg                                                                                        <  1:1104693/150‑1 (MQ=255)
                                                          aGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTgg                                                                                       >  2:962934/1‑150 (MQ=255)
                                                          aGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTgg                                                                                       <  1:1142907/150‑1 (MQ=255)
                                                          aGCGAGCCATGCATTCACCACGTAATACAAAGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTgg                                                                                       <  1:1197532/150‑1 (MQ=255)
                                                           gCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGTAAAAAAGAGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTggg                                                                                      <  2:1112621/150‑1 (MQ=255)
                                                           gCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTggg                                                                                      <  2:682126/150‑1 (MQ=255)
                                                               gCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGCCCGAACAGATCACGTTTTTTCGCTTTAAACGGCGTGGAAAt                                                                                  >  2:1784616/1‑150 (MQ=255)
                                                               gCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAAt                                                                                  >  1:1706945/1‑150 (MQ=255)
                                                                 cATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATca                                                                                <  2:1257963/150‑1 (MQ=255)
                                                                     cTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCTGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATcaccac                                                                            <  1:1189874/150‑1 (MQ=255)
                                                                      tttCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCTCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCACTTTAAACGGCTGGGAAATCaccacc                                                                           <  1:1656507/150‑1 (MQ=255)
                                                                      tttCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCaccacc                                                                           <  1:1723153/150‑1 (MQ=255)
                                                                         caccacGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAACTAAATATTGCGGCAGGAAAAAAGCACTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAAAATATCACGTTTATTCTCTTTAAACGGCTGTGAAATAACAACAAGc                                                                        >  1:362387/1‑150 (MQ=255)
                                                                          accacGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTGTAAACGGCGGGGAAATCACCACCAGCa                                                                       >  1:723525/1‑150 (MQ=255)
                                                                           ccacGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTAAAACGGCTGGGAAATCACCACCAGCAc                                                                      >  1:1673754/1‑150 (MQ=255)
                                                                            cacGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACt                                                                     >  2:1333599/1‑150 (MQ=255)
                                                                              cGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGc                                                                   >  2:226482/1‑150 (MQ=255)
                                                                              cGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCg                                                                  <  2:946438/150‑1 (MQ=255)
                                                                                     cAACGCGAAGTTCTGATTTACTGAGGTTCAAATAAATATAGCGCCAGGAAAACAGCGCTCCCGCAGGAGCGCAGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCAACACCAGCACTGCGACAATc                                                            <  1:905853/150‑1 (MQ=255)
                                                                                       aCGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCCCCAGCACTGCGACAATCAg                                                          >  2:1340276/1‑150 (MQ=255)
                                                                                        cgcgAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAgg                                                         <  1:587526/150‑1 (MQ=255)
                                                                                            aaGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAg                                                     <  2:1706945/150‑1 (MQ=255)
                                                                                             aGTTATGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAACCGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGc                                                    <  1:217329/150‑1 (MQ=255)
                                                                                              gTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCg                                                   >  2:1275292/1‑150 (MQ=255)
                                                                                              gTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCg                                                   >  2:379224/1‑150 (MQ=255)
                                                                                               ttCTGATATACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGa                                                  <  2:1766173/150‑1 (MQ=255)
                                                                                                tCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCATGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATAACCACCAGCACTGCGACAATCAGGTAAGCGAt                                                 >  2:1096286/1‑150 (MQ=255)
                                                                                                 cTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGGCAATCTGGTAAGCGATa                                                >  2:1166021/1‑150 (MQ=255)
                                                                                                  tGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATaa                                               >  1:682743/1‑150 (MQ=255)
                                                                                                  tGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATaa                                               >  2:1494998/1‑150 (MQ=255)
                                                                                                   gATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATaaa                                              >  1:789943/1‑150 (MQ=255)
                                                                                                   gATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATaaa                                              <  1:357266/150‑1 (MQ=255)
                                                                                                      ttACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAAt                                           <  2:677170/150‑1 (MQ=255)
                                                                                                         cTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATAcc                                        >  1:535053/1‑150 (MQ=255)
                                                                                                           gAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGACTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGa                                      <  2:1789956/150‑1 (MQ=255)
                                                                                                             gggTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGagc                                    <  2:457920/150‑1 (MQ=255)
                                                                                                              ggTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATCCCGagca                                   >  2:1463662/1‑150 (MQ=255)
                                                                                                                 cAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTGTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGcc                                >  1:672683/1‑150 (MQ=255)
                                                                                                                  aaataaatATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGccc                               >  1:664676/1‑149 (MQ=255)
                                                                                                                   aataaataTAGCGGCAGGAAAAAAGCGCACCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCAc                              >  2:1160585/1‑150 (MQ=255)
                                                                                                                       aataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACAACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCg                          >  1:771973/1‑150 (MQ=255)
                                                                                                                             gCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCAttt                    >  1:318213/1‑150 (MQ=255)
                                                                                                                               ggCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGACATTTcc                  <  1:1511794/150‑1 (MQ=255)
                                                                                                                                   ggAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGAc              <  1:1494998/150‑1 (MQ=255)
                                                                                                                                    gAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGAcc             <  2:1353967/150‑1 (MQ=255)
                                                                                                                                       ttcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGCACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAAGCACCACCAGCACTGCGGCAATCAGGTAAGCGATAATAATTCCGAGCCGCCACTGCGGCAGTTCCGGACCTaa          >  1:1638923/5‑150 (MQ=255)
                                                                                                                                         aaGCGCTCCCGCAGGAGCGCCGAATCGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCGACTGCGGCATTTCCAGACCTaaaa        >  2:71215/1‑150 (MQ=255)
                                                                                                                                          ggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTaaaaa       >  1:647311/2‑150 (MQ=255)
                                                                                                                                          aGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAATTCACCACCAGCACTGCGACAATCAGGTAAGCGATATAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTaacaa       >  2:639055/1‑150 (MQ=255)
                                                                                                                                             gcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATc    <  2:647311/150‑1 (MQ=255)
                                                                                                                                               tCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATccc  <  1:957887/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
GAGGGCAGAGGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCC  >  NC_000913/1232335‑1232627

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: