Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,389,852 G→T 91.8% R110S (CGC→AGC)  ynaI ← mechanosensitive channel protein, very small conductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,389,8520GT91.8% 613.2 / 47.5 231R110S (CGC→AGC) ynaImechanosensitive channel protein, very small conductance
Reads supporting (aligned to +/- strand):  ref base G (12/7);  new base T (100/112);  total (112/119)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 233 total aligned reads displayed.

CCATACCGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCGACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGTTTATAGTTTTCTATAACATAATT  >  CP009273/1389705‑1389996
                                                                                                                                                   |                                                                                                                                                
ccATACCGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATcgcg                                                                                                                                                  <  2:3176885/148‑1 (MQ=255)
  aTACCGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTAc                                                                                                                                                <  1:2243544/148‑1 (MQ=255)
  aTACCGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTAc                                                                                                                                                <  1:2177854/148‑1 (MQ=255)
  aTACCGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTAc                                                                                                                                                <  2:2459512/148‑1 (MQ=255)
   tACCGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTAcc                                                                                                                                               >  1:88228/1‑148 (MQ=255)
   tACCGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCGAcc                                                                                                                                               >  1:211369/1‑148 (MQ=255)
      cGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCttt                                                                                                                                            >  1:489765/1‑148 (MQ=255)
          aGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTa                                                                                                                                        >  2:2709185/1‑148 (MQ=255)
           gCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTaa                                                                                                                                       >  2:1267275/1‑148 (MQ=255)
               gACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATctg                                                                                                                                   >  1:2984654/1‑148 (MQ=255)
                accaccAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATctgc                                                                                                                                  >  1:3134081/1‑148 (MQ=255)
                accaccAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATctgc                                                                                                                                  >  1:2454623/1‑148 (MQ=255)
                   accaATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGt                                                                                                                               <  2:2939132/148‑1 (MQ=255)
                   accaATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGt                                                                                                                               >  2:1534724/1‑148 (MQ=255)
                    ccaATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGtt                                                                                                                              <  1:2397058/148‑1 (MQ=255)
                    ccaATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGtt                                                                                                                              >  1:784330/1‑148 (MQ=255)
                     caATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGttt                                                                                                                             <  1:1167649/148‑1 (MQ=255)
                     caATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCGACCTTTTTTAATCTGCTGttt                                                                                                                             >  1:3115016/1‑148 (MQ=255)
                       aTACCACCAACGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGttttt                                                                                                                           <  2:2783175/148‑1 (MQ=255)
                       aTACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCGACCTTTTTTAATCTGCTGttttt                                                                                                                           >  2:2599673/1‑148 (MQ=255)
                           caccaAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCaaa                                                                                                                       <  2:231967/148‑1 (MQ=255)
                              caAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTg                                                                                                                    <  2:2759005/148‑1 (MQ=255)
                               aaaGGTCAGCAAGCCAGAAAGGCTCCTGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTTCCTTTTTTAATCTGCTGTTTTTCAAACTGa                                                                                                                   <  1:2414902/148‑1 (MQ=255)
                                  ggTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTa                                                                                                                <  2:1991301/148‑1 (MQ=255)
                                  ggTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTa                                                                                                                <  1:2379447/148‑1 (MQ=255)
                                  ggTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTa                                                                                                                >  2:2158652/1‑148 (MQ=255)
                                  ggTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTa                                                                                                                <  1:1749782/148‑1 (MQ=255)
                                  ggTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTa                                                                                                                <  1:2911684/148‑1 (MQ=255)
                                     cAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATc                                                                                                             >  1:612661/1‑148 (MQ=255)
                                       gCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATcag                                                                                                           >  1:960867/1‑148 (MQ=255)
                                        cAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATcagc                                                                                                          <  2:2286301/148‑1 (MQ=255)
                                        cAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATcagc                                                                                                          >  2:3080493/1‑148 (MQ=255)
                                        cAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATcagc                                                                                                          <  2:2917540/148‑1 (MQ=255)
                                         aaGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATcagca                                                                                                         <  1:1577407/148‑1 (MQ=255)
                                          aGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATcagcag                                                                                                        >  2:2818202/1‑148 (MQ=255)
                                           gCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGa                                                                                                       >  1:1020322/1‑148 (MQ=255)
                                           gCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGa                                                                                                       >  1:2143519/1‑148 (MQ=255)
                                           gCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGa                                                                                                       >  1:2531767/1‑148 (MQ=255)
                                            ccAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGaa                                                                                                      >  1:2025393/1‑148 (MQ=255)
                                              aGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATACTGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAAt                                                                                                    <  2:3182901/148‑1 (MQ=255)
                                                 aaGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGc                                                                                                 >  2:2886154/1‑148 (MQ=255)
                                                 aaGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGc                                                                                                 <  2:1088395/148‑1 (MQ=255)
                                                   ggCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCtt                                                                                               >  2:1587056/1‑148 (MQ=255)
                                                    gCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTa                                                                                              >  2:2128963/1‑148 (MQ=255)
                                                     cTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTaa                                                                                             <  2:693320/148‑1 (MQ=255)
                                                     cTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTaa                                                                                             <  2:725728/148‑1 (MQ=255)
                                                      tCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAAt                                                                                            <  1:329503/148‑1 (MQ=255)
                                                      tCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAAt                                                                                            <  1:352946/148‑1 (MQ=255)
                                                         tGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGt                                                                                         <  1:2466304/148‑1 (MQ=255)
                                                         tGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGt                                                                                         <  1:1331341/148‑1 (MQ=255)
                                                                aTGTTCGCCATAAATCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATa                                                                                  <  1:1473627/148‑1 (MQ=255)
                                                                 tGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATaa                                                                                 <  1:2231780/148‑1 (MQ=255)
                                                                   ttCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGNTATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAAcc                                                                               >  1:1677401/1‑148 (MQ=255)
                                                                    tCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCg                                                                              >  2:2664263/1‑148 (MQ=255)
                                                                    tCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCGACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCg                                                                              >  1:919679/1‑148 (MQ=255)
                                                                        cATAAAGCAGAACAAGCACCACAATAATGGTGATTTTTATGATACGCGACATAATCCGCGCACTTGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAAt                                                                          <  2:889453/148‑1 (MQ=255)
                                                                          tAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCa                                                                        >  2:1501328/1‑148 (MQ=255)
                                                                          tAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCa                                                                        >  2:3027308/1‑148 (MQ=255)
                                                                             aGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAggg                                                                     <  1:903028/148‑1 (MQ=255)
                                                                                  aaCAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTa                                                                >  1:2631571/1‑148 (MQ=255)
                                                                                    cAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACa                                                              >  1:138660/1‑148 (MQ=255)
                                                                                    cAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCGACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACa                                                              >  2:2551138/1‑148 (MQ=255)
                                                                                       gCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGAt                                                           >  1:197491/1‑148 (MQ=255)
                                                                                        caccacAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATa                                                          >  1:1834628/1‑148 (MQ=255)
                                                                                          ccacAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATaaa                                                        <  1:1702132/148‑1 (MQ=255)
                                                                                          ccacAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCGACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATaaa                                                        <  2:3020207/148‑1 (MQ=255)
                                                                                           cacaATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATaaaa                                                       >  1:715799/1‑148 (MQ=255)
                                                                                             caATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATc                                                     >  2:1749162/1‑148 (MQ=255)
                                                                                              aataatGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCg                                                    >  2:1378352/1‑148 (MQ=255)
                                                                                               ataatGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGa                                                   <  2:2984654/148‑1 (MQ=255)
                                                                                                 actGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATa                                                 <  2:1825795/146‑1 (MQ=255)
                                                                                                    ggTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACg                                              <  1:2327534/148‑1 (MQ=255)
                                                                                                     gTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGc                                             <  2:492756/148‑1 (MQ=255)
                                                                                                       gATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCt                                           <  1:2664263/148‑1 (MQ=255)
                                                                                                        aTTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACAGCAATCAGGGACGTACAGATAAAATCGATAACGCCTg                                          >  2:2355748/1‑148 (MQ=255)
                                                                                                        aTTTTGATGATACGCGAAATAATCCGCGCACTGGTTATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTg                                          <  1:745798/148‑1 (MQ=255)
                                                                                                         ttttGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTgg                                         >  1:3085901/1‑148 (MQ=255)
                                                                                                           ttGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTgggg                                       <  2:3256291/148‑1 (MQ=255)
                                                                                                              atgatACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGa                                    <  2:1713833/148‑1 (MQ=255)
                                                                                                               tgatACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGAt                                   >  2:2464698/1‑148 (MQ=255)
                                                                                                               tgatACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGAt                                   >  2:42571/1‑148 (MQ=255)
                                                                                                                       gACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGa                           <  2:349312/148‑1 (MQ=255)
                                                                                                                         cATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGt                         >  2:748080/1‑148 (MQ=255)
                                                                                                                         cATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGt                         <  1:1105108/148‑1 (MQ=255)
                                                                                                                           tAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTATTCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGttt                       >  2:2407722/1‑148 (MQ=255)
                                                                                                                           tAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGttt                       <  1:535212/148‑1 (MQ=255)
                                                                                                                                  cgcACTGGTGCTCTCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGTTTATAGttt                >  1:841669/1‑148 (MQ=255)
                                                                                                                                    cACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGTTTATAGTTTTc              <  1:3246460/148‑1 (MQ=255)
                                                                                                                                     aCTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGTTTATAGTTTTCt             >  1:1263552/1‑148 (MQ=255)
                                                                                                                                      ctgtttatatCGATACATTTTTTAATCTGATTTTTTTCAAACTTATTAATCAGCATAAATAGCTTAATCGACATAATAACCGCAAACAGGGACGTACATATAAAATCGATAACGCCTGGGGTGATAAATTTGAGTTTATAGTTTTCta            <  1:561860/142‑1 (MQ=255)
                                                                                                                                      cTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGTTTATAGTTTTCta            >  2:1757547/1‑148 (MQ=255)
                                                                                                                                        ggTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGTTTATAGTTTTCtata          <  2:189193/148‑1 (MQ=255)
                                                                                                                                        ggTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGTTTATAGTTTTCtata          <  2:679493/148‑1 (MQ=255)
                                                                                                                                           gttatCGCTACCTTTTTTAATGTGATTTTTTTAAAACTTATTATTCAGCATAAATAGCTTAAACGTCATAATACCCGCATTCTGGGAAGTACAGTTAATATCGATAACGCCTGGGGTGATAAATTTGAGTTTTTAGTTTTCTATAACa       <  1:797919/146‑1 (MQ=255)
                                                                                                                                            atatCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGTTTATAGTTTTCTATAACAt      <  2:138904/148‑1 (MQ=255)
                                                                                                                                            atatCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGTTTATAGTTTTCTATAACAt      >  1:2759080/1‑148 (MQ=255)
                                                                                                                                             tatCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGTTTATAGTTTTCTATAACAta     >  2:1886080/1‑148 (MQ=255)
                                                                                                                                              atCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGTTTATAGTTTTCTATAACAtaa    >  2:2569931/1‑148 (MQ=255)
                                                                                                                                              atCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGTTTATAGTTTTCTATAACAtaa    >  2:3202407/1‑148 (MQ=255)
                                                                                                                                               tCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGTTTATAGTTTTCTATAACAtaat   <  2:1044670/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CCATACCGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCGACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGTTTATAGTTTTCTATAACATAATT  >  CP009273/1389705‑1389996

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