Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,249,909 +C intergenic (+47/‑176) yjbI → / → ubiC hypothetical protein/chorismate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,249,9051.C99.6% 753.8 / ‑4.2 235intergenic (+43/‑180)yjbI/ubiChypothetical protein/chorismate lyase
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (121/113);  total (122/113)
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.96e-01

Only 100 of 253 total aligned reads displayed.

TTTTTACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATG‑CCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTT  >  NC_000913/4249813‑4249991
                                                                                             |                                                                                      
tttttACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATg                                                                                         >  1:239999/1‑93 (MQ=255)
 ttttACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATG‑c                                                                                       >  1:49193/1‑93 (MQ=255)
  tttACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATG‑cc                                                                                      >  2:5855598/1‑93 (MQ=255)
     aCTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATG‑ccccc                                                                                   <  2:1969559/93‑2 (MQ=255)
     aCTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATG‑ccccc                                                                                   >  2:170523/1‑92 (MQ=255)
          gAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACgg                                                                               <  2:2271303/93‑1 (MQ=255)
           aGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGt                                                                              >  1:6270445/1‑93 (MQ=255)
              tGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGc                                                                           <  1:3807327/93‑1 (MQ=255)
              tGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCAGCCTGGACGGTGATGCCCCCTACGGTTGc                                                                           <  1:3868686/93‑1 (MQ=255)
                  gAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTc                                                                       <  1:742130/93‑1 (MQ=255)
                  gAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTc                                                                       >  1:3487891/1‑93 (MQ=255)
                   aGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCg                                                                      >  1:3192697/1‑93 (MQ=255)
                    gCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGc                                                                     <  2:2474637/93‑1 (MQ=255)
                    gCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGc                                                                     >  1:5525578/1‑93 (MQ=255)
                     cTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGcc                                                                    >  2:5755704/1‑93 (MQ=255)
                       cATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAg                                                                  >  2:296103/1‑93 (MQ=255)
                       cATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGACGCCCCCTACGGTTGCCCTCGCCAg                                                                  <  2:6074265/93‑1 (MQ=255)
                        aTCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCCCCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGc                                                                 <  1:4488800/93‑1 (MQ=255)
                         tCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCa                                                                <  2:1039190/93‑1 (MQ=255)
                         tCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCa                                                                >  2:218147/1‑93 (MQ=255)
                         tCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCa                                                                <  1:3022632/93‑1 (MQ=255)
                         tCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTAGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCACGCCAGCa                                                                <  1:4064536/93‑1 (MQ=255)
                           gTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACg                                                              >  1:1935559/1‑93 (MQ=255)
                            tataATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACgg                                                             >  1:1988103/1‑93 (MQ=255)
                            tataATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACgg                                                             <  2:5623862/93‑1 (MQ=255)
                             ataATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACggg                                                            <  2:3351256/93‑1 (MQ=255)
                             ataATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACggg                                                            <  2:4548882/93‑1 (MQ=255)
                              taATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGc                                                           >  2:4218452/1‑93 (MQ=255)
                                 tGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATc                                                        >  2:4965150/1‑93 (MQ=255)
                                  gAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCg                                                       <  2:4635379/93‑1 (MQ=255)
                                    ttAAGGGTATTTTAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTCCGGTTGCCCTCGCCAGCACGGGCATCGGt                                                     <  2:6300124/92‑1 (MQ=255)
                                     tAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTa                                                    >  1:5901960/1‑93 (MQ=255)
                                        gggTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAg                                                 <  1:131299/93‑1 (MQ=255)
                                        gggTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAg                                                 >  1:5371728/1‑93 (MQ=255)
                                        gggTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAg                                                 >  1:1787591/1‑93 (MQ=255)
                                        gggTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAg                                                 <  2:5990155/93‑1 (MQ=255)
                                        gggTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAg                                                 >  1:3263093/1‑93 (MQ=255)
                                            gtTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTa                                             >  2:5633064/1‑93 (MQ=255)
                                             tttAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTaa                                            >  1:5057067/1‑93 (MQ=255)
                                              ttAAGTAAAGGAAAACATCCCCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAg                                           <  1:1526050/93‑1 (MQ=255)
                                              ttAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAg                                           <  2:3778449/93‑1 (MQ=255)
                                              ttAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAg                                           >  2:3253388/1‑93 (MQ=255)
                                                 aGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGtt                                        >  1:4328485/1‑93 (MQ=255)
                                                 aGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGtt                                        <  2:3313909/93‑1 (MQ=255)
                                                 aGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGtt                                        <  1:776521/93‑1 (MQ=255)
                                                    aaaGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTCCGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCaa                                     <  1:4087810/93‑1 (MQ=255)
                                                    aaaGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCaa                                     <  1:5623706/93‑1 (MQ=255)
                                                    aaaGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCaa                                     >  2:4264795/1‑93 (MQ=255)
                                                    aaaGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCaa                                     >  1:4681762/1‑93 (MQ=255)
                                                      aGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACa                                   <  1:419148/93‑1 (MQ=255)
                                                       ggAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACAt                                  <  1:1023694/93‑1 (MQ=255)
                                                        gAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCaagctc                                 >  2:2105643/1‑89 (MQ=255)
                                                        gAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATc                                 <  1:3974349/93‑1 (MQ=255)
                                                        gAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATc                                 <  2:260280/93‑1 (MQ=255)
                                                         aaaaCATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCg                                >  2:4902787/1‑93 (MQ=255)
                                                         aaaaCATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCg                                >  2:5440776/1‑93 (MQ=255)
                                                          aaaCATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGt                               >  2:1444784/1‑93 (MQ=255)
                                                          aaaCATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGt                               <  2:1838063/93‑1 (MQ=255)
                                                          aaaCATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGt                               >  2:3914775/1‑93 (MQ=255)
                                                           aaCATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGtt                              >  2:3296386/1‑93 (MQ=255)
                                                           aaCATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGtt                              >  1:2886828/1‑93 (MQ=255)
                                                            aCATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGttt                             >  1:1196821/1‑93 (MQ=255)
                                                             cATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTCCAACATCGtgtt                            >  1:1872621/1‑93 (MQ=255)
                                                               tCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTAc                          >  2:3507557/1‑93 (MQ=255)
                                                               tCACCGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTAc                          >  2:3305627/1‑93 (MQ=255)
                                                                  ccGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCAc                       >  2:3166665/1‑93 (MQ=255)
                                                                   cGTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGCTCAACATCGTTTTACCACt                      >  2:1074598/1‑93 (MQ=255)
                                                                   cGTTCCTGGCATCCTGGACGGAGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGGTCGACATCGTTTTAcgagt                      >  1:1306981/1‑89 (MQ=255)
                                                                    gTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACtt                     >  1:5273880/1‑93 (MQ=255)
                                                                    gTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACtt                     >  1:3700036/1‑93 (MQ=255)
                                                                    gTTCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACtt                     <  1:3465627/93‑1 (MQ=255)
                                                                     ttCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTc                    >  2:3025546/1‑93 (MQ=255)
                                                                     ttCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTc                    >  1:559570/1‑93 (MQ=255)
                                                                      tCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCa                   <  1:66872/93‑1 (MQ=255)
                                                                      tCCTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCa                   <  2:5372002/93‑1 (MQ=255)
                                                                        cTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATg                 >  2:5251714/1‑93 (MQ=255)
                                                                        cTGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGCCATCGGTAAAGCGTAAGGTTCAACATCGGTTTACCACTTCAtc                 >  1:3039612/1‑92 (MQ=255)
                                                                         tGGCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGc                >  2:5144962/1‑93 (MQ=255)
                                                                         tGGCAACCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGc                >  2:4107342/1‑93 (MQ=255)
                                                                         cggCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGc                <  1:5856945/92‑1 (MQ=255)
                                                                          ggCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCg               <  2:3089815/93‑1 (MQ=255)
                                                                           gCATCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGa              >  2:2814485/1‑93 (MQ=255)
                                                                             aTCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGAtt            >  2:2124536/1‑93 (MQ=255)
                                                                             aTCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGAtt            >  2:1543875/1‑93 (MQ=255)
                                                                             aTCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGAtt            >  2:3416646/1‑93 (MQ=255)
                                                                              tCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGAttg           >  2:5774231/1‑93 (MQ=255)
                                                                              tCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGAttg           >  1:1042600/1‑93 (MQ=255)
                                                                              tCCTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGAttg           >  1:217270/1‑93 (MQ=255)
                                                                               ccTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGAttgt          <  2:2917137/93‑1 (MQ=255)
                                                                               ccTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGAttgt          <  1:252684/93‑1 (MQ=255)
                                                                                cTGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGAttgtt         >  1:2105682/1‑93 (MQ=255)
                                                                                 tGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGAttgttg        >  2:3525464/1‑93 (MQ=255)
                                                                                 tGGACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGAttgttg        <  2:2537987/93‑1 (MQ=255)
                                                                                  ggACGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGc       <  2:4478054/93‑1 (MQ=255)
                                                                                    aCGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGt     <  1:5123845/93‑1 (MQ=255)
                                                                                    aCGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGt     <  1:2722645/93‑1 (MQ=255)
                                                                                    aCGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGt     <  2:850767/93‑1 (MQ=255)
                                                                                     cGGTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGtt    <  1:1638238/93‑1 (MQ=255)
                                                                                       gTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGtttt  >  2:4061786/1‑93 (MQ=255)
                                                                                       gTGATGCCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGtttt  <  1:4831419/93‑1 (MQ=255)
                                                                                             |                                                                                      
TTTTTACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATG‑CCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTT  >  NC_000913/4249813‑4249991

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