Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 72,205 T→G 8.8% intergenic (+90/+24) yabI → / ← thiQ DedA family inner membrane protein/thiamine/thiamine pyrophosphate ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091372,2050TG 8.8% 322.0 / 15.2 111intergenic (+90/+24)yabI/thiQDedA family inner membrane protein/thiamine/thiamine pyrophosphate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (56/45);  new base G (4/6);  total (60/51)
Fisher's exact test for biased strand distribution p-value = 5.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 113 total aligned reads displayed.

ATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATC  >  NC_000913/72056‑72354
                                                                                                                                                     |                                                                                                                                                     
aTTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGtt                                                                                                                                                       >  2:941/1‑150 (MQ=255)
 ttCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATCGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTc                                                                                                                                                      <  2:741650/150‑1 (MQ=255)
 ttCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTc                                                                                                                                                      <  2:1543422/150‑1 (MQ=255)
  tCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCa                                                                                                                                                     <  2:615408/150‑1 (MQ=255)
   cGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCAc                                                                                                                                                    <  2:885769/150‑1 (MQ=255)
   cGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCAc                                                                                                                                                    >  1:1150163/1‑150 (MQ=255)
   cGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCAc                                                                                                                                                    >  1:1223351/1‑150 (MQ=255)
     ccACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCAcgc                                                                                                                                                  >  2:1058377/1‑150 (MQ=255)
     ccACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGAGAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCATTCAcgc                                                                                                                                                  >  1:223300/1‑150 (MQ=255)
      cACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCAcgcc                                                                                                                                                 >  2:7733/1‑150 (MQ=255)
          cGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCaa                                                                                                                                             >  2:168299/1‑149 (MQ=255)
           gTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATc                                                                                                                                            <  2:1504912/150‑1 (MQ=255)
           gTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTGAAcgccttatc                                                                                                                                            <  2:1680395/150‑6 (MQ=255)
            ttGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                           <  2:467130/150‑1 (MQ=38)
            ttGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATcc                                                                                                                                           <  1:1331930/150‑1 (MQ=255)
            tggATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                           <  2:349840/148‑1 (MQ=38)
            gtGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTAAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                           <  2:353093/149‑1 (MQ=38)
             tGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                          <  2:1653973/150‑1 (MQ=38)
             tGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCg                                                                                                                                          <  2:597469/150‑1 (MQ=255)
                 gCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATGCGGCGTGAACGCCTTATCCGAcct                                                                                                                                      <  2:1128671/150‑3 (MQ=38)
                 gCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAg                                                                                                                                      >  2:432848/1‑150 (MQ=255)
                 gCCGGTGTATATCGAAATTTTGCGTAAAGTGGTGGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAg                                                                                                                                      <  1:149001/150‑1 (MQ=255)
                   cGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTg                                                                                                                                    <  1:1373322/150‑1 (MQ=255)
                   cGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTg                                                                                                                                    <  1:1366643/150‑1 (MQ=255)
                   cGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTGTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTg                                                                                                                                    <  1:726986/150‑1 (MQ=255)
                    ggTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGc                                                                                                                                   <  2:1626985/150‑1 (MQ=255)
                        tataTCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATAc                                                                                                                               <  1:1468698/150‑1 (MQ=255)
                        tataTCGATAATTTGCCTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGGGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGCTCACGCCGCATCCGACAGTGCATAc                                                                                                                               <  1:303211/150‑1 (MQ=255)
                         atatCGATTTTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACt                                                                                                                              <  2:1714705/150‑1 (MQ=255)
                            tCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAAc                                                                                                                           >  1:1138086/1‑150 (MQ=255)
                             cGATATTTTGCCTAACGTGGTTGGGGTTAAGGAGATACTCTTGTGCGGGTTGCCTGATGCGACGCTGGCCGCTCCCTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGAGACGGCGTTCACGCCGCATCCGACAGTGCATACTAAcc                                                                                                                          <  1:1059113/150‑1 (MQ=38)
                              gATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAAccc                                                                                                                         >  1:246909/1‑150 (MQ=255)
                              gATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAAccc                                                                                                                         >  1:1360319/1‑150 (MQ=255)
                                aatTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGGTGCCTGATGCGCCGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGt                                                                                                                       <  2:1369277/149‑1 (MQ=255)
                                 atTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTa                                                                                                                      >  1:1131507/1‑150 (MQ=255)
                                    ttGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATc                                                                                                                   <  2:40796/150‑1 (MQ=255)
                                       cGTAAAGTGGTTGGGGTATAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTGCATACTAACCCGTAATcccc                                                                                                                <  1:529471/150‑1 (MQ=255)
                                        gTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCa                                                                                                               <  1:7965/150‑1 (MQ=255)
                                         tAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCaa                                                                                                              >  2:474359/1‑150 (MQ=255)
                                          aaaGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCAGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGTATAATGCGTTCACGCCGCTTCCGTCAGTGCATACTAACCCGTAATCCCCata                                                                                                             >  2:837938/1‑148 (MQ=38)
                                           aaGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCCCGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATa                                                                                                            >  2:1197263/1‑150 (MQ=255)
                                               ggTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAACCGCACTACCCGTTGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTACCCCGTAATCCCCAATAGTGc                                                                                                        >  1:192784/1‑150 (MQ=255)
                                               ggTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATTCTAACCCGTAATCCCCAATAGTGc                                                                                                        >  1:1590480/1‑150 (MQ=255)
                                               ggTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGc                                                                                                        <  1:717692/150‑1 (MQ=255)
                                                 gtGGGGTTTAGGAGATAGGCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCg                                                                                                      <  1:221183/149‑1 (MQ=255)
                                                     ggTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGc                                                                                                  <  2:710269/150‑1 (MQ=255)
                                                         tAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGACGCACTc                                                                                              >  2:837010/1‑150 (MQ=255)
                                                             agaTAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTGCACGCCGCAGCCGACAGTTCATACTAACCCGTAATCCACAATAGTTCCGAAGCACTCGCCt                                                                                          >  1:285289/1‑150 (MQ=38)
                                                                tAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTTCATAACGCACGACCCGTAGGTCGGATAAGGCGGTCACGCCGCATCCGACAGTGCATACTTACCCGGAATCCCCAGTAGTGCCGATGCACTCGCGTTAc                                                                                       >  1:882197/1‑150 (MQ=37)
                                                                tAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTAc                                                                                       >  1:1570860/1‑150 (MQ=255)
                                                                tAGTCTTGGGCGGGTTGCCTGCTGCGACGCTTGCCGCGTCTTATCTGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAGGCACTCGCCTTAc                                                                                       <  1:683143/150‑1 (MQ=255)
                                                                  gTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGACTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCg                                                                                     >  1:1470048/1‑150 (MQ=255)
                                                                    cTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCt                                                                                   <  1:103002/150‑1 (MQ=255)
                                                                      tgtgCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTca                                                                                 <  2:717347/150‑1 (MQ=255)
                                                                          cGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGGAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTccacca                                                                             >  1:182654/1‑148 (MQ=255)
                                                                                gCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCAtt                                                                       <  1:50678/150‑1 (MQ=255)
                                                                                gCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCAtt                                                                       <  2:1570860/150‑1 (MQ=255)
                                                                                gCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCAtt                                                                       >  1:1583328/1‑150 (MQ=255)
                                                                                  cTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAGAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTgg                                                                     >  2:1094543/1‑150 (MQ=255)
                                                                                  cTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACCACTCATTgg                                                                     >  2:885940/1‑150 (MQ=255)
                                                                                  cTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTgg                                                                     <  1:1276621/150‑1 (MQ=255)
                                                                                       ccGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATa                                                                >  1:383537/2‑150 (MQ=255)
                                                                                          aCGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATAccc                                                             >  2:286749/1‑150 (MQ=255)
                                                                                                gCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGATGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAg                                                       <  1:1686313/150‑1 (MQ=255)
                                                                                                 ccGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAgg                                                      <  1:7733/150‑1 (MQ=255)
                                                                                                  cgcgTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGc                                                     <  1:866750/150‑1 (MQ=255)
                                                                                                   gcgTCTTATCAGGCCTACAAAACGCACTACCAGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCATTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCg                                                    >  1:801556/1‑150 (MQ=255)
                                                                                                     gTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGAt                                                  <  1:400493/150‑1 (MQ=255)
                                                                                                         tATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATgcgc                                              >  1:456334/1‑150 (MQ=255)
                                                                                                         tATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACACATTGGTCATACCCTGCCAGGCGATgcgc                                              >  2:671606/1‑150 (MQ=255)
                                                                                                            cAGGCCTACAAAACTCACTACCCGTAGGTCGGATAAGGCGTTCACGCCTCAGCCGACAGTGCAAACGAACCCGGAATCCCCAATCGTGCCGAAGCACTCGCCTGAGCGCTCAACAACTCATTGGAGATACCCAGCCCGGCGAGGCGCCCg                                           >  2:820269/1‑150 (MQ=255)
                                                                                                            cAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCg                                           >  1:166869/1‑150 (MQ=255)
                                                                                                                 cTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCTAAGAACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGc                                      >  2:86016/1‑150 (MQ=255)
                                                                                                                 cTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGc                                      <  2:1360319/150‑1 (MQ=255)
                                                                                                                  tACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCGCGCCGCATCCGACAGTGCAGACTACCCCGTAATCCCCAAAAGTGCCGAAGCACTCGCCTTACCGCTCAACACCTCAGTGGGCAGACCCGGCCCGGCGATGGCCCCGTCGGCg                                     >  1:1068699/1‑150 (MQ=255)
                                                                                                                  tACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTGACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCg                                     >  2:476982/1‑150 (MQ=255)
                                                                                                                     aaaaCGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGact                                  >  2:1607328/1‑150 (MQ=255)
                                                                                                                     aaaaCGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGact                                  <  1:657054/150‑1 (MQ=255)
                                                                                                                      aaaCGCACTACCCGTAGGTCGGATAAGGCGTTCCCGCCGCAGCCGACAGTGCCTACTAACCCGGAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATCCCCTGCCCGGCGATTCGCCCGTCGGCGccta                                 >  1:262593/1‑150 (MQ=255)
                                                                                                                         cGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGActacta                              <  2:451598/150‑1 (MQ=255)
                                                                                                                           cACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTAcc                            <  1:1552279/150‑1 (MQ=255)
                                                                                                                               aCCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAgc                        >  2:1599589/1‑150 (MQ=255)
                                                                                                                                cccGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACcgcgcg                       >  2:178383/1‑150 (MQ=255)
                                                                                                                                  gcactgTCGGATGCGGCGTGAACGCCTTATCCGACCGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCGCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGCCTACTACCGAGcgcg                     >  2:371035/6‑150 (MQ=38)
                                                                                                                                   cactgTCGGATGCGGCGTGAACGCCTTATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGt                    >  1:579213/5‑150 (MQ=255)
                                                                                                                                    tAGGTCGGATAAGGCGTTAACGCCGCATCCGACAGTGCATACTAACCCGTAATCACCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTg                   >  1:919074/1‑150 (MQ=255)
                                                                                                                                     aGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTgg                  >  2:1673932/1‑150 (MQ=255)
                                                                                                                                      ggTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGc                 <  2:246909/150‑1 (MQ=255)
                                                                                                                                      ggTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGc                 <  1:1623136/150‑1 (MQ=255)
                                                                                                                                         cGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGAt              >  2:658263/1‑150 (MQ=255)
                                                                                                                                          ggatgcggCGTGAACGCCTTATCCGACCGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAAAAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGCATACTTCCGAGCGCGTGGCGATc             >  2:306620/7‑150 (MQ=255)
                                                                                                                                          ggATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATc             >  2:1627270/1‑150 (MQ=255)
                                                                                                                                            aTAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCcg           >  1:671501/1‑150 (MQ=255)
                                                                                                                                            aTAAGGCGTTCACGCCGCATACGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTCCCGCTCAACAACTCCTTGGTCATACACTGCCAGGCGATGCGCCCGTCGGCGACAACTACCGAGCGCGTGGCGATCcg           >  1:1156593/1‑150 (MQ=255)
                                                                                                                                              ttgccGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCcgcg         <  2:74396/146‑1 (MQ=255)
                                                                                                                                              aaGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCcgcg         >  2:1376365/1‑150 (MQ=255)
                                                                                                                                                ggCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCACCAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGcgcc       >  1:316397/1‑150 (MQ=255)
                                                                                                                                                ggCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCCTACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTAACGAGCGCGTGGCGATCCGcgcc       >  2:316408/1‑150 (MQ=255)
                                                                                                                                                ggCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGcgcc       >  1:1378151/1‑150 (MQ=255)
                                                                                                                                                     tCACGCCGCATCCGACAGTGCATACTAACCCGAAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCCGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATc  <  2:677153/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATC  >  NC_000913/72056‑72354

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: