Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,263,315 T→C 5.1% intergenic (+354/+260) yhaC → / ← rnpB pentapetide repeats‑related protein/RNase P, M1 RNA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,263,3150TC 5.1% 530.1 / 16.4 199intergenic (+354/+260)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base T (58/126);  new base C (7/3);  total (70/129)
Fisher's exact test for biased strand distribution p-value = 1.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 201 total aligned reads displayed.

CGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCA  >  CP009273/3263169‑3263462
                                                                                                                                                  |                                                                                                                                                   
cgcATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTa                                                                                                                                                    <  2:511358/148‑1 (MQ=255)
cgcATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTa                                                                                                                                                    <  1:342023/148‑1 (MQ=255)
  cATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTaaa                                                                                                                                                  >  1:163553/1‑148 (MQ=255)
  cATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGAGTCGTGGACAGATATGCATCTATGAGGGCCTCCAAAAGTaaa                                                                                                                                                  >  2:1161969/1‑148 (MQ=25)
   aTCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTaaaa                                                                                                                                                 <  2:2945804/148‑1 (MQ=255)
   aTCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTaaaa                                                                                                                                                 <  2:1202549/148‑1 (MQ=255)
     ccGGCATGAATAAAGCGCACTCAACAATCTGTAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGTTAAGCCGCTTTCTTTTGGGTTTAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAAc                                                                                                                                               <  2:3089430/148‑1 (MQ=255)
     ccGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATATTTCTGCCCCTCCAAAAGTAAAAAc                                                                                                                                               >  2:337219/1‑148 (MQ=255)
     ccGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCGTTCTTTTGGGTATAGTGTCGTTGTCAGTCATTCATATTTCTGCCCCACCAAAACTaaaaaa                                                                                                                                               >  1:1316296/1‑147 (MQ=21)
      gggCATGAATAAAGCGCACTCAACAATCTGAAACCCGACTGAACGGTTTATTACGTACATCAGGTAATACTGACAGATAAGCCGCTTTCTTTTTGGTATAGTGTCGTGGAGTGTCATTCATCTTTGTGCCCCTCCAAAAGTAAAAAcc                                                                                                                                              <  2:3557570/147‑1 (MQ=255)
        gCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACccg                                                                                                                                            <  2:1492648/148‑1 (MQ=255)
         cATGAATAAAGCGCACTCAACAATCTGAAACCCGGCTGAACGGTTTATTACGTACATCAGGTAAAACTGACCTATAAGCCGCTTTCTTTTGGGTATAGTTTCTTGGACTGTAATTCATCTTTCTGCCCCTCCAAAAGTAAAAACccgc                                                                                                                                           <  1:1256355/148‑1 (MQ=21)
         cATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACccgc                                                                                                                                           <  2:2848487/148‑1 (MQ=255)
         cATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACccgc                                                                                                                                           >  2:1388635/1‑148 (MQ=255)
         cATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACccgc                                                                                                                                           <  1:975743/148‑1 (MQ=255)
            gAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCGTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGa                                                                                                                                        >  2:1914541/1‑148 (MQ=255)
            gAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTATTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGa                                                                                                                                        <  2:1045784/148‑1 (MQ=255)
             aaTAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGaa                                                                                                                                       <  1:214476/148‑1 (MQ=255)
             aaTAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCaaa                                                                                                                                                  >  1:1276673/1‑137 (MQ=255)
             aaTAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATATTTCTGCCCCTCCAAAAGTAAAAACCCGCCGaa                                                                                                                                       >  1:114001/1‑148 (MQ=255)
              aTAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAg                                                                                                                                      >  1:2652164/1‑148 (MQ=255)
                  aGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCggg                                                                                                                                  <  1:3277942/148‑1 (MQ=255)
                  aGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAAATGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCggg                                                                                                                                  <  2:2089740/148‑1 (MQ=255)
                    cgcACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCCAAAGGAAAAACCCGCCGACGCGGGtg                                                                                                                                >  2:2553410/1‑147 (MQ=255)
                    cgcACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGtt                                                                                                                                <  1:459931/148‑1 (MQ=255)
                      cACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGtttt                                                                                                                              <  2:3004730/148‑1 (MQ=255)
                        cTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTa                                                                                                                            <  1:976046/148‑1 (MQ=255)
                        cTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTa                                                                                                                            <  1:2698349/148‑1 (MQ=255)
                        cTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTa                                                                                                                            <  2:658726/148‑1 (MQ=255)
                        cTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTa                                                                                                                            <  1:3099044/148‑1 (MQ=255)
                        cTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTa                                                                                                                            <  1:1450723/148‑1 (MQ=255)
                          caacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACg                                                                                                                          >  2:1130056/1‑148 (MQ=255)
                            acaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGt                                                                                                                         >  2:1832707/1‑147 (MQ=255)
                             caaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaa                                                                                                                       <  1:1130056/148‑1 (MQ=255)
                               aTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaaa                                                                                                                     <  1:2232768/148‑1 (MQ=255)
                               aTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGGaaaa                                                                                                                     >  2:3058739/1‑148 (MQ=255)
                                tCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAAc                                                                                                                    <  1:394979/148‑1 (MQ=255)
                                tCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATATTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAacac                                                                                                                    >  2:1171860/1‑148 (MQ=255)
                                     aaCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTGTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAAAAGGTg                                                                                                               >  1:1311716/1‑148 (MQ=255)
                                     aaCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATATTTCTTCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTg                                                                                                               >  1:488579/1‑148 (MQ=255)
                                        ccgccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaaa                                                                                                            <  1:719441/148‑1 (MQ=255)
                                        ccgccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaaa                                                                                                            >  2:3306762/1‑148 (MQ=255)
                                        ccgccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaaa                                                                                                            <  1:2250364/148‑1 (MQ=255)
                                           ccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTg                                                                                                         >  2:1041238/1‑148 (MQ=255)
                                             ggAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAc                                                                                                       <  2:2167778/148‑1 (MQ=255)
                                               aaCGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCg                                                                                                     >  2:524383/1‑148 (MQ=255)
                                               aaCGGTTTATTACGTACAGCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGCACAGTCATTCATCTTTCTGCCCCTCCAAACGTAAAAACCCGCCGAAGCGGGTTTTTACGTAACACAGGTGAAACTGACCg                                                                                                     <  2:318465/148‑1 (MQ=255)
                                                 cGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGTGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGAt                                                                                                   <  1:3156929/148‑1 (MQ=255)
                                                  ggTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATa                                                                                                  <  1:2445583/148‑1 (MQ=255)
                                                   gTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAAGTAAAACAGGTGAAACTGACCGCTaa                                                                                                 >  2:3119443/1‑148 (MQ=255)
                                                   gTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCAATCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTTCGTAAAAAAGGTGTAAATGACCGAGaa                                                                                                 >  2:2418047/1‑148 (MQ=14)
                                                     ttattaCGTACCTCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGc                                                                                               <  1:2346705/148‑1 (MQ=255)
                                                     ttattaCGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGc                                                                                               <  1:3475599/148‑1 (MQ=255)
                                                     ttattaCGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGc                                                                                               <  1:2881074/148‑1 (MQ=255)
                                                     ttattaCGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCCCCCAAGCGGGGTTTTACGTAAAACAGGTGAAACTGACCGATAAGc                                                                                               >  1:148421/1‑148 (MQ=255)
                                                       attaCGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCg                                                                                             <  2:401634/148‑1 (MQ=255)
                                                        ttaCGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGc                                                                                            <  1:658832/148‑1 (MQ=255)
                                                         taCGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGct                                                                                           <  2:1150951/148‑1 (MQ=255)
                                                         taCGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttt                                                                                                                               <  1:1474246/112‑1 (MQ=25)
                                                         taCGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCCGGTGAAACTGACCGATAAGCCGct                                                                                           >  2:3240888/1‑148 (MQ=12)
                                                           cGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGcttt                                                                                         <  2:163560/148‑1 (MQ=255)
                                                               cATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttcttt                                                                                     <  1:1818511/148‑1 (MQ=39)
                                                                  cAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                  <  2:1572397/148‑1 (MQ=37)
                                                                  cAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                  <  2:225291/148‑1 (MQ=37)
                                                                  cAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                  <  2:3340197/148‑1 (MQ=37)
                                                                  cAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                  <  1:3047945/148‑1 (MQ=37)
                                                                    ggTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGttt                                                                                                                               >  2:2014493/1‑101 (MQ=35)
                                                                       aaaaCTGACCGATAAGCCTCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGtata                                                                             <  2:1316296/148‑1 (MQ=25)
                                                                       aaaaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGtata                                                                             <  2:1192197/148‑1 (MQ=35)
                                                                        aaaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAg                                                                            <  2:1607863/148‑1 (MQ=34)
                                                                           cTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt                                                                         >  2:2850327/1‑148 (MQ=37)
                                                                           cTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt                                                                         <  1:516551/148‑1 (MQ=37)
                                                                                  aTAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGAc                                                                  <  2:1117095/148‑1 (MQ=37)
                                                                                  aTAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCATTCTGACCCAACAAAATTAAAAAACCGCCGAAGCGGGTTTTTACGTAAAAAAGGTTAAAATGACCGCAAATACGCTTTCTTTTGGGGATAGTGTCGttgcg                                                                  >  1:2013353/1‑144 (MQ=11)
                                                                                      gCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTAATATTTCAGCCCCGCCAAAAAGAAAAAACCCCCGAAGCCGGGTTTTTCGTAAAACAGGTGAAACAGACCGATAAGCCGCTTTGTGGTGGGGATTGTGTCTTGGACCGtg                                                              >  1:1099191/1‑147 (MQ=2)
                                                                                        cGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcat                                                            <  1:326059/148‑1 (MQ=37)
                                                                                          ctttctttTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcattc                                                          <  1:1388635/148‑1 (MQ=37)
                                                                                            ttctttTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcattcat                                                        <  2:2219832/148‑1 (MQ=37)
                                                                                             tctttTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcattcatt                                                       <  2:2652164/148‑2 (MQ=37)
                                                                                              ctttTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcattcattt                                                      <  1:67864/148‑3 (MQ=37)
                                                                                              ctttTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcattcattt                                                      <  1:1867131/148‑3 (MQ=37)
                                                                                                 ttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTc                                                   <  2:766537/148‑1 (MQ=37)
                                                                                                 ttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTc                                                   <  2:765609/148‑1 (MQ=37)
                                                                                                      tataGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAg                                                                                                <  1:589770/98‑1 (MQ=37)
                                                                                                      tataGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAg                                                                                                >  2:589770/1‑98 (MQ=37)
                                                                                                      tataGTGTCGTGGACAGTCATTCATATTTCTGCCCCTCCAAAAGTAAAAACCCGCCCAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCGTTATTTTGGGTATAGTGTCGTGGACAGTGAATAAATTTTCTGccc                                              >  1:253135/1‑148 (MQ=11)
                                                                                                         aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTc                                           <  2:804238/148‑1 (MQ=37)
                                                                                                         aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTc                                           <  2:972112/148‑1 (MQ=37)
                                                                                                         aGTGTCGTGGACAGTCAATCATCTTTATGCCCCTCCAAAAGTAAAAACCCCCCGAACCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCGTTCGTTTGGGTATAGTGTGGTGGACAGTCCTTCATTTTTGGGCCCCTc                                           >  1:3478105/1‑148 (MQ=11)
                                                                                                                tGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGc                                    <  1:619686/148‑1 (MQ=37)
                                                                                                                     aGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGGCCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCaaaaa                               <  1:1914541/148‑1 (MQ=255)
                                                                                                                      gTCATTCATCTTTCGGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAAc                              <  1:2418047/148‑1 (MQ=255)
                                                                                                                       tcattcatCTTTCTGCCCCTCCAAAATTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAAcc                             <  1:3058756/148‑1 (MQ=255)
                                                                                                                        cattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTg                                                                                                               >  2:3065911/1‑65 (MQ=40)
                                                                                                                           tcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                                                                                                                      >  2:683159/1‑55 (MQ=35)
                                                                                                                             atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGa                       <  2:148421/148‑1 (MQ=37)
                                                                                                                              tCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGaa                      >  1:2565004/1‑148 (MQ=37)
                                                                                                                                           ccAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACCGTCATTCATTTTTCTGCCCCTCCAAAAGCAAATACCCGCCGCCGCGGGTTTTTACg         >  2:2269427/1‑148 (MQ=11)
                                                                                                                                               aaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaaa     <  2:481939/148‑1 (MQ=37)
                                                                                                                                                  tAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGGCCGGCATTCATTTTTCTGCCCCGCCAAAAGCAAAAACCCGCCGCAGCGGGTGTTTACGTAAACCa  >  1:2460773/1‑148 (MQ=11)
                                                                                                                                                  |                                                                                                                                                   
CGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCA  >  CP009273/3263169‑3263462

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: