Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 72,199 A→C 6.9% intergenic (+84/+30) 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,1990AC 6.9% 313.7 / 9.9 112intergenic (+84/+30)yabI/thiQDedA family inner membrane protein/thiamine/thiamine pyrophosphate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (57/47);  new base C (2/6);  total (59/53)
Fisher's exact test for biased strand distribution p-value = 1.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 113 total aligned reads displayed.

TGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCG  >  NC_000913/72051‑72347
                                                                                                                                                    |                                                                                                                                                    
tgtTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAg                                                                                                                                                     <  2:940226/150‑1 (MQ=255)
 gtTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAgg                                                                                                                                                    >  2:1553423/1‑150 (MQ=255)
 gtTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAgg                                                                                                                                                    >  2:1546734/1‑150 (MQ=255)
    aaTTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGt                                                                                                                                                 <  1:1208034/150‑1 (MQ=255)
    aaTTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGCCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGt                                                                                                                                                 >  2:976626/1‑150 (MQ=255)
     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)
          ccACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCAcgc                                                                                                                                           >  2:1058377/1‑150 (MQ=255)
          ccACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCAcgc                                                                                                                                           <  2:217196/150‑1 (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)
                                 tCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAAc                                                                                                                    >  1:1138086/1‑150 (MQ=255)
                                  cGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAAcc                                                                                                                   >  2:1201469/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)
                                                        ggggTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGaa                                                                                             <  2:819658/150‑1 (MQ=255)
                                                          ggTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGc                                                                                           <  2:710269/150‑1 (MQ=255)
                                                                  agaTAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCt                                                                                   >  2:1623136/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:1239835/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)
                                                                                                                           aaaCGCACTACCCGTAGGTCGGATAAGGCGTTCCCGCCGCAGCCGACAGTGCCTACTAACCCGGAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATCCCCTGCCCGGCGATTCGCCCGTCGGCGccta                          >  1:262593/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)
                                                                                                                                     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          <  1:1623136/150‑1 (MQ=255)
                                                                                                                                           ggTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGc          <  2:246909/150‑1 (MQ=255)
                                                                                                                                            gTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACCACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCg         >  1:810684/1‑150 (MQ=255)
                                                                                                                                              cGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGAt       >  2:658263/1‑150 (MQ=255)
                                                                                                                                               ggATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATc      >  2:1627270/1‑150 (MQ=255)
                                                                                                                                                 aTAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCcg    >  1:671501/1‑150 (MQ=255)
                                                                                                                                                 aTAAGGCGTTCACGCCGCATACGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTCCCGCTCAACAACTCCTTGGTCATACACTGCCAGGCGATGCGCCCGTCGGCGACAACTACCGAGCGCGTGGCGATCcg    >  1:1156593/1‑150 (MQ=255)
                                                                                                                                                   aaGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCcgcg  >  2:1376365/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCG  >  NC_000913/72051‑72347

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