Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 176,979 C→T 6.5% intergenic (+25/+22) erpA → / ← yadS iron‑sulfur cluster insertion protein/UPF0126 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913176,9790CT 6.5% 331.7 / 8.0 109intergenic (+25/+22)erpA/yadSiron‑sulfur cluster insertion protein/UPF0126 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (49/51);  new base T (3/4);  total (52/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

Only 100 of 115 total aligned reads displayed.

TGAGCCTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATGTGT  >  NC_000913/176830‑177128
                                                                                                                                                     |                                                                                                                                                     
tGAGCCTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAAc                                                                                                                                                       >  1:1056717/1‑150 (MQ=255)
tGAGCCTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAAc                                                                                                                                                       <  2:1523491/150‑1 (MQ=255)
     cTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATcc                                                                                                                                                  >  2:755658/1‑150 (MQ=255)
        cAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGc                                                                                                                                               <  2:1673374/150‑1 (MQ=255)
        cAATATCGGGGCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTACCTAATCTGTTGTTCATGTGCCGGAGAAACTATCCGGc                                                                                                                                               <  1:907291/150‑1 (MQ=255)
         aaTATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCCGca                                                                                                                                              >  1:1031059/1‑150 (MQ=255)
          atatCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAGTTTATCCGGcac                                                                                                                                             <  1:17997/150‑1 (MQ=255)
           tatCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAGTTTATCCGGcaca                                                                                                                                            <  2:1127090/150‑1 (MQ=255)
            atCTGGTCGTCGGTTCCGTTGATTATACCGAATGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGATCATGTGCCGGATAGTTTATCCGGCACAt                                                                                                                                           <  1:129454/150‑1 (MQ=255)
              cTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACAtct                                                                                                                                         >  1:1365394/1‑150 (MQ=255)
                 gTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACAtcttct                                                                                                                                      >  1:1618830/1‑150 (MQ=255)
                  tCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTtctc                                                                                                                                     >  1:267669/1‑150 (MQ=255)
                    ggcggTTCCGTTGATTATACCTAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCAGGCACATCTtctatc                                                                                                                                   >  2:541975/1‑150 (MQ=255)
                       ggTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTc                                                                                                                                <  2:1512031/150‑1 (MQ=255)
                        gTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCa                                                                                                                               <  1:1478562/150‑1 (MQ=255)
                               tGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGTCACATCTTCTCTCTTCAACGCCCa                                                                                                                        >  2:378799/1‑150 (MQ=255)
                                gATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCAt                                                                                                                       >  1:1712062/1‑150 (MQ=255)
                                 aTTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCAtt                                                                                                                      <  2:680939/150‑1 (MQ=255)
                                 aTTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCAtt                                                                                                                      >  1:771767/1‑150 (MQ=255)
                                   tataCCTAAGGTCTGGAAGGTTCTAGTTTCATCTTGACCAACCCTAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATATGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATtct                                                                                                                    <  1:541975/150‑1 (MQ=255)
                                   tataCCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATtct                                                                                                                    <  2:1711929/150‑1 (MQ=255)
                                    ataCCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATtctc                                                                                                                   <  2:713144/150‑1 (MQ=255)
                                    ataCCGAAGGTCTGGAACGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTCTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATtcta                                                                                                                   >  2:1110844/1‑149 (MQ=255)
                                      aCCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCAt                                                                                                                 <  2:905929/150‑1 (MQ=255)
                                       ccGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATc                                                                                                                >  1:554893/1‑150 (MQ=255)
                                        cGAAGTGCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCAGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATcc                                                                                                               <  2:45357/150‑1 (MQ=255)
                                           aGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAgc                                                                                                            <  1:158282/150‑1 (MQ=255)
                                           aGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAgc                                                                                                            <  1:1137202/150‑1 (MQ=255)
                                            ggTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAGACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAgcg                                                                                                           >  1:172375/1‑150 (MQ=255)
                                            ggTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAgcg                                                                                                           >  1:878358/1‑150 (MQ=255)
                                               cTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTACTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATACGGCCCATCTTCTCTCTTCAACGCCCATTCTCAACCAGCGaaa                                                                                                        >  1:884404/1‑150 (MQ=255)
                                                 ggAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTTTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAAc                                                                                                      >  2:660468/1‑150 (MQ=255)
                                                   atggTTGTCGTTTCATCGTTACCAACCTGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGt                                                                                                    <  1:207131/148‑1 (MQ=255)
                                                    tggTTCTCGTTTCATCGTTACCAACCCGAACGCGAAAACCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTc                                                                                                   <  2:1084179/149‑1 (MQ=255)
                                                      gTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTGTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACAACTTCTCTCTTCAACGCCGATTCTCATCCAGAGCTAACGTCgg                                                                                                 >  1:323363/1‑150 (MQ=255)
                                                         ttcGTTTCATCGTGACCAACCCGAACGCTAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAg                                                                                              <  2:627698/149‑1 (MQ=255)
                                                         ctcGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCGTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCCTCCAGCGCAAACGGCGGTAg                                                                                              >  2:1270470/1‑150 (MQ=255)
                                                          tcGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGc                                                                                             <  1:1005943/150‑1 (MQ=255)
                                                           cGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCt                                                                                            >  2:1238769/1‑150 (MQ=255)
                                                             tttCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTa                                                                                          <  1:1341071/150‑1 (MQ=255)
                                                                cATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGa                                                                                       >  1:55538/1‑150 (MQ=255)
                                                                 aTCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGAt                                                                                      <  2:1618830/150‑1 (MQ=255)
                                                                 aTCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGAt                                                                                      >  1:695585/1‑150 (MQ=255)
                                                                  tCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATg                                                                                     <  1:755658/150‑1 (MQ=255)
                                                                   cGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGc                                                                                    <  1:1209784/150‑1 (MQ=255)
                                                                    gTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGcc                                                                                   <  2:1219131/150‑1 (MQ=255)
                                                                    gTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTGGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGcc                                                                                   <  1:530653/150‑1 (MQ=255)
                                                                       aCCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAAc                                                                                >  2:1430875/1‑150 (MQ=255)
                                                                        ccAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACg                                                                               <  2:969336/150‑1 (MQ=255)
                                                                        ccAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACg                                                                               >  1:456425/1‑150 (MQ=255)
                                                                        ccAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACg                                                                               <  2:878358/150‑1 (MQ=255)
                                                                          aaCCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCATCTGTCTTCAACGCCCATTCTCATCCAGAGCAAACGACGGTCTCTTAAGATGCCAACGaa                                                                             >  1:877338/1‑150 (MQ=255)
                                                                           aCCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAAt                                                                            <  2:267669/150‑1 (MQ=255)
                                                                           aCCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAAt                                                                            <  2:554893/150‑1 (MQ=255)
                                                                             ccGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAGTTTATCCGGCACAtgaac                                                                                                                                       <  2:1664738/89‑5 (MQ=37)
                                                                             ccGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAGTTTATCCGGCACAtgaac                                                                                                                                       >  1:1664738/1‑85 (MQ=37)
                                                                                aaCGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAtccg                                                                       <  2:10841/150‑1 (MQ=255)
                                                                                 aCGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAgccgc                                                                      >  1:851070/1‑150 (MQ=255)
                                                                                  cgcgAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAgccgcc                                                                     <  1:1146155/150‑1 (MQ=255)
                                                                                     gAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCATCGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGACAACGAATAGCCGCCAGc                                                                  >  1:1059353/1‑150 (MQ=255)
                                                                                         aGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGaa                                                              >  1:934208/1‑150 (MQ=255)
                                                                                         aGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGaa                                                              <  2:1419359/150‑1 (MQ=255)
                                                                                             ccTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCaa                                                          >  2:351999/1‑150 (MQ=255)
                                                                                                 cggttgcggttCTGCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCAAATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGc                                                      <  1:1568697/150‑1 (MQ=255)
                                                                                                      gcggttCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGAAAACGTCGGTAGCTTAAGATGCCAACGAATAGACGCCAGCCGAATCAATAGCGTCAc                                                 >  2:525903/1‑150 (MQ=255)
                                                                                                        ggttCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGa                                               >  1:700056/1‑150 (MQ=255)
                                                                                                         gttCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGAc                                              <  1:980857/150‑1 (MQ=255)
                                                                                                           tcttcCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCa                                            >  1:1209422/1‑150 (MQ=255)
                                                                                                           tcttcCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCa                                            <  1:1298252/150‑1 (MQ=255)
                                                                                                            cttcCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCAt                                           <  2:514512/150‑1 (MQ=255)
                                                                                                                  ttAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCcat                                     <  2:1209422/150‑1 (MQ=255)
                                                                                                                   tAGTATCTAATCTGTTGTTAATGTGCAGGATAAAATATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCcatc                                    <  2:1200984/150‑1 (MQ=255)
                                                                                                                    aGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACCAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCcatca                                   >  1:1167272/1‑150 (MQ=255)
                                                                                                                     ttATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCCTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCcatcat                                  <  2:1279240/149‑1 (MQ=255)
                                                                                                                         cTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTg                              >  1:1432417/1‑150 (MQ=255)
                                                                                                                         cTAAACAGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCACTTCAACGCGCATTCTCCTCCAGCGCTAACGTCGGTCGATTAAGATGCCAACGAATTGCCGCAAGCCGAATCAAAAGCGGCCCGCCCAGCCCCAGCATACTg                              >  2:6660/1‑150 (MQ=255)
                                                                                                                             tCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTg                          <  1:34216/150‑1 (MQ=255)
                                                                                                                                gttgttCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCTTCATACTGGCTGttt                       >  2:450485/1‑150 (MQ=255)
                                                                                                                                gttgttCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGACCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGAAGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGttt                       >  2:1402245/1‑150 (MQ=255)
                                                                                                                                  tgttCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTcc                     <  2:1432417/150‑1 (MQ=255)
                                                                                                                                    ttCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAg                   <  1:660468/150‑1 (MQ=255)
                                                                                                                                    tgcATGTTCCGGATAAAATATCCGGCACATCTTCTCTCTACAACGCCCATTCTCATCCAGCGCAACCGTCGGTAGCGTAAGATGCCAACTAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAg                   <  2:1076562/148‑1 (MQ=255)
                                                                                                                                       aTGTGCCGGATAGTTTATCCGGCACATCTACTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCACCGAAAAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTgg                >  1:328658/1‑150 (MQ=255)
                                                                                                                                         gtgCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTa              <  1:1732280/150‑1 (MQ=255)
                                                                                                                                         gtgCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGGAGCATAAGCTGCCAACGAATAGCCGCCAGCCGAATCAAGAGCGGCACGACCATGGCCAGCATACTGGCTGTTTCCAGGGGTa              >  1:1376435/1‑150 (MQ=255)
                                                                                                                                          tgCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTAc             >  2:893115/1‑150 (MQ=255)
                                                                                                                                           gCCGGATAGTTTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGCCCATGCCCATCATACTGGCTGTTTCCAGTGGTACg            >  2:1654392/1‑150 (MQ=255)
                                                                                                                                           gCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCCTCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACg            >  1:39562/1‑150 (MQ=255)
                                                                                                                                           gCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCGTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACg            >  1:251114/1‑150 (MQ=255)
                                                                                                                                           gCCGGATAAACTATCCGGCACATCTTATCTGTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAAGAGCCGCCAGCCGAATCAATAGCGTGACGACCATGGCCATCATACTGGCTGTTTCCAGTGGTACg            >  2:135096/1‑150 (MQ=255)
                                                                                                                                              ggATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAAAACCCGACAGCCGAATCAAAAGCGTCACCACCATGCCCATCATACTGGCTGATTCCACTGGTACGGaa         >  1:297814/1‑150 (MQ=255)
                                                                                                                                                aTAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGaaaa       >  1:562825/1‑150 (MQ=255)
                                                                                                                                                 tAAACTATCCGGAACATCTTCTCTCTTCAACGCCCATTCTCATCCCGCGCAAACGGCGGGAGCTTAAGATCCCATCGAATAGCCTCCAGCCGAATCAATAGCGGCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAAt      <  2:334444/150‑1 (MQ=255)
                                                                                                                                                   aaCTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACAATGCCAATCATACTGGCTTTTTCCAGTGGTACGGAAAAtgt    >  1:1593858/1‑150 (MQ=255)
                                                                                                                                                    aCTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGATTAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAAtgtg   <  2:599537/150‑1 (MQ=255)
                                                                                                                                                    aCTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAAtgtg   <  2:1058438/150‑1 (MQ=255)
                                                                                                                                                    aCTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAAtgtg   <  1:663635/150‑1 (MQ=255)
                                                                                                                                                    aCTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAAtgtg   <  2:695585/150‑1 (MQ=255)
                                                                                                                                                     cTATCCGGCTCATCTTCTCTCTTCAACGCGCATTCTCATCCAGCGCAAACGTCTGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGTAAAATgtgt  <  1:701015/150‑1 (MQ=255)
                                                                                                                                                     cTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATgtgt  <  1:1430875/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGAGCCTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATGTGT  >  NC_000913/176830‑177128

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