Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,150,903 A→C 8.7% L732R (CTG→CGG)  ppc ← phosphoenolpyruvate carboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,150,9030AC 8.7% 163.5 / 3.4 68L732R (CTG→CGG) ppcphosphoenolpyruvate carboxylase
Reads supporting (aligned to +/- strand):  ref base A (31/28);  new base C (2/5);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTCCGCCAGCCACAGGTCTGCTTTGGCGAAGACCATCTCCAGCATGCCGAGACGCGTCGAGAAGAATGGCCAATCGCGGCACATAGCCTCCAGCTCGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGAAGATCCACGGAATGGCGCGTAGTGACTCGACGCCGCCGGTTGGGCGACGTTTCGCCGGACGTGAACCCAACGGCAGTTTGCCCAGTTCTTGTTCCGGCGTAGCGGAGCG  >  NC_000913/4150756‑4151051
                                                                                                                                                   |                                                                                                                                                    
tATTCCGCTATCCACAGGTCTGCTTTGCCGTAGCCCACCGCCAGCATGCCGCGACGCGTCGAGAAGAATGGCCAATCGCGGCACATAGCCTCCAGCTCGCTCTGTTTGCCGTCTTCTATCACTTTTTGCAGCGCCGTACCTGCACccag                                                                                                                                                     <  2:622205/149‑1 (MQ=255)
tATTCCGCCAGCCACCGGTCTGCTTTGGCGAAGACCATCTCCAGCATGCCGAGACGCGTCGAGAAGAATGGCCAATCGCGGCACATAGCCTCCAGCTCGCTCTGTTTGCCGTCTTCGACCCCTTTTTGCAGCGCCGTACCTGCACccag                                                                                                                                                     <  2:452874/149‑1 (MQ=255)
  ttCCGCCAGCCACAGGTCTGCTTTGGCGAAGACCATCTCCAGCATGCCGAGACGAGTGGAGAAGAATGACCAAGCGCGTCACATGGCCTCCAGCTCGCTCTGGTCGGCGTCGTCGACCACTTTTGGGAGCGCCGTACCTGCACccagcc                                                                                                                                                   >  2:327958/1‑149 (MQ=255)
   tCCGCCAGCCTCAGGTCTCCTGTGGCGAATGCCATCTCCAGCATGCCGAGACGCGTCGAGAAGAATGGCCAATCGCGGCACATAGCCTCCAGCGCGCTCTGTTTGCCGTCTTAAAACAGATATTTCAGTGCCGTACCTGCACccagcca                                                                                                                                                  <  1:590562/149‑1 (MQ=255)
       ccagccaCAGGTCTGCTTTGGCGAAGACCATCTCCAGCATGCCGAGACGCGTCGAGAAGAATGGCCAATCGCGGCACATAGCCTCCAGCTCGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCg                                                                                                                                              >  1:622597/1‑149 (MQ=255)
        cagccaCAGGTCTGCTTTGGCGAAGACCATCTCCAGCATGCCGAGACGCGTCGAGAAGAATGGCCAATCGCGGCACATAGCCTCCAGCTCGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCgg                                                                                                                                             >  1:531972/1‑149 (MQ=255)
         agccaCAGGTCTGCTTTGGCGAAGACCATCTCCAGCATGCCGAGACGCGTCGAGAAGAATGGCCAATCGCGGCACATAGCCTCCAGCTCGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCggg                                                                                                                                            <  1:102692/149‑1 (MQ=255)
         agccaCAGGTCGGCTGTGGCGAAGACCATCTCCAGCATGCCGAGATGCGTCGAGAAGAATGGCCAATCGCGGCACATAGCCTCCAGCTAGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCAAGCGCGGCggg                                                                                                                                            <  2:336368/149‑1 (MQ=255)
           ccaCAGGTCTGCTTTGGCTAAGCCCATTTCCAGCATTCCGAGACGCGTCGAGAAGAATGGCCAATCGCGGCACATAGCCTCCAGCTCGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGa                                                                                                                                          <  1:184104/149‑1 (MQ=255)
            cacaGGTCTGCTTTGGCGAAGACCATCTCCAGCATGCCGAGACGCGGCGAGAAGAATGGCCAATCGCGGCACATAGCCTCCAGCTCGCGCTGTTTGCCGTCTTCGACCACGTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCggggcg                                                                                                                                         >  2:211882/1‑147 (MQ=255)
             acaGGTCTGCTTTGGCGAAGACCATCTCCAGCATGCCGAGACGCGTCGAGAAGAATGGCCAATCGCGGCACATAGCCTCCAGCTCGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCTCCCCGCCTGGCTGGGTGc                                                                                                                                        <  2:587038/149‑1 (MQ=255)
               cgttCTGCTTTGGCGAAGACCATTTCCCGCATGCCGAGACGCGTCGAGAAGAATGGCCAATCGCGGTACATAGCCTCCAGCTCGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCTCCCCGCCTGGCTGGGTGCag                                                                                                                                      <  2:225695/146‑2 (MQ=255)
               cggTCTGCTTTGGCGAGGACCATCTCCAGCATGCCGAGACGCGTCGAGTAGAATGGCCAATCGCGGCACATAGCCTCCCGCTCGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCTCCCCGCCTGGCTGGGTGCag                                                                                                                                      <  2:690965/148‑2 (MQ=255)
                ggTCTGCTTTGGCGAAGACCATCTCCAGCATGGCGAGACGCGGCCGGAAGAATGCCCAATCGCGGCACATAGCCTCCAGCTCGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATc                                                                                                                                     <  2:513542/149‑1 (MQ=255)
                 gTCTGCTTTGGCGAAGACCATCTCCAGCATGCCGAGACGCGTCGAGAAGAATGGCCAATCGCGGCACATAGCCTCCAGCTCGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCa                                                                                                                                    <  1:490087/149‑1 (MQ=255)
                   cTGCGTGGGCGAAGACCATCTCCTGCCTGCCGAGACGCGTCGAGAAGAATGGCCAATCGCGGCACATAGCCTCCAGCTCGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGa                                                                                                                                  <  2:121790/149‑1 (MQ=255)
                     gCTTTGGCGAAGACCATCTCCAGCAGGCCGGGACGCGTCGAGAAGAATGGCCAATCGCGGCACATAGCCTCCAGCGCGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACg                                                                                                                                <  2:13318/149‑1 (MQ=255)
                      cTTTGGCGAAGACCATCTCCAGCATGCCGAGACGCGTCGAGAAGAATGGCCAATCGCGGCACATAGCCTCCAGCTCGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACgg                                                                                                                               <  1:434217/149‑1 (MQ=255)
                      cTTTGGCGAAAACCATCCCCAGCAAACCCAGAGTGGGCGGGAAGCATGGCCAACCGCCGCACATTGCCTCCAGCTCGCTCTGTTTGCCGTCGCCGACCACCTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACgg                                                                                                                               <  2:327583/149‑1 (MQ=255)
                           gCGAAGACCATCTCCAGCATGCCGAGACGCGTCGAGAAGAATGGCCAATCGCGGCACATAGCCTCCAGCTCGCGCTGTTTGCCGTCTTCGACCAATATTAGCTTCTCCGGACCTGCCCCCCGCCAGTCGGGGAGCATCAGACGGTTTTg                                                                                                                          >  1:379622/1‑149 (MQ=255)
                               aGACCATCTCCAGCATGCCGAGACGCGTCGAGAAGAATGGCCAATCGCGGCACATAGCCTCCAGCTCGCTCTGTTTGCCGTCTTCGCCCACTTTTTGCCGCGCCGTACCTGCACCCCGCCAGGCGTGGAGCATCAGACGGTTTTGCGTc                                                                                                                      <  1:430525/149‑1 (MQ=255)
                                gACCATCTCCAGCATGCCGAGACGCGTCGAGAAGAATGGCCAATCGCGGCACATAGCCTCCAGCTCGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGGGCCGTACCTGCACCCCGCCAGGCGGGGAGCATCAGACGGTTTTGCGTcc                                                                                                                     >  1:494094/1‑149 (MQ=255)
                                gACCATCTCCAGCATGCCGAGACGCGTCGAGAAGAATGGCCAATCGCGGCACATAGCCTCCAGCTCGCTCTGTTTACCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTcc                                                                                                                     <  2:476776/149‑1 (MQ=255)
                                 aCCTTCTCCAGCATGCCGAGACGCGTCGAGAAGAATGGCCAATCGCGGCACATAGCCTCCAGCTCGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCa                                                                                                                    <  2:262021/149‑1 (MQ=255)
                                 aCCATCTCCAGCATGCCGAGACGCGTCGAGAAGAATGGCCAATCGCGGCCCATAGCCTCCAGCTCGCGCTGTTTGTCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCa                                                                                                                    >  2:135076/1‑149 (MQ=255)
                                  ccATCTCCAGCATGCCGAGACGCGTCGAGAAGAATGGCCAATCGCGGCACATAGCCTCCAGCTCGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCAg                                                                                                                   <  1:696048/149‑1 (MQ=255)
                                  ccATCTCCAGCATGCCGAGACGCGTCGAGAAGAATGGCCAATCGCGGCACATAGCCTCCAGCTCGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCAg                                                                                                                   <  1:357320/149‑1 (MQ=255)
                                       tcCACCCTGGCGAGACGCGTCGAGAAGAATGGCCAATCGCGCCACATAGCCTCCAGCTCGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGaa                                                                                                              <  1:211359/149‑1 (MQ=255)
                                        ccAGCATGCCGAGACGCGTCGAGAAGAATGGCCAATCGCGGCACATAGCCTCCAGCTCGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGAAg                                                                                                             >  2:431199/1‑149 (MQ=255)
                                           gCATGCCGAGACGCGTCGAGAAGAATGGCCAATCGCGGCACATAGCCTCCAGCTCGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACGGTTGTGCGTCCAGGCGaggttc                                                                                                          >  1:126453/1‑144 (MQ=255)
                                           gCAGGCCGAGACGCGTCGAGAAGAATGGCCAATCGCGGCACATAGCCTCCAGCTCGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGAAGATc                                                                                                          <  1:762773/149‑1 (MQ=255)
                                              tGCCGAGACGCGTCGAGAAGAAGGGCCAATCCCGGCACATAGCCTCCAGCTCGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGCGCATCAGACGGTTTCGCGTCCAGGCGAAGATCCAc                                                                                                       <  1:644048/149‑1 (MQ=255)
                                                  gagaCGCGTCGAGAAGAATGGCCAATCGCGGCACATAGCCTCCAGCTCGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCGGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGAAGATCCACGGaa                                                                                                   >  2:620650/1‑149 (MQ=255)
                                                    gaCGCGTCGAGAAGAATGGCCAACCGCGGCACATAGCCTCCAGCTCGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGAAGATCCACGGAATg                                                                                                 <  1:559494/149‑1 (MQ=255)
                                                     aCGCGTCGAGAAGAATGGCCAATCGCGGCACATAGCCTCCAGCTCGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCCGGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGAAGATCCACGGAATgg                                                                                                >  2:533376/1‑149 (MQ=255)
                                                          tCGAGAAGAATGGCCAATCGCGGCACATAGCCTCCAGCTCGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGAAGATCCACGGAATGGCGCGt                                                                                           >  1:797454/1‑149 (MQ=255)
                                                             agaagaATGGCCAATCGCGGCACATAGCCTCCAGCTCGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGAAGATCCACGGAATGGCGCGTAgg                                                                                        >  1:155431/1‑148 (MQ=255)
                                                             agaagaATGGCCAATCGCGGCACATAGCCTCCAGCTCGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGAAGATCCACGGAAAGGCGCGTAGt                                                                                        >  1:155435/1‑149 (MQ=255)
                                                                 gaaTGGCCAATCGCGGCACATAGCCTCCAGCTCGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGGGCATCAGACGGTTTTGCGACCAGGCGAAGATCCACGGACTGGCGCGTAGTGACt                                                                                    >  2:488854/1‑149 (MQ=255)
                                                                  aaTGGCCAATCGCGGCACATAGCCTACAGCTCGCTATGTTTGCCGCCTTCGACCACTTTTTGCAGCGCCGTACATGCACCGAGCCAGGCGGGGAGCAGCAGACGGTTTTGCGTCCAGGCGAAGATCCACGGAATGGCCCGTAGTTAcgc                                                                                   >  1:391140/1‑147 (MQ=255)
                                                                       ccAATCGCGGCACATAGCCTCCAGCTCGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGTAGATCCACGGGATGGCGCGGAGTGACTCGAcgc                                                                              >  1:342183/1‑149 (MQ=255)
                                                                        cAATCGCGGCACATAGCCTCCAGCTCGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGAAGATCCACGGAATGGCGCGTAGTGACTCGAcgcc                                                                             >  2:273746/1‑149 (MQ=255)
                                                                         aaTCGCGGCACATAGCCTCCAGCTCGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGAAGATCCACGGAATGGCGCGTAGTGACGCGAcgccg                                                                            >  1:202124/1‑149 (MQ=255)
                                                                          aTCGCGGCACATAGCCTCCAGCTCGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGAAGAGCCACGGAATGGCGCGTAGTGACTCGAcgccgc                                                                           <  2:32454/149‑1 (MQ=255)
                                                                            cgcgGCACGTGGCGTGCAGCTGGGTCGGTTTGCCGTCTTCGACCGCTGTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGAAGATCCACGGAATGGCGGGTAGTGACTCGACGccgccg                                                                         <  2:184961/149‑1 (MQ=255)
                                                                                  acaTAGCCTCCAGCTCGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGAAGATCCACGGAATGGCGCGTAGTGACTCGACGCCGCCGGTTggg                                                                   >  2:255549/1‑149 (MQ=255)
                                                                                  acaTAGCCTCCAGCTCGCTCTGTTTACAGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGAAGATCCACGGAATGGCGCGTAGTGACTCGACGCCGCCGGTTggg                                                                   >  1:283117/1‑149 (MQ=255)
                                                                                    aTAGCCTCCAGCTCGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGAAGATCCACGGAATGGCGCGTAGTGACTCGACGCCGCCGGTTGGGCg                                                                 >  1:172991/1‑149 (MQ=255)
                                                                                            cAGCTCGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGAAGATCCACGGAATGGCGCGTAGTGACTCGACGCCGCCGGTTGGGCGACGTTTCg                                                         >  1:409950/1‑149 (MQ=255)
                                                                                                tcgctcTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGAAGATCCACGGAATGGCGCGTAGTGACTCGACGCCGCCGGTTGGGCGACGTTTCGCCgg                                                     >  1:82702/1‑149 (MQ=255)
                                                                                                 cgctcTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGAAGATCCACGGAATGGCGCGTAGTGACTCGACGCCGCCGGTTGGGCGACGTTTCGCCGGa                                                    <  2:338222/149‑1 (MQ=255)
                                                                                                      tGTTTGCGGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCTGGTGGGCGGGGGGCGTCGGCCGGTTTTGCGTCCAGGCGAAGATCCAAGGGATGGCGCGTAGTGACTCGAGGCCGCCGGGTGGGGGCCGTTTCGCCGGGCGTGa                                               >  2:402119/1‑149 (MQ=255)
                                                                                                         gtGCCGTCCTCGACCACTTTTTGCAGCGCCGTACCTGCACCCCGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGAAGATCCACGGAATGGCGCGTAGTGACTCGACGCCGCCGGTTGGGCGACGTTTCGCCGGACGTGAAcc                                            <  2:391442/148‑1 (MQ=255)
                                                                                                            ccGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGAAGATCCACGGAATGGCGCGTAGTGACTCGACGCCGCCGGTTGGGCGACGTTTCGCCGGACGTGAACCCaa                                         >  1:318373/1‑149 (MQ=255)
                                                                                                             cGTCTTCGACCCATTTTTTCAGCGGCGTACCTGCCCCCAGCCAGGCGGGGAGGATCAGACGGTTTTGCGTCCAGGCGAAGATCCACGGAATGGCGCGTAGTGACTCGACGCCGCCGGTTGGGCGACGTTTCGCCGGACGTGAACCCAAc                                        <  2:126453/149‑1 (MQ=255)
                                                                                                                   cGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGGCCAGGCGAAGATCCACGGAATGGCACGTAGTGACTCGACGCCGCCGGGTGGGCGACGTTTCGACGGACGTGCACCCACCGGCAGt                                  >  1:599824/1‑149 (MQ=255)
                                                                                                                      ccACTTTTTGCCGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGAAGATCCACGGAATGGCGCGTAGTGACTCGACGCCGCCGGTTGGGCGACGTTTCGCCGGACGTGAACCCAACGGCAGTTTg                               <  2:52706/149‑1 (MQ=255)
                                                                                                                       cACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGAAGATCCACGGAATGGCGCGTAGTGACTCGACGCCGCCGGTTGGGCGACGTTTCGCCGGACGTGAACCCAACGGCAGTTTGc                              >  1:130422/1‑149 (MQ=255)
                                                                                                                       cACTTTTCGCAGCGCCGTACCTGCTCCCAGCCAGGCGGGGAGCATCAGCCGGTTTTGCGTCCAGGCGAAGATCCACGGAATGGCGCGTAGTGACTCGACGCCGCCGGTTGGGCGACGTTTCGCCGGACGTGAACCCAACGGCAGTTTGc                              <  2:612021/149‑1 (MQ=255)
                                                                                                                        aCTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGAAGATCCACGGAATGGCGCGTAGTGACTCGACGCCGCCGGTTGGGCGACGTTTCGCCGGACGTGAACCCAACGGCAGTTTGcc                             <  1:620650/149‑1 (MQ=255)
                                                                                                                         gttgttGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGAAGATCCACGGAATGGCGCGTAGTGACTCGACGCCGCCGGTTGGGCGACGTTTCGCCGGACGTGAACCCAACGGCCGTTTGccc                            >  2:534873/5‑149 (MQ=255)
                                                                                                                         cTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGAAGATCCACGGAATGGCGCGTAGTGACTCGACGCCGCCGGTTGGGCGACGTTTCGCCGGACGTGAACCCAACGGCAGTTTGccc                            <  2:417897/149‑1 (MQ=255)
                                                                                                                          tttttGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGAAGATCCACGGAATGGCGCGTAGTGACTCGACGCCGCCGGTTGGGCGACGTTTAGCCGGACGTGAACCCAACGGCAGTTTGCCCa                           >  1:537865/1‑149 (MQ=255)
                                                                                                                                cAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGAAGATCCACGGAATGGCGCGTAGTGACTCGACGCCGCCGGTTGGGCGACGTTTCGCCGGACGTGAACCCAACGGCAGTTTGCCCAGTTCtt                     <  1:740004/149‑1 (MQ=255)
                                                                                                                                 aGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGAAGATCCACGGAATGGCGCGTAGTGACTCGACGCCGCCGGTTGGGCGACGTTTCGCCGGACGTGAACCCAACGGCAGTTTGCCCAGTTCTTg                    >  1:451505/1‑149 (MQ=255)
                                                                                                                                  gcgcCGTACCTGCCCCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGAAGATCCACGGAATGGCGCGTAGTGACTCGACGCCGCCGGTTGGGCGACGTTTCGCCGGACGTGAACCCAACGGCAGTTTGCCCAGTTCTTGt                   <  2:391140/149‑1 (MQ=255)
                                                                                                                                    gcCGTACCTGCACCCAGCCAAGCGGGGGGCATCAGACGGCTTTGCGGCCAGGCTAGGATTCACGCAATGGCGCGAAGTGACTAGACGCCGCCGGCTGAGCCGAGTCTCCCCGGACGTGAACCCAACGGCAGTTTGCCCCGttgtctgtc                 >  2:281501/1‑142 (MQ=255)
                                                                                                                                           cTGCACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGAAGATCCACGGAATGGCGCGTAGTGACTCGACGCCGCCGGTTGGGCGACGTTTCGCCGGACGTGAACCCAACGGCAGTTTGCCCAGTTCTTGTTCCGGCGtg          >  2:89790/1‑148 (MQ=255)
                                                                                                                                           cTGCACCCAGCCAGGCGGGGAGCAGCAGACGGACCTGCGTCGAGGAGAACAACCGCGGGGGGTCGCGTCATGACCCGCCGCCGCCGGGCGGGCGCCGTTTCGCCGGGCGCGAACCCAACGGCAGGCCGCCCCGTTCTTCTTCCGGCGTa          >  2:331145/1‑149 (MQ=255)
                                                                                                                                            tGCACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGAAGATCCACGGAATGGCGCGTAGTGACTCGACGCCGCCGGTTGGGCGACGTTTCGCCGGACGTGAACCCAACGGCAGTTTGGCCAGTTCTTGTTCCGGCGTag         >  1:61679/1‑149 (MQ=255)
                                                                                                                                              cACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGAAGCTCCCCGGACTGGCGCGTAGCGACTCGACGCCGCCGGTTGGGCGACGTTTCGCCGGGCGTGAACCCAACGGCAGTTTGCCCCGTGCTTGCTCCGGGGTggcg       >  1:94093/1‑149 (MQ=255)
                                                                                                                                                cccAGCCAGGCGGGGAGCATCAGACGGTTTCGCGTCCAGGCGAAGATCCACGGAATGGCGCGTAGTGACTCGACGCCGCCGGTTGGGCGACGTTTCGCCGGACGTGAACCCAACGGCAGTTTGCCCAGTTCTTGTTCCGGCGTagcgga     <  2:479111/149‑1 (MQ=255)
                                                                                                                                                 ccagccagGCGGGGAGCATCAGACGGTTTTGCGTCCGGGCGAAGATCCACGGAATGGCGCGTAGTGACTCGACGCCGCCGGTTGGGCGACGTTTCGCCGGACGTGAACCCAACGGCAGTTTGCCCAGTTCTTGTTCCGGCGTagcggag    >  1:414352/1‑149 (MQ=255)
                                                                                                                                                 ccagccagGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGCAGATCCACGGAATGGCGCGTAGTGACTCGACGCCGCCGGTTGGGCGACGTTTCGCCGGACGTGAACCCAACGGCAGTTTGCCCAGTTCTTGTTCCGGCGTagcggag    >  2:747534/1‑149 (MQ=255)
                                                                                                                                                 ccagccagGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGAAGATCCACGGAATGGCGCGTAGTGACTCGACGCCGCCGGTTGGGCGACGTTTCGCCGGACGTGAACCCAACGGCAGTTTGCCCAGTTCTTGTTCCGGCGTcgcggag    >  1:341655/1‑149 (MQ=255)
                                                                                                                                                  cagccagGCGGGGCGCATCAGACGGTTTTGCGTCCAGGCGAAGATCCACGGAATGGCGCGTAGTGACTCGACGCCGCCGGTTGGGCGACGTTTCGCCGGACGTGAACCCAACGGCAGTTTGCCCAGTTCCTGTTCCGGCGTagcggagc   <  1:211882/149‑1 (MQ=255)
                                                                                                                                                   agccagGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGAAGATCCACGGAATGGCGCGTAGTGACTCGACGCCGCCGGTTGGGCGACGTTTCGCCGGACGTGAACCCAACGGCAGTTTGCCCAGTTCTTGTTCCGGCGTagcggagcg  <  1:273746/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TATTCCGCCAGCCACAGGTCTGCTTTGGCGAAGACCATCTCCAGCATGCCGAGACGCGTCGAGAAGAATGGCCAATCGCGGCACATAGCCTCCAGCTCGCTCTGTTTGCCGTCTTCGACCACTTTTTGCAGCGCCGTACCTGCACCCAGCCAGGCGGGGAGCATCAGACGGTTTTGCGTCCAGGCGAAGATCCACGGAATGGCGCGTAGTGACTCGACGCCGCCGGTTGGGCGACGTTTCGCCGGACGTGAACCCAACGGCAGTTTGCCCAGTTCTTGTTCCGGCGTAGCGGAGCG  >  NC_000913/4150756‑4151051

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: