Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,322,705 C→T 6.0% intergenic (‑200/+28) secG ← / ← glmM preprotein translocase membrane subunit/phosphoglucosamine mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,322,7050CT 6.0% 367.2 / 7.8 117intergenic (‑200/+28)secG/glmMpreprotein translocase membrane subunit/phosphoglucosamine mutase
Reads supporting (aligned to +/- strand):  ref base C (49/61);  new base T (5/2);  total (54/63)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.26e-01

Only 100 of 126 total aligned reads displayed.

AAGCGGGAGTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACGCAGCAACACGCGTC  >  NC_000913/3322556‑3322854
                                                                                                                                                     |                                                                                                                                                     
aaGCGGGAGTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATcgc                                                                                                                                                       <  1:869507/150‑1 (MQ=255)
 agatggAGTTTTTAATCTTTCCCACTGAATCGGGTTTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATcgcc                                                                                                                                                      <  1:1415946/146‑1 (MQ=255)
 aGCGGGAGTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATcgcc                                                                                                                                                      >  2:278778/1‑150 (MQ=255)
          ttttAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAACCAGCAAAACCGATGTGTTAGCTGAGAAAAACGGCGACAATCGCCGCCTTTTTa                                                                                                                                             >  1:361256/1‑150 (MQ=255)
          ttttAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTa                                                                                                                                             >  1:1234251/1‑150 (MQ=255)
           tttAATCTTTACCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTat                                                                                                                                            >  2:770506/1‑149 (MQ=255)
             tAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGcc                                                                                                                                          <  1:1164445/150‑1 (MQ=255)
             tAATCTATCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCttttttctc                                                                                                                                          <  2:689195/150‑5 (MQ=255)
               aTCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAg                                                                                                                                        <  1:794666/150‑1 (MQ=255)
                     cccACTGAAGCGGGTTTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCTATTTTATCGCTTTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCt                                                                                                                                  <  1:1635265/150‑1 (MQ=255)
                     cccACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCt                                                                                                                                  >  2:899156/1‑150 (MQ=255)
                        aCTGAAGCGGGTGTGAATTCTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAAc                                                                                                                               <  2:1681163/150‑1 (MQ=255)
                                gggTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTaaa                                                                                                                       >  1:769608/1‑150 (MQ=255)
                                gggTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTaaa                                                                                                                       >  2:483019/1‑150 (MQ=255)
                                    gtgAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCCTTTTAGCCAGTTATCTAACGCTTTAAACGGc                                                                                                                   >  1:1138829/1‑150 (MQ=255)
                                     tgAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCtt                                                                                                                 <  1:1725642/150‑1 (MQ=255)
                                      gAATACTACCCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCTAGTGCGGAAAAAAACTGCAAAACCGATGTGTTCGCTGAGAAAAACGGCGACAATCGCGGCGTTTTTAGCCAGTTATCTTACGCTTGAAGCGGCtt                                                                                                                 >  2:1856700/1‑150 (MQ=255)
                                      gAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCtt                                                                                                                 >  2:689022/1‑150 (MQ=255)
                                      gAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCtt                                                                                                                 >  2:669369/1‑150 (MQ=255)
                                      gAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCtt                                                                                                                 >  2:26605/1‑150 (MQ=255)
                                        aTACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCtttt                                                                                                               <  1:588412/150‑1 (MQ=255)
                                         tACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTa                                                                                                              <  1:1864824/150‑1 (MQ=255)
                                          aCTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTAc                                                                                                             <  1:1614383/150‑1 (MQ=255)
                                            tAACCAAAGCGAATGACCTTCGCAAGGGCATTTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTg                                                                                                           <  1:1190384/150‑1 (MQ=255)
                                            tAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTg                                                                                                           <  2:122441/150‑1 (MQ=255)
                                              aCCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCa                                                                                                         >  2:1503152/1‑150 (MQ=255)
                                               ccAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAACCGGCTTTTACTGCAt                                                                                                        <  1:290235/150‑1 (MQ=255)
                                               ccAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCAt                                                                                                        >  2:899716/1‑150 (MQ=255)
                                               ccAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCAt                                                                                                        >  2:1536266/1‑150 (MQ=255)
                                               ccAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCAt                                                                                                        >  2:2013631/1‑150 (MQ=255)
                                                cAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATc                                                                                                       <  1:877849/150‑1 (MQ=255)
                                                  aaGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCgg                                                                                                     <  1:326423/150‑1 (MQ=255)
                                                   aGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGc                                                                                                    <  1:1005253/150‑1 (MQ=255)
                                                     cGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTATTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGcta                                                                                                  >  2:1248555/1‑148 (MQ=255)
                                                       actGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGTAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTAATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATg                                                                                                <  2:1385910/148‑1 (MQ=255)
                                                       aaTGACCTTCGCAAGGGCAATTTTAACGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTTCATCGGCGATg                                                                                                <  1:330559/150‑1 (MQ=255)
                                                        aTGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGc                                                                                               <  1:1916766/150‑1 (MQ=255)
                                                         tGACCTTCGCAAGGTCAATTTTATCGCATTGTATCAACTGCGGAAAAACTCCGCAAAACCGATGTGTTCGATGTGAAAAAAGGCGACAAACGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCg                                                                                              <  1:363022/150‑1 (MQ=255)
                                                            ccTTCGGAATGGAAATTTTATCGCATTTTTTCAACTTCGGTAAAAAACAGGAAAACCGTTGGGTTTGCTTTGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGgtg                                                                                           <  2:721484/150‑1 (MQ=255)
                                                              ttCGCAAGGGCAATTTTATCGCTTTGTATCACCTGCGGAAAAAACCAGCAAAACCGATGTGTTCGCTGAGAAAAAACGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGtgtg                                                                                         <  2:293816/150‑1 (MQ=255)
                                                              ttCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAACCAGCAAAACCGATGGGTTCGCTGAGAAAAACGGCGCCAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGtgtg                                                                                         <  2:1689072/150‑1 (MQ=255)
                                                              ttCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGtgtg                                                                                         <  1:435554/150‑1 (MQ=255)
                                                              ttCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGtgtg                                                                                         <  2:1059050/150‑1 (MQ=255)
                                                                 gCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCaa                                                                                      >  2:539905/1‑150 (MQ=255)
                                                                 gCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCaa                                                                                      <  1:1306277/150‑1 (MQ=255)
                                                                   aaGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAAt                                                                                    >  2:1104797/1‑150 (MQ=255)
                                                                    aGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAAtt                                                                                   <  1:992043/150‑1 (MQ=255)
                                                                        cAATATTATCGCCTTGTATCACCTGCGGAAAAAAACCGCAAAACCGATGTTTTAGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGt                                                                               <  2:304156/150‑1 (MQ=255)
                                                                         aaTTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTc                                                                              <  1:76983/150‑1 (MQ=255)
                                                                              tATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTGTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGGAAATTCAGTCACCTg                                                                         >  2:800007/1‑150 (MQ=255)
                                                                               aTCGCAATGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAATAAAGGCGCCAATCGCCGCCTTATTAGCCAGTTATCGAACGCTTTAAACTGCATATACTGCATCGGCGATCCGGTGGGCAAATTCAGTCACCTgc                                                                        >  1:313156/1‑150 (MQ=255)
                                                                                tCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgcg                                                                       <  1:1464762/150‑1 (MQ=255)
                                                                                 cGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgcgc                                                                      >  1:943452/1‑150 (MQ=255)
                                                                                     ttGTATCAACTGCGGATAAAAACAGCAATACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCg                                                                  <  2:707232/150‑1 (MQ=255)
                                                                                     ttGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCg                                                                  <  2:1283216/150‑1 (MQ=255)
                                                                                       gTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTACAAATTCAGTCACCTGCGCTTCGTc                                                                >  1:1996040/1‑150 (MQ=255)
                                                                                        tATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACTCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCCCCTGCGCTTCGTCt                                                               >  1:1229929/1‑150 (MQ=255)
                                                                                           cAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCg                                                            <  1:1998780/150‑1 (MQ=255)
                                                                                           cAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCg                                                            <  2:1234251/150‑1 (MQ=255)
                                                                                           cAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCg                                                            >  1:424965/1‑150 (MQ=255)
                                                                                            aaCTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGc                                                           <  2:1364127/150‑1 (MQ=255)
                                                                                             aCTGCGGAAAAAAACCGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGcc                                                          <  2:984786/150‑1 (MQ=255)
                                                                                             aCTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGcc                                                          <  1:1874243/150‑1 (MQ=255)
                                                                                                 cGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTcc                                                      <  1:1554210/150‑1 (MQ=255)
                                                                                                    aaaaaaaCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTccacc                                                   >  2:806552/1‑150 (MQ=255)
                                                                                                      aaaaaCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCAt                                                 <  1:1080550/150‑1 (MQ=255)
                                                                                                       aaaaCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATc                                                >  1:485962/1‑150 (MQ=255)
                                                                                                        aaaCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCa                                               >  2:1450903/1‑150 (MQ=255)
                                                                                                           cAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCTCCTTCCACCATCAcgc                                            >  1:1268358/1‑150 (MQ=255)
                                                                                                               aaaaCCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAAt                                        >  1:632674/1‑150 (MQ=255)
                                                                                                               aaaaCCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAAt                                        <  1:483019/150‑1 (MQ=255)
                                                                                                               aaaaCCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGAGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAAt                                        >  1:1879044/1‑150 (MQ=255)
                                                                                                                        gtgtTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTc                               <  2:1174787/150‑1 (MQ=255)
                                                                                                                            tCGCTGAGAATAAAGGCGACAATCGCCGCCTTTTTAGGCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTg                           <  2:1006001/150‑1 (MQ=255)
                                                                                                                            tCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAATTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCATGCGAATTAACGGTTCGGTg                           <  1:384713/150‑1 (MQ=255)
                                                                                                                             cGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGc                          >  1:1808368/1‑150 (MQ=255)
                                                                                                                             cGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGc                          >  2:199261/1‑150 (MQ=255)
                                                                                                                               ctggctaaaaaGGCGGCGATTGTCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCg                        >  2:1147409/7‑150 (MQ=255)
                                                                                                                               cTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCg                        >  1:1429479/1‑150 (MQ=255)
                                                                                                                                 ggctaaaaaGGCGGCGATTGTCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGa                      >  2:1844740/5‑150 (MQ=255)
                                                                                                                                 gagaAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGa                      <  1:1020123/150‑1 (MQ=255)
                                                                                                                                  gctaaaaaGGCGGCGATTGTCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGAt                     >  1:1354899/4‑150 (MQ=255)
                                                                                                                                  agaAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGAt                     >  2:1984524/1‑150 (MQ=255)
                                                                                                                                  agaAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGAt                     <  2:663460/150‑1 (MQ=255)
                                                                                                                                   ctaaaaaGGCGGCGATTGTCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGAtt                    >  1:1757721/3‑150 (MQ=255)
                                                                                                                                      aaaaGGCGGCGATTGTCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGt                                 >  1:741415/1‑134 (MQ=255)
                                                                                                                                      aaaaGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTAc                 >  1:1052664/1‑150 (MQ=255)
                                                                                                                                        aaGGCGACACTCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACgc               <  2:1582393/150‑1 (MQ=255)
                                                                                                                                         aGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACgcc              >  1:496506/1‑149 (MQ=255)
                                                                                                                                         aGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACgca              <  1:2013631/150‑1 (MQ=255)
                                                                                                                                         aGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACgca              <  1:899716/150‑1 (MQ=255)
                                                                                                                                            cGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACgcagca           >  1:473238/1‑150 (MQ=255)
                                                                                                                                            cGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACgcagca           >  2:1123640/1‑150 (MQ=255)
                                                                                                                                            cGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGACTTCGCCATCCACCATCACGCGATTTAACGGTTCGGTGCCGCATTTACgcagca           >  1:117753/1‑150 (MQ=255)
                                                                                                                                              aCAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACGCAGCAac         <  1:1845627/150‑1 (MQ=255)
                                                                                                                                                    gccgcgTTTTTAGCCAGTTATCTAACGCTTTAATCGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCAGCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACGCAGCAACAcgcg  >  2:1322288/1‑149 (MQ=255)
                                                                                                                                                    gccgccTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGGAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCCTCACGCGAAGTACCGGTTCGGTGCCGGATTGCCGCCGCAACCcgcg  >  2:227362/1‑149 (MQ=255)
                                                                                                                                                    gccgccTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACGCAGCAACACGCGt   >  1:645837/1‑150 (MQ=255)
                                                                                                                                                     ccgccTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACGCAGCAACACGCGTc  >  1:1604702/1‑150 (MQ=255)
                                                                                                                                                     ccgccTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACGCAGCAACACGCGTc  <  2:816749/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAGCGGGAGTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACGCAGCAACACGCGTC  >  NC_000913/3322556‑3322854

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: