Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,070,463 T→A 5.9% intergenic (‑62/+52) 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,4630TA 5.9% 199.6 / 2.6 78intergenic (‑62/+52)yihR/yihSputative sulphoquinovose mutarotase/sulphoquinovose isomerase
Reads supporting (aligned to +/- strand):  ref base T (41/31);  new base A (3/2);  total (45/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

CGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGT  >  NC_000913/4070314‑4070609
                                                                                                                                                     |                                                                                                                                                  
cgTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCTATTTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGACAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCtt                                                                                                                                                    >  1:1602683/1‑150 (MQ=255)
cgTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCtt                                                                                                                                                    <  2:1441889/150‑1 (MQ=255)
 gTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCttt                                                                                                                                                   <  2:1131240/150‑1 (MQ=255)
  tGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                                                                  <  2:1524755/150‑1 (MQ=255)
    aTAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg                                                                                                                                                <  2:935292/150‑1 (MQ=255)
      aaTCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggggt                                                                                                                                              >  1:1261069/1‑149 (MQ=255)
      aaTCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGga                                                                                                                                              >  2:463868/1‑150 (MQ=255)
      aaTCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGga                                                                                                                                              <  2:1075837/150‑1 (MQ=255)
         cccccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGagag                                                                                                                                           <  1:1028346/150‑1 (MQ=255)
            ccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGt                                                                                                                                        <  2:451903/150‑1 (MQ=255)
            ccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGt                                                                                                                                        <  2:444738/150‑1 (MQ=255)
             cgcgGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGtt                                                                                                                                       <  1:1591283/150‑1 (MQ=255)
                ggCAAGTGAAACGGTCTGTCCACTTCAATGCATATTGGTTATTTGCATACAAGTCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTACCCTCTCCCCAAAGGGGCGAggggaccg                                                                                                                                    <  2:739839/150‑6 (MQ=255)
                ggCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTACCCTCTCCCCAAAGGGGCGAggggaccg                                                                                                                                    <  2:733086/150‑6 (MQ=255)
                 gCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAggg                                                                                                                                   <  1:27792/150‑1 (MQ=255)
                   aaGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGgag                                                                                                                                 <  2:1511317/150‑1 (MQ=255)
                   aaGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGgag                                                                                                                                 <  2:1515279/150‑1 (MQ=255)
                     gTGACACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGagag                                                                                                                               <  1:1896438/150‑1 (MQ=255)
                     gTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGagag                                                                                                                               >  1:309184/1‑150 (MQ=255)
                       gAAACGGTCTGTCCACTGCAATTCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAggg                                                                                                                             <  2:1062394/150‑1 (MQ=255)
                        aaaCGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGTCACCTAAATTAACCACATACAGCCACCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGGCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAgggg                                                                                                                            <  1:1892676/150‑1 (MQ=39)
                         aaCGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGATAGGGTTAGGGAGAGGGGa                                                                                                                           >  2:917623/1‑150 (MQ=255)
                            ggTCTGTCCTCTGCAATGCATATTGTTGTTTTGCATACAAGGCACCTTAATTAACGACCTACAGCCGCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGTGAGGGTTAGGGAGAGGGGaaaa                                                                                                                        <  2:914321/150‑1 (MQ=38)
                             gTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACCGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAc                                                                                                                       <  1:467439/150‑1 (MQ=255)
                                   ccACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACAGCACCGGACGCTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCAc                                                                                                                 >  2:1059408/1‑150 (MQ=255)
                                    cACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCAcc                                                                                                                >  2:1041058/1‑150 (MQ=255)
                                     aCTTCAATTAATATTGGTGATTTCCATACAAGGCACCTTATATAACTACATACAGCCTCCCGGTAGGCTGATGTAAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTTGTGAGAGGGGAAAACCCGCAccc                                                                                                               <  2:20131/150‑1 (MQ=255)
                                          aaTGCAAATTGGTGATTTGCATACAAGGCACCTTAATTAACGACTTCCAGCCTCCCGTGAGGCTGATGTCAAACCAAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAttt                                                                                                          <  2:1344550/150‑1 (MQ=255)
                                           aTGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtttt                                                                                                         >  2:1771482/1‑150 (MQ=255)
                                              cATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAg                                                                                                      >  2:1089918/1‑150 (MQ=255)
                                                tatTGGTGATTTGAATAAAAGGCACCTTAATTAACCACATACACCCTCCCGTCAGGCTTATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGAt                                                                                                    <  1:985335/150‑1 (MQ=38)
                                                  ttGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGAACCCTATTTTCAGATCa                                                                                                  >  2:553638/1‑150 (MQ=255)
                                                   tGGTGATTTGCATACAAGGCACCTTAATTAACGCCATACAGCCTCCCGTGAGGCTGATGTCACACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAg                                                                                                 >  2:323348/1‑150 (MQ=255)
                                                       gATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTAGGCGAGACGGTTAGGGAGCGGGGAAAACCCGCCCCAGATTTTCAGAGAAGTTAt                                                                                             >  2:1667925/1‑150 (MQ=255)
                                                        ttttGCTTACAAGGCACAGTAATTAACGTCATACATGCTCCCGTGTGGATTATGTCAAACCCAAAAACTTCACCGGCCGGTCCCCTCGCCCCTTTTGGGTGTTGGTTTTGTAGAGGGGAAAACCCGCAACCTATTTTCAGATCAGTTAtt                                                                                            <  2:232988/149‑1 (MQ=255)
                                                        aTTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTAtt                                                                                            >  1:614279/1‑150 (MQ=255)
                                                         tttGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTAttt                                                                                           >  2:1799974/1‑150 (MQ=255)
                                                           tGCCTACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTGTGGGGAGCGGGTTAGGGAGCGGGGAAAACCCGCCCCCGATTTGCAGATCCGTTAGTTCg                                                                                         >  2:1443822/1‑150 (MQ=255)
                                                             cATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGGGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCa                                                                                       <  2:1659085/150‑1 (MQ=255)
                                                              aTACACGGCACCTTAATTAACGACAAACAGCCTCCCGTGAGCCTGATGTCAAACCCAAAAACTGCACCTGACGGTCCCCTCGCCCCTTTGGGGAGATGGTTAGGTAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCAt                                                                                      <  1:289747/150‑1 (MQ=255)
                                                                 cAAGGCACCTTAATAAACGACAAACCGCCGCCCGTGTGGCTTATGTAAAACCAAAAAACTCCACCGGACGGTCCCCTCGCCCCTTTGGGGAGTGGGTTTGGTAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTaa                                                                                   <  2:561851/150‑1 (MQ=255)
                                                                   aGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCTTTAatt                                                                                 >  2:1516223/1‑149 (MQ=255)
                                                                   aGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTAAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAata                                                                                 <  2:5653/150‑1 (MQ=255)
                                                                       aCCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACTGTCCCCTCGCCCATTTGGGGAGAGGGTTAGGGAGAGGGGAAACCCCGCACCCGATGTTCAGATCAGTTATTTCGCATTAATATCca                                                                             >  1:164716/1‑150 (MQ=255)
                                                                            aattaaCGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAACACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGa                                                                        >  1:213242/1‑150 (MQ=255)
                                                                             attaaCGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTAGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGAc                                                                       >  1:1484336/1‑150 (MQ=255)
                                                                                aaCGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCg                                                                    >  2:615919/1‑150 (MQ=255)
                                                                                  cGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCt                                                                  >  2:842021/1‑150 (MQ=255)
                                                                                     caaaCAGCCACCCGTGAGGCTGATGTAAAACCAAAAAACGGCACCGGACGGTCCCCTCGCCCCTTTGGGGTGAGGGTTCGGGAGCGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCa                                                               <  2:708713/147‑1 (MQ=255)
                                                                                      aTACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAACCCCGCCCCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCaa                                                              >  1:1561861/1‑150 (MQ=255)
                                                                                       atcAGCCTCACGGGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGCAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGCCCCGCTGCAAc                                                             >  1:247216/3‑150 (MQ=255)
                                                                                        aCAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGTGGGTTAGGGAGAGGGGAAAACCCGCACCCAATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAccc                                                            >  2:1168358/1‑150 (MQ=255)
                                                                                        aCAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCCCCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAAcc                                                            >  1:78980/1‑150 (MQ=255)
                                                                                         cAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCg                                                           >  2:1485793/1‑150 (MQ=255)
                                                                                           gCCTCACGGGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCCCCCTATTTTCAGATCAGTTTTTTCGCATTAATATCCAGCAGACCCGCTGCACCCGCt                                                         >  1:1090065/1‑150 (MQ=39)
                                                                                             cTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTAATTCGCATTAAAATCCAGCAGCCCCGCTGCAACCGCTgg                                                       >  1:1706830/1‑150 (MQ=255)
                                                                                                ccGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCACCCGCTGGAGc                                                    >  2:455769/1‑150 (MQ=255)
                                                                                                 cGTGAGGCTTATGTAAAACCCAAAAACAGCACCGGACGGTCCCCTCGCCCGTTGGGGGAGTGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGcc                                                   <  1:556187/150‑1 (MQ=255)
                                                                                                        cTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCgg                                            >  1:522593/1‑150 (MQ=255)
                                                                                                         tGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCGCCTTTGGGGAGAGGGTTAGGGATAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCggg                                           <  2:544627/150‑1 (MQ=255)
                                                                                                           aTGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGgcg                                         <  1:1053847/150‑1 (MQ=255)
                                                                                                             gTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAACCCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCt                                       >  2:484135/1‑150 (MQ=255)
                                                                                                                 aaCCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACg                                   >  2:1352110/1‑150 (MQ=255)
                                                                                                                   cccAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAACCCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCCACCGCTGGAGCCATGCCGGGCGCTAACggg                                 >  2:902829/1‑150 (MQ=255)
                                                                                                                   cccAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACggg                                 >  2:16104/1‑150 (MQ=255)
                                                                                                                    ccAAAAACTGCACCGGCCGGTCCCCTCGCCCCTATGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATATACAGATCAGTTATTTCGCAGTAATAACCAGCAGACCCGCTGCAACCGCTGGAGCCAGTCCGGGCTCTAACgggt                                >  2:1152935/1‑149 (MQ=255)
                                                                                                                          aCTCCACCGTACTTTCCCCTCGCCCATTTGGGTAGTGTTTTTGTGATTGGTGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATAcgcg                          <  1:917623/150‑1 (MQ=255)
                                                                                                                           cTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCAGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCgg                         >  2:1788561/1‑150 (MQ=255)
                                                                                                                            ggCACCGGCCGGTCCCCTCCCCCCTTTGGGGCGTGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCggg                        <  2:1706830/149‑1 (MQ=255)
                                                                                                                                     cctaacccTCTCCCCAAAGGGGCGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCCCCAgg               >  2:103162/6‑150 (MQ=255)
                                                                                                                                     cGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGCGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAgg               >  2:1621693/1‑150 (MQ=255)
                                                                                                                                     cGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAgg               <  1:1240051/150‑1 (MQ=255)
                                                                                                                                        tCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGTGATCACCAGGCaa            >  2:1909699/1‑150 (MQ=255)
                                                                                                                                         ccccTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAAt           >  1:45715/1‑150 (MQ=255)
                                                                                                                                         acccTCTCCCCAAAGGGGCGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAAt           >  1:1522147/2‑150 (MQ=255)
                                                                                                                                           ccTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATgc         <  2:614279/150‑1 (MQ=255)
                                                                                                                                              cGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATgcagc      >  2:1910318/1‑150 (MQ=255)
                                                                                                                                                ccccTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGcagcag    >  2:240604/1‑150 (MQ=255)
                                                                                                                                                 cccTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAATGGGATACGCGGGATCACCAGGCAATGCAGCAgg   >  2:456963/1‑150 (MQ=255)
                                                                                                                                                  ccTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCTGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGt  >  1:1719630/1‑150 (MQ=255)
                                                                                                                                                  ccTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGt  >  1:1425181/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGT  >  NC_000913/4070314‑4070609

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: