Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 72,206 C→A 8.7% intergenic (+91/+23) 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,2060CA 8.7% 310.4 / 13.6 110intergenic (+91/+23)yabI/thiQDedA family inner membrane protein/thiamine/thiamine pyrophosphate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (56/44);  new base A (4/6);  total (60/50)
Fisher's exact test for biased strand distribution p-value = 5.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 114 total aligned reads displayed.

TTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCT  >  NC_000913/72057‑72355
                                                                                                                                                     |                                                                                                                                                     
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)
    ccACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCAcgc                                                                                                                                                   <  2:217196/150‑1 (MQ=255)
    ccACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGAGAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCATTCAcgc                                                                                                                                                   >  1:223300/1‑150 (MQ=255)
         cGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCaa                                                                                                                                              >  2:168299/1‑149 (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)
                            cGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAAcc                                                                                                                           >  2:400493/1‑150 (MQ=255)
                            cGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAAcc                                                                                                                           >  2:1201469/1‑150 (MQ=255)
                            cGATATTTTGCCTAACGTGGTTGGGGTTAAGGAGATACTCTTGTGCGGGTTGCCTGATGCGACGCTGGCCGCTCCCTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGAGACGGCGTTCACGCCGCATCCGACAGTGCATACTAAcc                                                                                                                           <  1:1059113/150‑1 (MQ=38)
                             gATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAAccc                                                                                                                          >  1:1360319/1‑150 (MQ=255)
                             gATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAAccc                                                                                                                          >  1:246909/1‑150 (MQ=255)
                               aatTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGGTGCCTGATGCGCCGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGt                                                                                                                        <  2:1369277/149‑1 (MQ=255)
                                atTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTa                                                                                                                       >  1:1131507/1‑150 (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)
                                                gtGGGGTTTAGGAGATAGGCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCg                                                                                                       <  1:221183/149‑1 (MQ=255)
                                                  ggggTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGaa                                                                                                     <  2:819658/150‑1 (MQ=255)
                                                    ggTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGc                                                                                                   <  2:710269/150‑1 (MQ=255)
                                                        tAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGACGCACTc                                                                                               >  2:837010/1‑150 (MQ=255)
                                                            agaTAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCt                                                                                           >  2:1623136/1‑150 (MQ=255)
                                                            agaTAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTGCACGCCGCAGCCGACAGTTCATACTAACCCGTAATCCACAATAGTTCCGAAGCACTCGCCt                                                                                           >  1:285289/1‑150 (MQ=38)
                                                            agaTAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGGAGGTCGGATCAGGCGTTCACGCCGACTCCGACAGTGCATACTAACCCTGAATCCCCCATAGTGCCGCAGCCCTCACCt                                                                                           >  1:289235/1‑150 (MQ=37)
                                                               tAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTTCATAACGCACGACCCGTAGGTCGGATAAGGCGGTCACGCCGCATCCGACAGTGCATACTTACCCGGAATCCCCAGTAGTGCCGATGCACTCGCGTTAc                                                                                        >  1:882197/1‑150 (MQ=37)
                                                               tAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTAc                                                                                        >  1:1570860/1‑150 (MQ=255)
                                                               tAGTCTTGGGCGGGTTGCCTGCTGCGACGCTTGCCGCGTCTTATCTGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAGGCACTCGCCTTAc                                                                                        <  1:683143/150‑1 (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)
                                                                                 cTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTgg                                                                      >  1:1487254/1‑150 (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:1715600/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                                   <  1:657054/150‑1 (MQ=255)
                                                                                                                    aaaaCGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGact                                   >  2:1607328/1‑150 (MQ=255)
                                                                                                                        cGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGActacta                               <  2:451598/150‑1 (MQ=255)
                                                                                                                         gCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACtactac                              >  1:651161/1‑150 (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)
                                                                                                                                    aGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTgg                   >  2:1673932/1‑150 (MQ=255)
                                                                                                                                     ggTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGc                  <  1:1623136/150‑1 (MQ=255)
                                                                                                                                      gTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACCACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCg                 >  1:810684/1‑150 (MQ=255)
                                                                                                                                        cGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGAt               >  2:658263/1‑150 (MQ=255)
                                                                                                                                         ggatgcggCGTGAACGCCTTATCCGACCGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAAAAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGCATACTTCCGAGCGCGTGGCGATc              >  2:306620/7‑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)
                                                                                                                                                     cACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCt  <  2:1583328/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCT  >  NC_000913/72057‑72355

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: