Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,389,852 G→T 100% 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,8520GT98.3% 726.7 / 2.9 233R110S (CGC→AGC) ynaImechanosensitive channel protein, very small conductance
Reads supporting (aligned to +/- strand):  ref base G (3/1);  new base T (116/113);  total (119/114)
Fisher's exact test for biased strand distribution p-value = 6.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 234 total aligned reads displayed.

ATACCGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCGACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGTTTATAGTTTTCTATAACATAATT  >  CP009273/1389707‑1389996
                                                                                                                                                 |                                                                                                                                                
aTACCGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTAc                                                                                                                                                <  1:386879/148‑1 (MQ=255)
  aCCGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCt                                                                                                                                              >  2:193403/1‑148 (MQ=255)
     gACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCtttt                                                                                                                                           >  1:1718284/1‑148 (MQ=255)
      aCAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCttttt                                                                                                                                          <  1:258519/148‑1 (MQ=255)
      aCAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCttttt                                                                                                                                          >  1:444691/1‑148 (MQ=255)
       cAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCtttttt                                                                                                                                         >  2:1898186/1‑148 (MQ=255)
        aGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTa                                                                                                                                        <  2:209643/148‑1 (MQ=255)
         gCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTaa                                                                                                                                       <  2:1861383/148‑1 (MQ=255)
          ccAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAAt                                                                                                                                      <  1:2805614/148‑1 (MQ=255)
          ccAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAAt                                                                                                                                      <  1:3205373/148‑1 (MQ=255)
              accaccAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATctgc                                                                                                                                  <  2:653423/148‑1 (MQ=255)
              accaccAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATctgc                                                                                                                                  >  1:613503/1‑148 (MQ=255)
                  ccaATACCACCAAAGGTCAGCACGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGtt                                                                                                                              <  2:3325644/148‑1 (MQ=255)
                   caATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGttt                                                                                                                             <  1:3478257/148‑1 (MQ=255)
                        ccaccaAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCaa                                                                                                                        <  2:959897/148‑1 (MQ=255)
                        ccaccaAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCaa                                                                                                                        >  2:1177532/1‑148 (MQ=255)
                        ccaccaAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCaa                                                                                                                        <  1:1183842/148‑1 (MQ=255)
                        ccaccaAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCaa                                                                                                                        <  1:1181800/148‑1 (MQ=255)
                        ccaccaAAGGTCAGCAAGCCAGAAACGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCaa                                                                                                                        <  2:3442825/148‑1 (MQ=255)
                           ccaAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACt                                                                                                                     >  1:301961/1‑148 (MQ=255)
                             aaaGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGa                                                                                                                   <  1:260759/148‑1 (MQ=255)
                              aaGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGAt                                                                                                                  >  1:3212521/1‑148 (MQ=255)
                              aaGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGAt                                                                                                                  >  1:1441160/1‑148 (MQ=255)
                               aGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGAtt                                                                                                                 >  2:819060/1‑148 (MQ=255)
                                ggTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTa                                                                                                                >  2:3505883/1‑148 (MQ=255)
                                ggTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTa                                                                                                                >  2:736750/1‑148 (MQ=255)
                                ggTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAAATGATTa                                                                                                                >  2:2499279/1‑148 (MQ=255)
                                 gTCAGCAAGCCAGAAGGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTaa                                                                                                               <  1:2174321/148‑1 (MQ=255)
                                   cAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATc                                                                                                             <  2:1515065/148‑1 (MQ=255)
                                   cAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATc                                                                                                             <  2:2871893/148‑1 (MQ=255)
                                      cAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATcagc                                                                                                          >  1:801877/1‑148 (MQ=255)
                                      cAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATcagc                                                                                                          >  1:1359934/1‑148 (MQ=255)
                                            aGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAAt                                                                                                    >  1:893405/1‑148 (MQ=255)
                                             gAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATa                                                                                                   <  2:3075930/148‑1 (MQ=255)
                                              aaaGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAg                                                                                                  >  1:696319/1‑148 (MQ=255)
                                              aaaGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAg                                                                                                  >  1:1258703/1‑148 (MQ=255)
                                              aaaGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAg                                                                                                  >  1:1717969/1‑148 (MQ=255)
                                                agccTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCt                                                                                                >  2:2989527/4‑148 (MQ=255)
                                                aGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCt                                                                                                >  2:550301/1‑148 (MQ=255)
                                                 ggCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCtt                                                                                               <  2:263085/148‑1 (MQ=255)
                                                 ggCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCtt                                                                                               <  2:262536/148‑1 (MQ=255)
                                                  gCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTa                                                                                              <  2:1815701/148‑1 (MQ=255)
                                                         ccGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCa                                                                                       >  1:2471219/1‑148 (MQ=255)
                                                          cGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAg                                                                                      >  1:942428/1‑148 (MQ=255)
                                                           gAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGa                                                                                     >  1:2031952/1‑148 (MQ=255)
                                                            aaaTGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGaa                                                                                    >  1:910988/1‑148 (MQ=255)
                                                             aaTGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAAt                                                                                   >  1:2539628/1‑148 (MQ=255)
                                                                gTTCGCCATAACGCAGAACAAGCACCACAATACTGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTACTCGTCAGAATAAc                                                                                <  2:1062539/148‑1 (MQ=255)
                                                                 ttCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAAcc                                                                               >  1:2933471/1‑148 (MQ=255)
                                                                 ttCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAAcc                                                                               >  1:1010198/1‑148 (MQ=255)
                                                                     ccATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCaa                                                                           >  1:3335906/1‑148 (MQ=255)
                                                                       aTAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATc                                                                         >  1:151361/1‑148 (MQ=255)
                                                                          aaGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAgg                                                                      >  1:1869686/1‑148 (MQ=255)
                                                                            gCAGAACAAGCACCACACTACTGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGa                                                                    <  2:1028632/148‑1 (MQ=255)
                                                                            gCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGa                                                                    >  1:1179807/1‑148 (MQ=255)
                                                                             cAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGAc                                                                   >  2:1807153/1‑148 (MQ=255)
                                                                             cAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGAc                                                                   >  1:1034596/1‑148 (MQ=255)
                                                                              aGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACg                                                                  <  1:1167322/148‑1 (MQ=255)
                                                                              aGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACg                                                                  <  1:3399438/148‑1 (MQ=255)
                                                                                  cAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACa                                                              >  1:2009222/1‑148 (MQ=255)
                                                                                   aaGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCCCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAg                                                             >  2:141259/1‑148 (MQ=255)
                                                                                    aGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGa                                                            >  2:2822638/1‑148 (MQ=255)
                                                                                     gCACCACTATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCGACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGAt                                                           >  2:1293190/1‑148 (MQ=255)
                                                                                     gCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGAt                                                           >  2:2559897/1‑148 (MQ=255)
                                                                                     gCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGAt                                                           >  2:3392436/1‑148 (MQ=255)
                                                                                      caccacAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATa                                                          >  1:632007/1‑148 (MQ=255)
                                                                                         cacaATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATaaaa                                                       <  2:511108/148‑1 (MQ=255)
                                                                                            aataatGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTTAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCg                                                    <  2:2339508/148‑1 (MQ=255)
                                                                                            aataatGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCg                                                    <  2:24118/148‑1 (MQ=255)
                                                                                               aatGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATa                                                 <  1:1679994/148‑1 (MQ=255)
                                                                                               aatGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATa                                                 <  2:1319860/148‑1 (MQ=255)
                                                                                                atGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATaa                                                >  2:463641/1‑148 (MQ=255)
                                                                                                 tGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAAc                                               >  2:266776/1‑148 (MQ=255)
                                                                                                    tGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGcc                                            >  1:876966/1‑148 (MQ=255)
                                                                                                    tGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGcc                                            <  1:2679175/148‑1 (MQ=255)
                                                                                                     gATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCt                                           <  1:771323/148‑1 (MQ=255)
                                                                                                     gATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCt                                           <  2:914857/148‑1 (MQ=255)
                                                                                                        tttGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCAATAACGCCTggg                                        <  2:444691/148‑1 (MQ=255)
                                                                                                         ttGATGATACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTgggg                                       >  1:1201960/1‑148 (MQ=255)
                                                                                                               atACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATaa                                 <  1:3149215/148‑1 (MQ=255)
                                                                                                               atACGCGACATAATCCGCGCACTGGTGATATCGCGACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATaa                                 >  1:3109523/1‑148 (MQ=255)
                                                                                                                tACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATaaa                                <  2:784727/148‑1 (MQ=255)
                                                                                                                tACGCGACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATaaa                                <  2:3284849/148‑1 (MQ=255)
                                                                                                                    cgACATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTg                            >  1:2277057/1‑148 (MQ=255)
                                                                                                                       cATAATCCGCGCACTGGTGATTTCGTTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTTATCAGCTGAAATAGCTTAATCGTCAGAATAACCGCAATCTGGGACGTACAGATAAAATCGTTAACGCCTGGGGTGATAAATTTGAGt                         <  2:3046411/148‑1 (MQ=255)
                                                                                                                       cATAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGt                         >  1:2965546/1‑148 (MQ=255)
                                                                                                                         tAATCCGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGttt                       >  2:19840/1‑148 (MQ=255)
                                                                                                                           aTCCGCGCACTGGTGATATCGCGACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGTTtat                     <  2:3171099/148‑1 (MQ=255)
                                                                                                                            tCCGCGCACTGGTGAAATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGTTtata                    >  1:2047949/1‑148 (MQ=255)
                                                                                                                             ccGCGCACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATAAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGTTTATAg                   >  2:2959185/1‑148 (MQ=255)
                                                                                                                              cgcgCACGGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGTTTATAGt                  <  1:509737/148‑1 (MQ=255)
                                                                                                                                  cACTGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGTTTATAGTTTTc              <  2:62557/148‑1 (MQ=255)
                                                                                                                                     tGGTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGTTTATAGTTTTCtat           <  2:2463520/148‑1 (MQ=255)
                                                                                                                                      ggTGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGTTTATAGTTTTCtata          <  2:2277057/148‑1 (MQ=255)
                                                                                                                                       gTGATATCGCTACCTTTTTTAATCTGCTTTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGTTTATAGTTTTCTATaa         <  1:1797162/148‑1 (MQ=255)
                                                                                                                                        tGATATCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACTCCTGGGGTGATAAATTTGAGTTTATTGTTTTCTATAAc        >  2:1483046/1‑148 (MQ=255)
                                                                                                                                          atatCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGTTTATAGTTTTCTATAACAt      <  1:1917721/148‑1 (MQ=255)
                                                                                                                                           tatCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCTTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTAGAGTTTATAGTTTTCTATAACAta     >  1:2340795/1‑148 (MQ=255)
                                                                                                                                            atCGCTACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGTTTATAGTTTTCTATAACAtaa    <  2:3235497/148‑1 (MQ=255)
                                                                                                                                              cgctACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGTTTATAGTTTTCTATAACAtaatt  <  2:632007/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
ATACCGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCGACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGTTTATAGTTTTCTATAACATAATT  >  CP009273/1389707‑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: