Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,872,913 A→G 6.7% intergenic (‑93/+101) ftsB ← / ← ygbE cell division protein FtsB/conserved inner membrane protein YgbE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,872,9130AG 6.7% 223.2 / 5.4 75intergenic (‑93/+101)ftsB/ygbEcell division protein FtsB/conserved inner membrane protein YgbE
Reads supporting (aligned to +/- strand):  ref base A (33/37);  new base G (2/3);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

GCGAATACTGTAGCCAGACCAGAATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTAC  >  NC_000913/2872766‑2873058
                                                                                                                                                   |                                                                                                                                                 
gCGACTACTGTCGCCAGACCAGAATAGCCAGCACCAGCAGCGTTAGTTTCCCCATCCTGCCCCCTGAAAAACGGCATCATCATCCCCTGCATCCGACGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTAcc                                                                                                                                                 <  1:1537579/150‑1 (MQ=255)
gCGAATACTGTAGCCAGACCAGAATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTAcc                                                                                                                                                 >  1:7311/1‑150 (MQ=255)
  gAATACTGTAGCCAGACCAGAATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACcac                                                                                                                                               <  1:1242182/150‑1 (MQ=255)
        tGTAGCCAGACCAGAATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGt                                                                                                                                         >  2:27236/1‑150 (MQ=255)
           aGCCAGACCAGAATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCa                                                                                                                                      >  1:928677/1‑150 (MQ=255)
           aGCCAGACCAGAATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCa                                                                                                                                      >  1:1039494/1‑150 (MQ=255)
           aGCCAGACCAGAATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGACAAATGTGGTACATTTGcccg                                                                                                                                      <  1:1387294/150‑2 (MQ=255)
           aGCCAGACCAGAATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGACAAATGTGGTACATTTGcccg                                                                                                                                      <  1:1282974/150‑2 (MQ=255)
            gCCAGACCAGAATAGCCAGCAACAGCAGCGTTAGTTTACCAATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCAt                                                                                                                                     <  1:1608910/150‑1 (MQ=255)
             ccagaccagaATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGACAAATGTGGTACATTTGccc                                                                                                                                       <  1:1078735/147‑1 (MQ=255)
             ccagaccagaATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGACAAATGTGGTACATTTGccc                                                                                                                                       >  2:1078735/1‑147 (MQ=255)
              cagaccagaATAGCCAGCATCAGCAGCGTTAGTTTACCCATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTg                                                                                                                                   <  2:516889/150‑1 (MQ=255)
                 accagaATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTa                                                                                                                                <  1:1236352/150‑1 (MQ=255)
                 accagaATAGCCAGCAACACCAGCGTTAGTTTACCCATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTa                                                                                                                                <  1:1446623/150‑1 (MQ=255)
                         aGCCAGCAACAGCAGCGTTAGTTTACCCATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATAcc                                                                                                                        <  1:1182634/150‑1 (MQ=255)
                          gCCAGCAACAGCAGCGTTAGTTTACCCATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATAccc                                                                                                                       >  2:466470/1‑150 (MQ=255)
                          gCCAGCAACAGCAGCGTTAGTTTACCCATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATAccc                                                                                                                       >  2:308327/1‑150 (MQ=255)
                          gCCAGCAACAGCAGCGTTAGTTTACCCATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATAccc                                                                                                                       >  1:154209/1‑150 (MQ=255)
                           ccAGCAACAGCAGCGTTAGTTTACCCATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCa                                                                                                                      >  2:1118970/1‑150 (MQ=255)
                             aGCAACAGCAGCGTTAGTTTACCCATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAgg                                                                                                                    <  1:308327/150‑1 (MQ=255)
                               cAACAGCAGCGTTAGTTTACCCATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGc                                                                                                                  <  2:897538/150‑1 (MQ=255)
                                   cccaccGTTAGTTTACCCATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAACACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGc                                                                                                              <  1:712337/145‑1 (MQ=255)
                                       gcGTTAGTTTACCCATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGaa                                                                                                          >  2:707606/1‑150 (MQ=255)
                                          ttAGTTTACCCATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACAt                                                                                                       <  1:1249286/150‑1 (MQ=255)
                                            aGTTTACCCATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATaa                                                                                                     >  2:431588/1‑150 (MQ=255)
                                            aGTTTACCCATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATaa                                                                                                     <  2:459141/150‑1 (MQ=255)
                                                 aCCCATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCa                                                                                                <  1:1089260/150‑1 (MQ=255)
                                                  cccATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGCCAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCat                                                                                               >  2:297628/1‑150 (MQ=255)
                                                   ccATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCatt                                                                                              <  2:493430/150‑1 (MQ=255)
                                                    cATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCatta                                                                                             <  2:186848/150‑1 (MQ=255)
                                                     aTCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCattat                                                                                            <  1:509863/150‑1 (MQ=255)
                                                       ccTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCAttatta                                                                                          <  2:165457/150‑1 (MQ=255)
                                                       ccTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCAttatta                                                                                          >  1:485808/1‑150 (MQ=255)
                                                           cccccTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTa                                                                                      >  1:301103/1‑150 (MQ=255)
                                                               cTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGt                                                                                  >  2:1318741/1‑150 (MQ=255)
                                                               cTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGt                                                                                  >  2:806857/1‑150 (MQ=255)
                                                                   aaaaCGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGa                                                                              <  2:928677/150‑1 (MQ=255)
                                                                      aCGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAAtt                                                                           <  1:137237/150‑1 (MQ=255)
                                                                         gCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATg                                                                        <  2:285364/150‑1 (MQ=255)
                                                                         gCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATg                                                                        >  1:732386/1‑150 (MQ=255)
                                                                            tcatcatcCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAAc                                                                     <  1:806857/150‑1 (MQ=255)
                                                                            tcatcatcCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAAc                                                                     <  1:304559/150‑1 (MQ=255)
                                                                            tcatcatcCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAAc                                                                     >  2:761161/1‑150 (MQ=255)
                                                                            tcatcatcCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAAc                                                                     <  1:466470/150‑1 (MQ=255)
                                                                                catcCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACagag                                                                 >  2:381826/1‑150 (MQ=255)
                                                                                aatcCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACagag                                                                 <  2:178793/149‑1 (MQ=255)
                                                                                 atcCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAgagg                                                                <  2:1341528/150‑1 (MQ=255)
                                                                                   cccATGCCTCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAgaggag                                                              <  2:1568337/150‑1 (MQ=255)
                                                                                   cccATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAgaggag                                                              <  2:114128/150‑1 (MQ=255)
                                                                                   cccATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAgaggag                                                              <  2:1422090/150‑1 (MQ=255)
                                                                                        gCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACggg                                                         >  2:1023059/1‑150 (MQ=255)
                                                                                            ccGAAGACGCCTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAg                                                     >  1:66753/1‑150 (MQ=255)
                                                                                            ccGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACCTTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAg                                                     >  2:829986/1‑150 (MQ=255)
                                                                                            ccGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAg                                                     <  2:1039494/150‑1 (MQ=255)
                                                                                             cGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGt                                                    >  1:1537666/1‑150 (MQ=255)
                                                                                                 gacgacTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACaa                                                >  1:1480172/1‑150 (MQ=255)
                                                                                                 gacgacTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACaa                                                <  1:27236/150‑1 (MQ=255)
                                                                                                  acgacTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACaaa                                               <  1:91756/150‑1 (MQ=255)
                                                                                                        ctACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGcc                                         >  1:1112827/1‑150 (MQ=255)
                                                                                                           cATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACAt                                                                                                       <  2:1486626/85‑1 (MQ=255)
                                                                                                           cATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACAt                                                                                                       >  1:1486626/1‑85 (MQ=255)
                                                                                                           cATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAg                                      <  1:761161/150‑1 (MQ=255)
                                                                                                            aTCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGt                                     >  1:800575/1‑150 (MQ=255)
                                                                                                             tCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTa                                    <  1:1318741/150‑1 (MQ=255)
                                                                                                               ctctGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGc                                  >  2:645905/1‑150 (MQ=255)
                                                                                                                      ggggATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATaa                           >  1:708779/1‑150 (MQ=255)
                                                                                                                             ccGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTgc                    >  1:388741/1‑150 (MQ=255)
                                                                                                                                cgACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGcc                 >  1:403072/1‑150 (MQ=255)
                                                                                                                                    aaCGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAAc             <  1:645905/150‑1 (MQ=255)
                                                                                                                                     aCGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACa            >  2:276837/1‑150 (MQ=255)
                                                                                                                                       gcgGGCAAATGTACCACATTTGTCGATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATa          <  2:344866/150‑1 (MQ=255)
                                                                                                                                        tggacAAATGTGGTACATTTGCCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATaa         >  1:1110540/5‑150 (MQ=255)
                                                                                                                                        cgGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATaa         >  1:674147/1‑150 (MQ=255)
                                                                                                                                            cAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGcc     >  1:1213323/1‑150 (MQ=255)
                                                                                                                                               aTGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTAc  <  1:431588/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GCGAATACTGTAGCCAGACCAGAATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTAC  >  NC_000913/2872766‑2873058

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