Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,872,915 C→T 6.9% intergenic (‑95/+99) 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,9150CT 6.9% 223.9 / 4.7 72intergenic (‑95/+99)ftsB/ygbEcell division protein FtsB/conserved inner membrane protein YgbE
Reads supporting (aligned to +/- strand):  ref base C (32/35);  new base T (2/3);  total (34/38)
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.78e-01

GCGAATACTGTAGCCAGACCAGAATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCA  >  NC_000913/2872766‑2873063
                                                                                                                                                     |                                                                                                                                                    
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)
                                                                                                                                                    ccACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCa  >  2:623136/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCGAATACTGTAGCCAGACCAGAATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACATCCTCTGTTGGGGATACCGCGACAACGCGGGCAAATGTACCACATTTGTCCATTGTTACGTATACCCAGGGCGTGCAGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCA  >  NC_000913/2872766‑2873063

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: