Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,070,413 T→G 7.9% intergenic (‑12/+102) yihR ← / ← yihS putative sulphoquinovose mutarotase/sulphoquinovose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,070,4130TG 7.9% 184.8 / 5.9 73intergenic (‑12/+102)yihR/yihSputative sulphoquinovose mutarotase/sulphoquinovose isomerase
Reads supporting (aligned to +/- strand):  ref base T (31/36);  new base G (2/4);  total (33/40)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.16e-01

CCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAG  >  NC_000913/4070265‑4070558
                                                                                                                                                    |                                                                                                                                                  
cccTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCATGTGAAACGGTCTGTCCACTGCAATGCATTTTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTg                                                                                                                                                   >  2:66997/1‑150 (MQ=255)
cccTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTg                                                                                                                                                   <  2:99792/150‑1 (MQ=255)
cccTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTg                                                                                                                                                   <  2:1028284/150‑1 (MQ=255)
    taaaGGTTAATTCAGCCAGACCAGCACCGACGGTCACATTGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCTCTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGc                                                                                                                                               <  1:354395/149‑1 (MQ=255)
         gTTAATTCAGCGAGACCAGCACCGACGGTCACAATGCGCTCGTGATACTCCCCCGCGGCAAGTGAAACGGTCAGTCCACTGCACTGCATATTGGTGATTTGCATACTAGGCACCTTAATTAACGACATACAGCCCCCCGTGAGGCTGATg                                                                                                                                         <  1:1362121/150‑1 (MQ=255)
         gTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATg                                                                                                                                         <  2:840132/150‑1 (MQ=255)
          ttAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTTATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGt                                                                                                                                        <  1:500095/150‑1 (MQ=255)
          ttAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCACGGGAGGCTgtatg                                                                                                                                         <  1:1439709/150‑5 (MQ=255)
            aaTTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCa                                                                                                                                      <  2:556264/150‑1 (MQ=255)
            aaTTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCACGGGAGGCTGTATGTc                                                                                                                                       <  1:783184/150‑1 (MQ=255)
            aaTTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTACAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCa                                                                                                                                      <  2:195497/150‑1 (MQ=255)
              ttCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATCAGCCTCACGGGAGGCTGTATGTcgtt                                                                                                                                    <  2:882584/150‑4 (MQ=255)
                 aGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAACCGGTCTGTCCCCTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGTCTGATGTCAATcca                                                                                                                                 >  1:280023/1‑149 (MQ=255)
                 aGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCACTGCATATTGGTGATTTGCAAACAAGGCACCTTAAATAACGACATCAGCCTCACGGGAGGCTGTATGTcgttaat                                                                                                                                 <  1:1136936/150‑7 (MQ=255)
                 aGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAAccc                                                                                                                                 >  1:1477864/1‑150 (MQ=255)
                    cAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCaaa                                                                                                                              >  2:22865/1‑150 (MQ=255)
                       aCCAGCACCTACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAAc                                                                                                                           <  2:1031995/150‑1 (MQ=255)
                       aCCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAAc                                                                                                                           <  2:1341344/150‑1 (MQ=255)
                        ccAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTGAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACt                                                                                                                          >  1:888996/1‑150 (MQ=255)
                                      cacaATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACACGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAACCCCAACCACTGCACCGGACGGTcc                                                                                                            >  2:1638986/1‑150 (MQ=255)
                                        caATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTcccc                                                                                                          >  2:1439278/1‑150 (MQ=255)
                                         aaTGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGGCTGGCCACTGCAATGCTTATTGGTGATTTGCAGTCAAGGCACCTTAATTATCGCCATACAGCCTCCCGGGAGGCTGATGGGCAACCCGAAAACTGCACCGGCCGGGccccg                                                                                                         >  2:598832/1‑149 (MQ=255)
                                          aTGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTc                                                                                                        <  1:718566/150‑1 (MQ=255)
                                          aTGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAACGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGAACCGGACGGTCCCCta                                                                                                        >  2:823567/1‑149 (MQ=255)
                                                  gTGATAATCCCCCGCGGCAAGTGAAACGATCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCttt                                                                                                >  2:343811/1‑150 (MQ=255)
                                                    gATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCTTATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgg                                                                                              <  1:1683654/150‑1 (MQ=255)
                                                     aTAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg                                                                                             <  2:1685131/150‑1 (MQ=255)
                                                     aTAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg                                                                                             <  2:55744/150‑1 (MQ=255)
                                                       aaTCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAACCCCAAAACCTGCACCGGACGGTCCCCTCGCCCCGTTGGGga                                                                                           >  2:43837/1‑150 (MQ=255)
                                                       aaTCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGga                                                                                           <  1:1626076/150‑1 (MQ=255)
                                                       aaTCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCCCCCCGTTGGGga                                                                                           >  2:1584011/1‑150 (MQ=255)
                                                          ccccccCGGCAAGGGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACATTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGagag                                                                                        <  2:1072175/150‑1 (MQ=255)
                                                          cccccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGagag                                                                                        <  2:530807/150‑1 (MQ=255)
                                                           ccccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCTCCCCAAAGGGGCGAgg                                                                                       <  1:1480400/150‑1 (MQ=255)
                                                           ccccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAgg                                                                                       <  1:922221/150‑1 (MQ=255)
                                                           ccccGCGGCAAGTGAAACGGTCTGACCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAATGACATACAGCCTCCCGTTAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAgg                                                                                       <  1:734446/150‑1 (MQ=255)
                                                             ccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGCCGGTCCCCTCGCCCCTTTGGGGAGAGGGt                                                                                     <  2:893290/150‑1 (MQ=255)
                                                             ccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGt                                                                                     <  2:1019905/150‑1 (MQ=255)
                                                               gcgGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTAGCATACAAGGCACCTTAATTAACGACATACAGCCACCCGTGAGGCGGATGTCAAACCCAAAAACTGCCCCGGACGGTCCCCTCTCCCCAAAGGGGCGAggggac                                                                                   <  1:1703546/150‑4 (MQ=255)
                                                                cgGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAg                                                                                  <  1:399861/150‑1 (MQ=255)
                                                                 ggCAAGGGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                                                                                 <  1:1567288/150‑1 (MQ=255)
                                                                  gCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAggg                                                                                <  1:468500/150‑1 (MQ=255)
                                                                    aaGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGgag                                                                              <  2:1449598/150‑1 (MQ=255)
                                                                      gTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGagag                                                                            >  1:1030065/1‑150 (MQ=255)
                                                                      gTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCGTTGGGGAGAGGGTTAGGGagag                                                                            >  1:1733590/1‑150 (MQ=255)
                                                                        gAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAggg                                                                          >  2:189004/1‑150 (MQ=255)
                                                                         aaaCGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGCGAgggg                                                                         >  2:752570/1‑150 (MQ=255)
                                                                         aaaCGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAgggg                                                                         <  2:39635/150‑1 (MQ=255)
                                                                          aaCGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGa                                                                        <  1:22865/150‑1 (MQ=255)
                                                                              gTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAc                                                                    >  2:1618545/1‑150 (MQ=255)
                                                                              gTCTGTCCACAGCAATGCATATTGGTGATTTGCATACAAGGCACCATAATTAACGACATACAGCCACCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCTGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAc                                                                    <  1:164372/150‑1 (MQ=255)
                                                                                  gTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCTTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCTCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGc                                                                <  2:195734/150‑1 (MQ=255)
                                                                                    ccACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCAc                                                              >  2:763757/1‑150 (MQ=255)
                                                                                    ccACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCCCAc                                                              >  2:1120072/1‑150 (MQ=255)
                                                                                        tGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCATAAACTGCACCCGACGGACCCCTCGCCCCTTTGGGGAGAGGGATAGGGAGAGGGGAAAACCCGCACCCTa                                                          >  1:396362/1‑150 (MQ=255)
                                                                                         gCAATGCATATTGGTGATTTGAATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGGTGTCACACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAt                                                         >  1:1165955/1‑150 (MQ=255)
                                                                                          cAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGAAAACCCGCACCCTAttt                                                       <  1:1472258/150‑1 (MQ=255)
                                                                                           aaTGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTGTGGTGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAttt                                                       >  1:255223/1‑150 (MQ=255)
                                                                                            aTGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtttt                                                      >  2:189871/1‑150 (MQ=255)
                                                                                                    tGGTGACTTGCATACAAGCCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAACCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAg                                              <  1:275701/150‑1 (MQ=255)
                                                                                                              cATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGACAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCGATTTTCAGATCAGTTATTTCGCa                                    >  2:852872/1‑150 (MQ=255)
                                                                                                                 aCAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCCCCCTATTTTCAGATCAGTTATTTCGCATTa                                 >  2:139558/1‑150 (MQ=255)
                                                                                                                       cACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATcc                           >  1:1128623/1‑150 (MQ=255)
                                                                                                                           ctaattaaCGACATACAGCCTCCCGTGAGGCTGATGGCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCcagca                       <  1:291026/149‑1 (MQ=255)
                                                                                                                            tatttaaCGAAATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTTCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCcagcag                      <  2:280023/147‑1 (MQ=255)
                                                                                                                            taattacCGACATCCAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACCGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCcagcag                      <  1:23842/150‑1 (MQ=255)
                                                                                                                                taaCGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTAGTTCGCATTAATATCCAGCAGAccc                  >  2:1446720/1‑150 (MQ=255)
                                                                                                                                 aaCGACATACAGCCTCCCGTGGGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCg                 >  2:997404/1‑150 (MQ=255)
                                                                                                                                      cATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCa            >  2:1476287/1‑150 (MQ=255)
                                                                                                                                      cATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGAAGGTCCCCTCGCCCCTTTGGGGAGAGGGATAGGGAGAGGGGAAAACCCGCACCCAATATTCAGCTCAGTAATTTCGCAGTAATATCCAGCAGACCCGCTGCa            >  1:1152458/1‑150 (MQ=39)
                                                                                                                                      acatcAGCCTCACGGGAGGCTGTATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGc             >  2:432306/5‑150 (MQ=39)
                                                                                                                                       aTACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGCTCAGTGAGTTCGCATTAAGAGCCAGCAGACCCGCTGCaa           >  2:89442/1‑150 (MQ=39)
                                                                                                                                                cccGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTAATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAg  >  1:872091/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAG  >  NC_000913/4070265‑4070558

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: