Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,425,975 A→G 5.8% intergenic (‑59/+31) hisQ ← / ← hisJ histidine ABC transporter permease/histidine ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,425,9750AG 5.8% 447.2 / 12.8 138intergenic (‑59/+31)hisQ/hisJhistidine ABC transporter permease/histidine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (68/62);  new base G (5/3);  total (73/65)
Fisher's exact test for biased strand distribution p-value = 7.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.10e-02

Only 100 of 139 total aligned reads displayed.

GATTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATCTTCTT  >  NC_000913/2425830‑2426121
                                                                                                                                                 |                                                                                                                                                  
gatTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAgg                                                                                                                                                 <  2:169821/149‑1 (MQ=255)
  ttACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTTCGGGTAACCGGAGGGa                                                                                                                                               <  1:578644/149‑1 (MQ=255)
  ttACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGa                                                                                                                                               >  2:62797/1‑149 (MQ=255)
       gCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAACCCCATACAACGTGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCggg                                                                                                                                          <  1:1025712/149‑1 (MQ=255)
       gCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCggg                                                                                                                                          <  1:717708/149‑1 (MQ=255)
             aCTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGctttc                                                                                                                                     >  4:760021/1‑148 (MQ=255)
                 cAGAGCTGATAGCCTGCTCCAGCGTGACGAGCGCACCCTGGAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCtttctttc                                                                                                                                 >  4:781279/1‑148 (MQ=255)
                   gagCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCAc                                                                                                                               >  3:224369/1‑148 (MQ=255)
                         aTAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATgg                                                                                                                         >  3:526269/1‑148 (MQ=255)
                           aGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGa                                                                                                                      <  1:278645/149‑1 (MQ=255)
                            gCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGa                                                                                                                      <  4:526269/148‑1 (MQ=255)
                            gCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGAt                                                                                                                     <  2:1014063/149‑1 (MQ=255)
                               aGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAg                                                                                                                  >  1:625314/1‑149 (MQ=255)
                               aGCTCCAGCATGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTCCACCTATGGCGAtttg                                                                                                                  >  2:971507/1‑147 (MQ=255)
                                gCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAg                                                                                                                  <  4:137976/148‑1 (MQ=255)
                                gCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGc                                                                                                                 >  1:955817/1‑149 (MQ=255)
                                gCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGc                                                                                                                 >  2:532565/1‑149 (MQ=255)
                                 cTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCtttctttc                                                                                                                                 >  1:509643/1‑132 (MQ=255)
                                  tCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAgg                                                                                                                                                 >  2:807198/1‑115 (MQ=255)
                                  tCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGcca                                                                                                               <  2:593629/149‑1 (MQ=255)
                                  tCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGcca                                                                                                               <  1:388797/149‑1 (MQ=255)
                                  tCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGcca                                                                                                               <  1:284530/149‑1 (MQ=255)
                                  tCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGcc                                                                                                                >  4:424080/1‑148 (MQ=255)
                                     aGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGccacca                                                                                                            <  1:794621/149‑1 (MQ=255)
                                     aGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGccacc                                                                                                             >  4:32870/1‑148 (MQ=255)
                                           aCGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAAc                                                                                                       <  4:85520/148‑1 (MQ=255)
                                           aCGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACa                                                                                                      <  1:539861/149‑1 (MQ=255)
                                           aCGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACa                                                                                                      <  1:1023709/149‑1 (MQ=255)
                                                   aCCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATc                                                                                              >  1:385047/1‑149 (MQ=255)
                                                     ccTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCCCCGATGGCGATTAGCCACCATAAACATCAAAATCg                                                                                             <  4:123155/148‑1 (MQ=255)
                                                     ccTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGa                                                                                            >  1:963326/1‑149 (MQ=255)
                                                      cTGTAAAATAACACCTGAAAACCCATACAACCTGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCATGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGaa                                                                                           <  1:368963/149‑1 (MQ=255)
                                                      cTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGa                                                                                            <  3:781279/148‑1 (MQ=255)
                                                       tGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAg                                                                                          >  2:420606/1‑149 (MQ=255)
                                                       tGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAg                                                                                          >  1:216402/1‑149 (MQ=255)
                                                       tGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAg                                                                                          >  2:580876/1‑149 (MQ=255)
                                                        gTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAg                                                                                          >  3:37806/1‑148 (MQ=255)
                                                        gTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAg                                                                                          >  3:109488/1‑148 (MQ=255)
                                                        gTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAg                                                                                          >  3:1046614/1‑148 (MQ=255)
                                                        gTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGt                                                                                         >  2:726064/1‑149 (MQ=255)
                                                           aaaTAACACCTGAAAACCCATTCCACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTAc                                                                                       <  3:762290/148‑1 (MQ=255)
                                                            aaTAACACCTGAAAACCCATACAACATGTTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTGTTGTCTCCCTCCGGTTGCGCGAAGGCAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACtt                                                                                     >  2:768608/1‑149 (MQ=255)
                                                                     tGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTaa                                                                            >  2:865597/1‑149 (MQ=255)
                                                                     tGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTaa                                                                            <  2:622710/149‑1 (MQ=255)
                                                                     tGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTa                                                                             <  3:824570/148‑1 (MQ=255)
                                                                      gAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTaa                                                                            >  4:374751/1‑148 (MQ=255)
                                                                       aaaaCCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAAtt                                                                          >  1:82527/1‑149 (MQ=255)
                                                                       aaaaCCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAAt                                                                           >  4:975431/1‑148 (MQ=255)
                                                                        aaaCCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAAttt                                                                         <  1:420606/149‑1 (MQ=255)
                                                                           cccATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTtct                                                                       <  3:30113/148‑1 (MQ=255)
                                                                              aTACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTa                                                                   >  1:811796/1‑149 (MQ=255)
                                                                                  aacaTGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGt                                                                <  4:719353/148‑1 (MQ=255)
                                                                                     aTGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGaa                                                                                           >  2:516990/1‑118 (MQ=255)
                                                                                     aTGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGaa                                                                                           <  1:516990/118‑1 (MQ=255)
                                                                                     aTGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTAcc                                                             >  3:577715/1‑148 (MQ=255)
                                                                                          gcCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCggg                                                                                                                                          >  1:782986/1‑66 (MQ=255)
                                                                                          gcCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCggg                                                                                                                                          <  2:782986/66‑1 (MQ=255)
                                                                                           cctgtcctgtCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTAcc                                                             <  1:497274/142‑1 (MQ=255)
                                                                                           cctgtcctgtCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTAcc                                                             >  2:497274/1‑142 (MQ=255)
                                                                                            ctgtcctgtcGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAgcg                                                      >  4:88628/1‑148 (MQ=255)
                                                                                              gtcctgtcGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAgcgcg                                                    <  4:698840/148‑1 (MQ=255)
                                                                                              gtcctgtcGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAgcgcg                                                    <  4:1048080/148‑1 (MQ=255)
                                                                                              gtcctgtcGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAgcgcg                                                    >  3:55518/1‑148 (MQ=255)
                                                                                              gtcctgtcGGGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAGACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAgcgcg                                                    >  3:245775/1‑148 (MQ=255)
                                                                                               tcctgtcGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCa                                                  >  2:944249/1‑149 (MQ=255)
                                                                                               tcctgtcGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCa                                                  >  2:1044186/1‑149 (MQ=255)
                                                                                                     cGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCg                                             >  4:48025/1‑148 (MQ=255)
                                                                                                       tgtgtgGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGctttc                                                                                                                                     <  2:899911/58‑1 (MQ=255)
                                                                                                       tgtgtgGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGctttc                                                                                                                                     >  1:899911/1‑58 (MQ=255)
                                                                                                            ggtgAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAgg                                      >  4:239344/1‑148 (MQ=255)
                                                                                                             gtgAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCt                                    <  2:666011/149‑1 (MQ=255)
                                                                                                              tgAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGa                                                                                            <  4:94133/92‑1 (MQ=255)
                                                                                                              tgAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGa                                                                                            >  3:94133/1‑92 (MQ=255)
                                                                                                              tgAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCt                                    >  3:1622/1‑148 (MQ=255)
                                                                                                               gAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCATCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCtt                                   <  3:417134/148‑1 (MQ=255)
                                                                                                               gAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCtt                                   >  3:616524/1‑148 (MQ=255)
                                                                                                               gAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCtt                                   <  4:921765/148‑1 (MQ=255)
                                                                                                                  ttACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATCTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCt                                    >  2:236249/1‑144 (MQ=255)
                                                                                                                  ttACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATCTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCt                                    <  1:236249/144‑1 (MQ=255)
                                                                                                                     ccGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAg                            >  2:463274/1‑149 (MQ=255)
                                                                                                                      cGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGt                           <  1:944249/149‑1 (MQ=255)
                                                                                                                           ctTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCtt                       <  3:88628/148‑1 (MQ=255)
                                                                                                                           ctTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCtt                       <  4:646064/148‑1 (MQ=255)
                                                                                                                                ctcccttcGGGTAACCGGAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTa                                                                   <  3:6475/99‑1 (MQ=255)
                                                                                                                                 tcccttcGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAg                >  1:1064697/1‑149 (MQ=255)
                                                                                                                                 tcccttcGGGTAACCGGAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTtctc                                                                      >  1:715106/1‑95 (MQ=255)
                                                                                                                                  ccctccGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAgc                                                       <  2:248580/109‑1 (MQ=255)
                                                                                                                                   cctccGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGtt              <  2:980730/149‑1 (MQ=255)
                                                                                                                                    ctccGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTc             >  1:877859/1‑149 (MQ=255)
                                                                                                                                    ctccGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTc             >  1:78010/1‑149 (MQ=255)
                                                                                                                                     tccGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTc             >  1:230250/1‑148 (MQ=255)
                                                                                                                                      ccGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCg            >  4:4693/1‑148 (MQ=255)
                                                                                                                                      ccGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGt           <  1:981304/149‑1 (MQ=255)
                                                                                                                                       cGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGtt          <  2:749490/149‑1 (MQ=255)
                                                                                                                                             cccGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTAtctt     <  4:109488/148‑1 (MQ=255)
                                                                                                                                              ccGAAGGGAGCGGGCTTTCTTTCACCGATGGCGGTTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGATAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTAtcttct   >  2:1063957/1‑149 (MQ=255)
                                                                                                                                              ccGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAg                                       >  2:333627/1‑113 (MQ=255)
                                                                                                                                              ccGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTAtcttct   >  1:432080/1‑149 (MQ=255)
                                                                                                                                               cGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCg            >  4:615073/1‑139 (MQ=255)
                                                                                                                                               cGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTAtcttct   >  4:777878/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GATTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATCTTCTT  >  NC_000913/2425830‑2426121

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: