Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,652,133 A→C 8.2% V52G (GTG→GGG) ‡ yfhM ← alpha2‑macroglobulin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,652,1330AC 8.2% 268.3 / 13.9 96V52G (GTG→GGG) ‡yfhMalpha2‑macroglobulin
Reads supporting (aligned to +/- strand):  ref base A (45/43);  new base C (6/2);  total (51/45)
Fisher's exact test for biased strand distribution p-value = 2.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

AATAACGCGTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGT  >  NC_000913/2651985‑2652276
                                                                                                                                                    |                                                                                                                                               
aaTAACGCGTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCAc                                                                                                                                                <  2:229719/150‑1 (MQ=255)
  tAACGCGTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACgg                                                                                                                                              >  2:281802/1‑150 (MQ=255)
     cgcgTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGgaga                                                                                                                                           >  2:857763/1‑150 (MQ=255)
        gtgtGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGTGCCGGagagac                                                                                                                                        <  1:1447495/150‑3 (MQ=255)
         tGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGTGCCGGagagaca                                                                                                                                       <  1:331645/150‑4 (MQ=255)
            gaAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCGATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGc                                                                                                                                    >  2:804592/1‑150 (MQ=255)
            gaAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGc                                                                                                                                    >  2:606935/1‑150 (MQ=255)
              aaTCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTg                                                                                                                                  <  2:908145/150‑1 (MQ=255)
                tCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGaa                                                                                                                                >  1:1397462/1‑150 (MQ=255)
                 ccTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAAc                                                                                                                               <  2:1127787/150‑1 (MQ=255)
                       ccGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGt                                                                                                                         >  1:1154656/1‑150 (MQ=255)
                          ggTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGtttt                                                                                                                      <  1:1280671/150‑1 (MQ=255)
                             cGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTg                                                                                                                   >  1:1261067/1‑150 (MQ=255)
                               agagGGATGGAGAACGTCAGCACCAGCGTGGCGGCACCATCAAGTGGGACTTCAGAGACATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGCAGAGTTTTGCTGAACTCGCGCTTTCTGaa                                                                                                                 <  2:1284506/150‑1 (MQ=255)
                                gagGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAag                                                                                                                >  2:1389983/1‑150 (MQ=255)
                                 agGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAaga                                                                                                               <  1:606935/150‑1 (MQ=255)
                                  gggATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAagag                                                                                                              >  1:1301751/1‑150 (MQ=255)
                                   ggATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGTGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAgg                                                                                                             >  2:268828/1‑150 (MQ=255)
                                   ggATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTAAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAgg                                                                                                             >  1:1317353/1‑150 (MQ=255)
                                    gATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGc                                                                                                            >  1:955224/1‑150 (MQ=255)
                                    gATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGc                                                                                                            >  1:1157101/1‑150 (MQ=255)
                                      tGGATAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCg                                                                                                          >  1:627927/1‑150 (MQ=255)
                                        gagaACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCt                                                                                                        <  1:1328111/150‑1 (MQ=255)
                                        gagaACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCt                                                                                                        >  2:1016587/1‑150 (MQ=255)
                                         agaACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCtt                                                                                                       >  1:379115/1‑150 (MQ=255)
                                         agaACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCtt                                                                                                       >  1:810024/1‑150 (MQ=255)
                                          taaCGTCATCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCttt                                                                                                      <  2:379845/149‑1 (MQ=255)
                                             cGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGt                                                                                                   <  2:1350693/150‑1 (MQ=255)
                                              gTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTa                                                                                                  >  2:106059/1‑150 (MQ=255)
                                              gTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTa                                                                                                  >  1:1276449/1‑150 (MQ=255)
                                               tCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTaa                                                                                                 >  1:662670/1‑150 (MQ=255)
                                                  gCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACtt                                                                                              >  2:112204/1‑150 (MQ=255)
                                                   cACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTCCCGGCACTCTGGTGGTCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGTACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTc                                                                                             <  1:508812/150‑1 (MQ=255)
                                                    aCCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCa                                                                                            <  2:934022/150‑1 (MQ=255)
                                                     ccAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAg                                                                                           >  1:493151/1‑150 (MQ=255)
                                                     ccAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAg                                                                                           >  1:670472/1‑150 (MQ=255)
                                                        gCGTGGCTGCACCATCAAGTTTGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAg                                                                                        <  2:1216273/150‑1 (MQ=255)
                                                        gCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAg                                                                                        <  2:1545486/150‑1 (MQ=255)
                                                        gCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAg                                                                                        >  1:831886/1‑150 (MQ=255)
                                                        gCGTGGCTGCAACATCAAGTTGGACTTCAGAGAGATCCAGCAATTTAATAAAATTAGCTGCACTCTGTTTGGAGATTTTATGTCTCTACGTCACTTAGAGTTTTGCTGTAATATTGTTTTCTTAATAGTCCGATTTAGTAACTTaagaag                                                                                        <  2:1205255/150‑1 (MQ=255)
                                                         cGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAgg                                                                                       >  2:385166/1‑150 (MQ=255)
                                                           tGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGgcg                                                                                     <  1:1015335/150‑1 (MQ=255)
                                                            ggCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGgcgc                                                                                    <  1:1016587/150‑1 (MQ=255)
                                                                gCACCATCAAGTTGGACTTCAGAGAGTTCCAGGAATTTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCTCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCt                                                                                <  2:1110576/150‑1 (MQ=255)
                                                                 cACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCCTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCtt                                                                               >  2:1623780/1‑150 (MQ=255)
                                                                     aTCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCtttttt                                                                           <  2:21660/150‑1 (MQ=255)
                                                                      tCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTc                                                                          <  1:857763/150‑1 (MQ=255)
                                                                           ttGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGc                                                                     <  1:1503474/150‑1 (MQ=255)
                                                                            tGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCt                                                                    >  2:1466100/1‑150 (MQ=255)
                                                                                 ttCAGAGAGATCCAGCGATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTgg                                                               <  1:804592/150‑1 (MQ=255)
                                                                                  tCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGgc                                                              <  2:379115/150‑1 (MQ=255)
                                                                                   cAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGGCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGgcg                                                             >  2:327524/1‑150 (MQ=255)
                                                                                   cAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGgcg                                                             >  2:114779/1‑150 (MQ=255)
                                                                                     gagagaTCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGcgcg                                                           >  2:991642/1‑150 (MQ=255)
                                                                                     gagagaTCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGcgcg                                                           >  1:774335/1‑150 (MQ=255)
                                                                                      agagaTCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGt                                                          <  2:1175424/150‑1 (MQ=255)
                                                                                         gaTCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTAt                                                       >  1:964179/1‑150 (MQ=255)
                                                                                           tCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCg                                                     <  1:340188/150‑1 (MQ=255)
                                                                                            ccAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCTTCTTTTTTCACCGCTGTTGGCGCGTTATCgt                                                    <  1:733645/150‑1 (MQ=255)
                                                                                             cAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCgtt                                                   <  2:604868/150‑1 (MQ=255)
                                                                                             cAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCgtt                                                   >  1:270995/1‑150 (MQ=255)
                                                                                              aGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCgttg                                                  <  1:1083165/150‑1 (MQ=255)
                                                                                              aGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCgttg                                                  >  2:1376272/1‑150 (MQ=255)
                                                                                               gCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCgttgt                                                 >  1:1261847/1‑150 (MQ=255)
                                                                                               gCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCgttgt                                                 <  1:1506541/150‑1 (MQ=255)
                                                                                                     tCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGc                                           <  1:1244223/150‑1 (MQ=255)
                                                                                                     tCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGc                                           <  1:281802/150‑1 (MQ=255)
                                                                                                      cAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCa                                          >  1:290261/1‑150 (MQ=255)
                                                                                                       aGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCAc                                         >  2:86437/1‑150 (MQ=255)
                                                                                                       aGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCAc                                         >  2:884196/1‑150 (MQ=255)
                                                                                                            cTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCACATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTgc                                    >  1:914009/1‑150 (MQ=255)
                                                                                                                   gCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTgccagcgcc                             >  2:945796/1‑150 (MQ=255)
                                                                                                                     aCTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTgccagcgccag                           <  1:1623780/150‑1 (MQ=255)
                                                                                                                       tctGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCa                         <  2:831886/150‑1 (MQ=255)
                                                                                                                         tGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCAtt                       <  2:1207366/150‑1 (MQ=255)
                                                                                                                          gTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTa                      <  1:1031468/150‑1 (MQ=255)
                                                                                                                           ttGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTGGCACCCTGGCAGCGCCAGCATTAg                     <  2:893979/150‑1 (MQ=255)
                                                                                                                           ttGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAg                     <  1:1528834/150‑1 (MQ=255)
                                                                                                                             gggCCAGTTTTTGTCTCTCCGGCACGGATAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCGTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGTCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCa                   <  1:1089094/150‑1 (MQ=255)
                                                                                                                              ggCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTa                      >  1:717259/1‑146 (MQ=255)
                                                                                                                              ggCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTa                      <  2:717259/146‑1 (MQ=255)
                                                                                                                                ccAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGc                <  1:106059/150‑1 (MQ=255)
                                                                                                                                ccAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCACTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGc                >  1:128834/1‑150 (MQ=255)
                                                                                                                                     gcaaaactctcCGTGCCGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAggcg           >  2:41808/7‑150 (MQ=255)
                                                                                                                                     gcaaaactctcCGTGCCGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAggcg           >  2:1083051/7‑150 (MQ=255)
                                                                                                                                     gcaaaactctcCGTGCCGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAggcg           >  1:524184/7‑150 (MQ=255)
                                                                                                                                      caaaactctcCGTGCCGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAggcgg          >  2:872688/6‑150 (MQ=255)
                                                                                                                                      caaaactctcCGTGCCGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAggcgg          >  1:323206/6‑150 (MQ=255)
                                                                                                                                       ttGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAggcggc         <  1:1376272/150‑1 (MQ=255)
                                                                                                                                         gTCTCTCCGGCACGGAGAGTTTTGGTGAACTGGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATGGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTa       <  2:1317353/150‑1 (MQ=255)
                                                                                                                                         gTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTa       >  1:1036278/1‑150 (MQ=255)
                                                                                                                                         aactctcCGTGCCGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTa       >  2:1092255/3‑150 (MQ=255)
                                                                                                                                           ctctcCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTAcg     >  2:205815/1‑150 (MQ=255)
                                                                                                                                             ctcCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTAcgcg   >  1:1098959/1‑150 (MQ=255)
                                                                                                                                             ctcCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTAcgcg   <  2:1276449/150‑1 (MQ=255)
                                                                                                                                              tcCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGt  <  2:579042/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
AATAACGCGTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGT  >  NC_000913/2651985‑2652276

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