Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,772,237 A→G 9.5% pseudogene (946/1016 nt) yfjV ← CP4‑57 prophage; putative arsenite transporter, N‑terminal fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,772,2370AG 9.5% 299.4 / 17.6 105pseudogene (946/1016 nt)yfjVCP4‑57 prophage; putative arsenite transporter, N‑terminal fragment
Reads supporting (aligned to +/- strand):  ref base A (43/52);  new base G (5/5);  total (48/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.49e-01

Only 100 of 108 total aligned reads displayed.

GAAGCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTG  >  NC_000913/2772088‑2772386
                                                                                                                                                     |                                                                                                                                                      
gAAGCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACaa                                                                                                                                                        >  1:422816/1‑150 (MQ=255)
    cATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGcc                                                                                                                                                    <  1:1186678/150‑1 (MQ=255)
     aTCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCg                                                                                                                                                   <  2:11955/150‑1 (MQ=255)
       ctcCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCGTGACAATGCCGGt                                                                                                                                                 >  1:1519624/1‑150 (MQ=255)
        tcCAGCGTGTTGCGTGAGGTGCCTCAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCATTGTCATGACTGTGCCGGtg                                                                                                                                                <  2:316235/150‑2 (MQ=21)
          cAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCATTGTCATGACTGTGCCGGtgct                                                                                                                                              <  1:1458121/150‑4 (MQ=21)
          cAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACg                                                                                                                                              <  2:897235/150‑1 (MQ=255)
              gtgtTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACTCCAGCGCGGCCAGAGTGACAAACAGCACCGGCATTGTCATGACTGTGCCGGtgctgttt                                                                                                                                          <  2:907138/150‑8 (MQ=21)
              gtgtTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAg                                                                                                                                          <  1:173062/150‑1 (MQ=255)
                 ttGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTaa                                                                                                                                       >  1:330695/1‑150 (MQ=255)
                    cGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAatat                                                                                                                                    >  1:1049759/1‑150 (MQ=255)
                    cGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAatat                                                                                                                                    >  2:952130/1‑150 (MQ=255)
                         ggTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCa                                                                                                                               >  1:1155116/1‑150 (MQ=255)
                         ggTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCa                                                                                                                               >  2:1502236/1‑150 (MQ=255)
                         ggTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCa                                                                                                                               <  2:1626482/150‑1 (MQ=255)
                         ggTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCa                                                                                                                               >  2:363583/1‑150 (MQ=255)
                            gCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGt                                                                                                                            >  1:1297518/1‑150 (MQ=255)
                             ccACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCATTGTCATGACTGTGCCGGtgctg                                                                                                                                             >  2:1561535/1‑128 (MQ=18)
                             ccACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCATTGTCATGACTGTGCCGGtgctg                                                                                                                                             <  1:1561535/132‑5 (MQ=18)
                             ccACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTg                                                                                                                           >  1:1142834/1‑150 (MQ=255)
                                caGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGAtt                                                                                                                        >  2:46231/1‑150 (MQ=255)
                                caGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGAtt                                                                                                                        <  2:422816/150‑1 (MQ=255)
                                  ggCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACgg                                                                                                                                             <  2:723032/127‑1 (MQ=255)
                                  ggCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACgg                                                                                                                                             >  1:723032/1‑127 (MQ=255)
                                   gCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTc                                                                                                                     <  1:46231/150‑1 (MQ=255)
                                       ggTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCatat                                                                                                                 >  1:732462/1‑150 (MQ=255)
                                        gTTGTCATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATAtt                                                                                                                <  1:443987/150‑1 (MQ=255)
                                           gtgATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATAttttt                                                                                                             >  1:542757/1‑150 (MQ=255)
                                              aTAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGTCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTg                                                                                                          <  1:1125589/150‑1 (MQ=255)
                                              aTAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTg                                                                                                          <  2:158002/150‑1 (MQ=255)
                                               tAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTgg                                                                                                         <  1:351784/150‑1 (MQ=255)
                                                aGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTggg                                                                                                        <  1:1208549/150‑1 (MQ=255)
                                                aGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTggg                                                                                                        <  1:972824/150‑1 (MQ=255)
                                                aGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTggg                                                                                                        <  2:1297518/150‑1 (MQ=255)
                                                    ggTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAg                                                                                                    <  2:330695/150‑1 (MQ=255)
                                                     gTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGc                                                                                                   >  2:1166680/1‑150 (MQ=255)
                                                      tAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGca                                                                                                  <  1:912543/150‑1 (MQ=255)
                                                          gTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATg                                                                                              <  2:567523/150‑1 (MQ=255)
                                                           ttGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGc                                                                                             <  2:890739/150‑1 (MQ=255)
                                                            tGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCa                                                                                            <  1:123804/150‑1 (MQ=255)
                                                             gCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAg                                                                                           >  2:332578/1‑150 (MQ=255)
                                                                cATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCca                                                                                        >  1:433395/1‑150 (MQ=255)
                                                                cATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCca                                                                                        <  2:1175165/150‑1 (MQ=255)
                                                                 atatCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCcac                                                                                       <  2:560638/150‑1 (MQ=255)
                                                                  tatCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCcaca                                                                                      <  2:1341377/150‑1 (MQ=255)
                                                                  tatCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCcaca                                                                                      >  2:175965/1‑150 (MQ=255)
                                                                  tatCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCcaca                                                                                      >  1:1273329/1‑150 (MQ=255)
                                                                    tCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCGTGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCAcagc                                                                                    <  2:1519624/150‑1 (MQ=255)
                                                                    tCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCAcagc                                                                                    <  2:1412121/150‑1 (MQ=255)
                                                                     cAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCAcagca                                                                                   >  2:1134194/1‑150 (MQ=255)
                                                                       gTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCaga                                                                                 >  2:173568/1‑150 (MQ=255)
                                                                        tATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCagag                                                                                <  1:332028/150‑1 (MQ=255)
                                                                           cTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTg                                                                             <  2:937987/150‑1 (MQ=255)
                                                                           cTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTg                                                                             >  1:490265/1‑150 (MQ=255)
                                                                           cTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTg                                                                             >  1:1308276/1‑150 (MQ=255)
                                                                             aaTTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGcc                                                                           >  1:1030298/1‑150 (MQ=255)
                                                                             aaTTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGcc                                                                           >  2:539655/1‑150 (MQ=255)
                                                                              aTTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCa                                                                          <  2:396373/150‑1 (MQ=255)
                                                                               ttACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAg                                                                         <  2:732462/150‑1 (MQ=255)
                                                                                 aCAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGAc                                                                       >  1:843728/1‑150 (MQ=255)
                                                                                    aaGTGAGAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAAGAGCACCGGCAGAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGc                                                                    <  1:1409519/150‑1 (MQ=255)
                                                                                     aGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGcc                                                                   <  2:1142834/150‑1 (MQ=255)
                                                                                     aGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGcc                                                                   <  1:746752/150‑1 (MQ=255)
                                                                                      gTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCGCAAGTGATTGTGATATTTTTCTGGGAAAGCATATGCAGCGACAGCAGAGTTGCCAGACTGCCa                                                                  <  2:1049759/150‑1 (MQ=255)
                                                                                          cAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCg                                                              >  1:282244/1‑150 (MQ=255)
                                                                                          cAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCg                                                              >  1:1162023/1‑150 (MQ=255)
                                                                                          cAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCACCCACAGCAGAGTTGCCAGACTGCCAATCg                                                              >  1:703086/1‑150 (MQ=255)
                                                                                              gagCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGt                                                         >  1:980527/1‑150 (MQ=255)
                                                                                               agCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGTCAATCGGGGTg                                                        <  2:680838/150‑1 (MQ=255)
                                                                                               agCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTg                                                        <  1:18215/150‑1 (MQ=255)
                                                                                               agCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTg                                                        >  2:295732/1‑150 (MQ=255)
                                                                                                   gccacgccaGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGAttt                                                    >  1:483824/1‑150 (MQ=255)
                                                                                                       cgccaGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCgg                                                <  2:422883/150‑1 (MQ=255)
                                                                                                       cgccaGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCgg                                                <  1:1380492/150‑1 (MQ=255)
                                                                                                          caGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGcc                                             <  1:1137671/150‑1 (MQ=255)
                                                                                                          caGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGcc                                             <  2:296388/150‑1 (MQ=255)
                                                                                                                 gCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATca                                      <  1:332578/150‑1 (MQ=255)
                                                                                                                 gCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATca                                      <  1:1502236/150‑1 (MQ=255)
                                                                                                                     gagTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGcc                                             <  2:1609913/139‑1 (MQ=255)
                                                                                                                     gagTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGcc                                             >  1:1609913/1‑139 (MQ=255)
                                                                                                                     gagTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGc                                  <  2:1273329/150‑1 (MQ=255)
                                                                                                                         gACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAAt                              >  2:121130/1‑150 (MQ=255)
                                                                                                                             aaCAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGCTGTGATTTTCGGGCCTAAATCACAGCCAATcac                           >  2:443322/1‑150 (MQ=255)
                                                                                                                               ccgtaCCGGCATTGTCATGACTGTGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACAtt                        >  1:921551/5‑150 (MQ=37)
                                                                                                                                 gCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTgg                      <  2:483824/150‑1 (MQ=255)
                                                                                                                                  cACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGc                                                                                             <  2:924650/79‑1 (MQ=255)
                                                                                                                                  cACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGc                                                                                             >  1:924650/1‑79 (MQ=255)
                                                                                                                                   aCCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCa                    <  1:175965/150‑1 (MQ=255)
                                                                                                                                    ccGGCATTGTCATGACTGTGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCAt                   >  1:953667/1‑150 (MQ=255)
                                                                                                                                    ccGGCATTGTCATGACTGTGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCAt                   >  1:439210/1‑150 (MQ=255)
                                                                                                                                    ccGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCAt                   <  2:282244/150‑1 (MQ=255)
                                                                                                                                     cGGCATTGTCATGACTGTGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGTGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATca                            <  2:1605013/140‑1 (MQ=37)
                                                                                                                                     cGGCATTGTCATGACTGTGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGGGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATca                            >  1:1605013/1‑140 (MQ=255)
                                                                                                                                      ggCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATGTGATTGGctgtg                 <  2:433395/150‑5 (MQ=255)
                                                                                                                                            gTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCAtt           >  1:1605738/1‑150 (MQ=255)
                                                                                                                                             tCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTg          >  2:868088/1‑150 (MQ=255)
                                                                                                                                                tGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCct       >  2:586016/1‑150 (MQ=255)
                                                                                                                                                  aCAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCctct     <  2:1155116/150‑1 (MQ=255)
                                                                                                                                                     aTGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTg  <  2:542757/150‑1 (MQ=255)
                                                                                                                                                     aTGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTg  <  1:295732/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GAAGCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTG  >  NC_000913/2772088‑2772386

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