Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,425,968 T→G 7.0% intergenic (‑52/+38) 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,9680TG 7.0% 595.6 / 29.9 167intergenic (‑52/+38)hisQ/hisJhistidine ABC transporter permease/histidine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (76/79);  new base G (7/5);  total (83/84)
Fisher's exact test for biased strand distribution p-value = 5.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.35e-01

Only 100 of 168 total aligned reads displayed.

TTAAACCGATGATTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTC  >  NC_000913/2425820‑2426110
                                                                                                                                                    |                                                                                                                                              
ttAAACCGATGATTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGt                                                                                                                                                >  1:451528/1‑149 (MQ=255)
    aCCGATGATTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTAcc                                                                                                                                            <  2:150871/149‑1 (MQ=255)
     ccGATGATTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTAccc                                                                                                                                           >  4:597298/1‑149 (MQ=255)
        atgatTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGaa                                                                                                                                        <  4:487438/149‑1 (MQ=255)
            ttACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAgg                                                                                                                                      >  4:1859735/1‑147 (MQ=255)
            ttACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGa                                                                                                                                    <  1:29672/149‑1 (MQ=255)
            ttACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGa                                                                                                                                    >  2:1450365/1‑149 (MQ=255)
                aGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCgg                                                                                                                                <  1:1155608/149‑1 (MQ=255)
                 gCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGaa                                                                                                                                        <  2:1791302/140‑1 (MQ=255)
                 gCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGaa                                                                                                                                        >  1:657113/1‑140 (MQ=255)
                 gCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGaa                                                                                                                                        <  2:657113/140‑1 (MQ=255)
                 gCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGaa                                                                                                                                        >  1:1791302/1‑140 (MQ=255)
                         tACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAg                                                                                                                                       >  4:622590/1‑133 (MQ=255)
                         tACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGctttcttt                                                                                                                       >  2:1473083/1‑149 (MQ=255)
                         tACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGctttcttt                                                                                                                       >  1:1188809/1‑149 (MQ=255)
                            agagCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGctttc                                                                                                                          <  3:1682734/143‑1 (MQ=255)
                            agagCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGctttc                                                                                                                          >  4:1682734/1‑143 (MQ=255)
                             gagCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCAcc                                                                                                                   <  2:451528/149‑1 (MQ=255)
                              agCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGaa                                                                                                                                        <  1:1782695/127‑1 (MQ=255)
                              agCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGaa                                                                                                                                        >  2:1782695/1‑127 (MQ=255)
                               gCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTTCGGGTAACCGGAGGGa                                                                                                                                    >  4:1294696/1‑130 (MQ=255)
                                      gCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGAt                                                                                                          >  1:1148207/1‑149 (MQ=255)
                                       ccAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCa                                                                                                                     >  3:336815/1‑137 (MQ=255)
                                       ccAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGAtt                                                                                                         <  1:490156/149‑1 (MQ=255)
                                          gCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGc                                                                                                      <  1:939671/149‑1 (MQ=255)
                                            tCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGcca                                                                                                    <  4:588409/149‑1 (MQ=255)
                                              cAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCATGTGTGGTGAATTACCGCCTCTTGTCTCCCCTCGGGTAACCGGAGGGa                                                                                                                                    <  4:1187878/115‑1 (MQ=255)
                                              cAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCATGTGTGGTGAATTACCGCCTCTTGTCTCCCCTCGGGTAACCGGAGGGa                                                                                                                                    >  3:1187878/1‑115 (MQ=255)
                                                gCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGctttcttt                                                                                                                       <  4:1675252/126‑1 (MQ=255)
                                                gCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGctttcttt                                                                                                                       >  3:1675252/1‑126 (MQ=255)
                                                   tGACGAGCGCACCCTGTAAAATAACACCTGAATACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATaaa                                                                                             <  3:597298/149‑1 (MQ=255)
                                                   tGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATaaa                                                                                             <  2:1148207/149‑1 (MQ=255)
                                                            cACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTACTTTCACCGATGGCGATTAGCCACCATAAACATCAAAAt                                                                                    >  1:1465278/1‑149 (MQ=255)
                                                             aCCCTGTAAAATAACACCTGAAAACCGATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGctttc                                                                                                                          <  1:544598/110‑1 (MQ=255)
                                                             aCCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGctttc                                                                                                                          <  3:319458/110‑1 (MQ=255)
                                                             aCCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGctttc                                                                                                                          >  2:544598/1‑110 (MQ=255)
                                                             aCCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGctttc                                                                                                                          >  4:319458/1‑110 (MQ=255)
                                                             aCCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTa                                                                                                        >  1:1042399/1‑128 (MQ=255)
                                                             aCCCTGTAAAATAACACCGGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTa                                                                                                        <  2:1042399/128‑1 (MQ=255)
                                                                cTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGa                                                                                                           <  4:517678/122‑1 (MQ=255)
                                                                 tGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTTACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAg                                                                               >  3:1132421/1‑149 (MQ=255)
                                                                 tGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGaa                                                                                >  2:1356957/1‑148 (MQ=255)
                                                                  gTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGaa                                                                                <  2:939377/147‑1 (MQ=255)
                                                                  gTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGaa                                                                                >  1:939377/1‑147 (MQ=255)
                                                                   tAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGcc                                                                                                     >  4:1454885/1‑125 (MQ=255)
                                                                   tAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGcc                                                                                                     <  3:1454885/125‑1 (MQ=255)
                                                                   tAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTa                                                                             <  3:1758906/149‑1 (MQ=255)
                                                                      aaTAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACtt                                                                          >  3:1271709/1‑149 (MQ=255)
                                                                       aTAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGaa                                                                                >  3:383625/1‑142 (MQ=255)
                                                                       aTAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACttt                                                                         >  1:1591454/1‑149 (MQ=255)
                                                                          acacCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATg                                                                                                               <  1:1543287/108‑1 (MQ=255)
                                                                          acacCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATg                                                                                                               >  2:1543287/1‑108 (MQ=255)
                                                                             ccTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCt                                                                   >  3:499457/1‑149 (MQ=255)
                                                                              cTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGctttct                                                                                                                         <  4:1166164/94‑1 (MQ=255)
                                                                              cTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGctttct                                                                                                                         >  3:1166164/1‑94 (MQ=255)
                                                                                gAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAAt                                                                >  4:1769906/1‑149 (MQ=255)
                                                                                gAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAAt                                                                >  2:658465/1‑149 (MQ=255)
                                                                                     cccATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTtctc                                                           <  4:1249079/149‑1 (MQ=255)
                                                                                      ccATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGa                                                                                                           <  4:376504/100‑1 (MQ=255)
                                                                                      ccATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCg                                                          >  2:79493/1‑149 (MQ=255)
                                                                                      ccATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCg                                                          >  4:1842168/1‑149 (MQ=255)
                                                                                       cATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCCCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGaa                                                                                >  3:230983/1‑126 (MQ=255)
                                                                                        aTACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTa                                                                                                        >  1:442776/1‑101 (MQ=255)
                                                                                        aTACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTa                                                                                                        <  2:442776/101‑1 (MQ=255)
                                                                                         tACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCtttctttc                                                                                                                      >  3:1888304/1‑86 (MQ=255)
                                                                                         tACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTaa                                                       <  3:1842168/149‑1 (MQ=255)
                                                                                         tACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTaa                                                       <  1:79493/149‑1 (MQ=255)
                                                                                         tACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTaa                                                       <  4:1132421/149‑1 (MQ=255)
                                                                                          acaacaTGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACttttt                                                                       >  1:132024/1‑132 (MQ=255)
                                                                                            aacaTGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGAt                                                                                                                >  1:1439317/1‑89 (MQ=255)
                                                                                            aacaTGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGAt                                                                                                                <  2:1439317/89‑1 (MQ=255)
                                                                                              caTGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGt                                                         >  2:1149980/1‑142 (MQ=255)
                                                                                                 gctgcCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCaaa                                                                                      >  2:1306168/1‑110 (MQ=255)
                                                                                                   tgcCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTa                                                    <  1:1503070/142‑1 (MQ=255)
                                                                                                     cctgtcctgtCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTtctc                                                           <  4:1613488/133‑1 (MQ=255)
                                                                                                     cctgtcctgtCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTtctc                                                           >  3:1613488/1‑133 (MQ=255)
                                                                                                      ctgtcctgtcGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCATCATAAACATCAAAATCGAAGTACTTTTTCGc                                                                    <  4:229486/123‑1 (MQ=255)
                                                                                                          cctgtcGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCaa                                      <  2:1591454/149‑2 (MQ=255)
                                                                                                          cctgtcGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCAt                                      <  1:1257363/149‑1 (MQ=255)
                                                                                                                 tgtgtgGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGaa                                                                                >  3:528460/1‑100 (MQ=255)
                                                                                                                    gtggtgAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGt                                                >  3:619941/1‑129 (MQ=255)
                                                                                                                     tggtgAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAgg                           >  2:1218578/1‑149 (MQ=255)
                                                                                                                      ggtgAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGc                          >  4:425691/1‑149 (MQ=255)
                                                                                                                        tgAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCtt                        <  3:1359186/149‑1 (MQ=255)
                                                                                                                         gAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCttt                       <  1:1218578/149‑1 (MQ=255)
                                                                                                                          aaTTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTg                      >  3:659251/1‑149 (MQ=255)
                                                                                                                            ttACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTtct                                                            <  1:1274996/109‑1 (MQ=255)
                                                                                                                            ttACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTtct                                                            >  2:1274996/1‑109 (MQ=255)
                                                                                                                               ccGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAg                 >  1:1409696/1‑149 (MQ=255)
                                                                                                                                     ctTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTc           <  1:1006886/149‑1 (MQ=255)
                                                                                                                                        gTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTcgcg        >  3:1438463/1‑149 (MQ=255)
                                                                                                                                         gctcccttcGGGTAACCGGAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCgcgc       >  2:898953/2‑149 (MQ=255)
                                                                                                                                          ctcccttcGGGTAACCGGAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGaa                                                                                >  4:691654/1‑75 (MQ=37)
                                                                                                                                          ctcccttcGGGTAACCGGAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGaa                                                                                <  3:691654/75‑1 (MQ=37)
                                                                                                                                          ctccctccGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCa      >  1:564177/1‑149 (MQ=255)
                                                                                                                                           tcccttcGGGTAACCGGAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAg                 <  4:1406771/137‑1 (MQ=255)
                                                                                                                                           tcccttcGGGTAACCGGAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAg     >  2:337249/1‑149 (MQ=255)
                                                                                                                                           tccctccGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCttt                       <  4:328647/131‑1 (MQ=255)
                                                                                                                                            cccttcGGGTAACCGGAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCttt                       >  3:1775046/1‑130 (MQ=255)
                                                                                                                                              ctccGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTc  >  3:824641/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
TTAAACCGATGATTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTC  >  NC_000913/2425820‑2426110

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