Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,459,790 C→T 5.5% pseudogene (142/2513 nt) ydbA → pseudogene, autotransporter homolog; interrupted by IS2 and IS30

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,459,7900CT 5.5% 1064.6 / 29.3 290pseudogene (142/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base C (170/104);  new base T (6/10);  total (176/114)
Fisher's exact test for biased strand distribution p-value = 6.53e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 293 total aligned reads displayed.

ATTTATTATGCAAAGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCG  >  CP009273/1459642‑1459938
                                                                                                                                                    |                                                                                                                                                    
 tttATTATGCAAAGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACt                                                                                                                                                     >  1:2351371/1‑149 (MQ=255)
      tatGCAAAGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTc                                                                                                                                                <  2:3837723/149‑1 (MQ=255)
      tatGCAAAGGAAAACTCTATTGTCGGCCTGTAGTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTc                                                                                                                                                <  2:49685/149‑1 (MQ=255)
        tGCAAAGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCg                                                                                                                                              <  1:3888695/149‑1 (MQ=255)
          cAAAGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGa                                                                                                                                            <  1:4830044/149‑1 (MQ=255)
            aaGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAag                                                                                                                                          >  1:1894222/1‑149 (MQ=255)
            aaGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAag                                                                                                                                          >  2:482779/1‑149 (MQ=255)
             aGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAaga                                                                                                                                         >  1:1621255/1‑149 (MQ=255)
               gAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGc                                                                                                                                       <  1:770168/149‑1 (MQ=255)
                  aaCTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTa                                                                                                                                    <  2:3992158/149‑1 (MQ=255)
                  aaCTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTa                                                                                                                                    >  1:1004267/1‑149 (MQ=255)
                     tctATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAAc                                                                                                                                 >  1:4340728/1‑149 (MQ=255)
                       tATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGc                                                                                                                               >  1:479162/1‑149 (MQ=255)
                         ttGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCtt                                                                                                                             >  1:397125/1‑149 (MQ=255)
                          tGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTa                                                                                                                            <  1:4075428/149‑1 (MQ=255)
                          tGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTa                                                                                                                            <  1:3737919/149‑1 (MQ=255)
                           gTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTAc                                                                                                                           <  2:4623033/149‑1 (MQ=255)
                            tCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTAcc                                                                                                                          <  2:282890/149‑1 (MQ=255)
                             cGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTAccc                                                                                                                         >  1:839867/1‑149 (MQ=255)
                              ggCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCct                                                                                                                        <  1:3253106/149‑1 (MQ=255)
                               gCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCctc                                                                                                                       >  1:5069024/1‑149 (MQ=255)
                               gCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCctc                                                                                                                       >  2:720631/1‑149 (MQ=255)
                                ccTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCctct                                                                                                                      <  1:602120/149‑1 (MQ=255)
                                  tGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGa                                                                                                                    >  1:3344923/1‑149 (MQ=255)
                                  tGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGa                                                                                                                    <  1:1085664/149‑1 (MQ=255)
                                    tATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAAt                                                                                                                  <  1:963165/149‑1 (MQ=255)
                                    tATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAAt                                                                                                                  >  1:5138016/1‑149 (MQ=255)
                                     aTTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATg                                                                                                                 >  2:1297559/1‑149 (MQ=255)
                                      ttGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGc                                                                                                                >  1:3243354/1‑149 (MQ=255)
                                        gCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCAGGAAGAGCTAAAACGCTTACCCTCTGAATGctc                                                                                                              >  1:5194/1‑149 (MQ=255)
                                         cATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGctct                                                                                                             >  1:4625987/1‑149 (MQ=255)
                                          aTTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCtctc                                                                                                            >  1:1574160/1‑149 (MQ=255)
                                            tAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCt                                                                                                          <  2:4894229/149‑1 (MQ=255)
                                               ctctGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTa                                                                                                       <  2:3340080/149‑1 (MQ=255)
                                                tctGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAg                                                                                                      >  1:1812822/1‑149 (MQ=255)
                                                   gAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTtcgtcg                                                                                                   >  2:3796961/1‑149 (MQ=255)
                                                   gAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCg                                                                                                   >  1:1370867/1‑149 (MQ=255)
                                                          cAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACaaaa                                                                                            <  1:2623984/149‑1 (MQ=255)
                                                           aGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGACTGCTCTCCTTTAGTCGAACAAAAc                                                                                           >  1:1908150/1‑149 (MQ=255)
                                                           aGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAAc                                                                                           <  1:4243843/149‑1 (MQ=255)
                                                            gggTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAATTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAAcc                                                                                          <  1:4749088/149‑1 (MQ=255)
                                                            gggTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAAcc                                                                                          >  1:5167240/1‑149 (MQ=255)
                                                             ggTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCt                                                                                         <  1:1770848/149‑1 (MQ=255)
                                                                tGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGAt                                                                                      <  1:4583515/149‑1 (MQ=255)
                                                                  ggcggcAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGc                                                                                    <  2:3221725/149‑1 (MQ=255)
                                                                    cggcAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGATAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCa                                                                                  <  2:818221/149‑1 (MQ=255)
                                                                     ggcAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCAt                                                                                 >  2:4564696/1‑149 (MQ=255)
                                                                       cAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATgg                                                                               >  1:3612562/1‑149 (MQ=255)
                                                                            atCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTc                                                                          >  2:2575287/1‑149 (MQ=255)
                                                                                caGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCcaca                                                                      <  1:4983023/149‑1 (MQ=255)
                                                                                  gagGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAgg                                                                    >  1:2990357/1‑149 (MQ=255)
                                                                                     gTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGgcgc                                                                 <  1:286452/149‑1 (MQ=255)
                                                                                        gAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCgctgc                                                              >  2:5141611/1‑149 (MQ=255)
                                                                                        gAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCgctgc                                                              >  2:4080618/1‑149 (MQ=255)
                                                                                         aaaCCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCgctgct                                                             >  2:2131154/1‑149 (MQ=255)
                                                                                         aaaCCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAAACTGATGCCATGGCTTTCCACAGGCgctgct                                                             >  1:4767136/1‑149 (MQ=255)
                                                                                          aaCCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGctgctg                                                            >  2:1232946/1‑149 (MQ=255)
                                                                                             caccacAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGt                                                         >  1:1983558/1‑149 (MQ=255)
                                                                                             caccacAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGt                                                         >  1:3111604/1‑149 (MQ=255)
                                                                                                  caGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATc                                                    >  2:4908538/1‑149 (MQ=255)
                                                                                                   aGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCa                                                   <  1:2432163/149‑1 (MQ=255)
                                                                                                   aGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCa                                                   <  1:3916230/149‑1 (MQ=255)
                                                                                                   aGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCa                                                   <  1:4167743/149‑1 (MQ=255)
                                                                                                       gaaaagaaaaATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGc                                               <  2:3412611/149‑1 (MQ=255)
                                                                                                       gaaaagaaaaATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGc                                               <  2:271093/149‑1 (MQ=255)
                                                                                                       gaaaagaaaaATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGc                                               <  1:1237256/149‑1 (MQ=255)
                                                                                                          aagaaaaATACCAATGTGACTTGTCCGGCAGACCCAGGAAAATTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCtt                                            <  1:3556975/149‑1 (MQ=255)
                                                                                                           agaaaaATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTa                                           >  2:15465/1‑149 (MQ=255)
                                                                                                             aaaaaTACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGc                                         <  2:1621255/149‑1 (MQ=255)
                                                                                                              aaaaTACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGcc                                        <  2:4404322/149‑1 (MQ=255)
                                                                                                              aaaaTACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGcc                                        >  2:430833/1‑149 (MQ=255)
                                                                                                              aaaaTACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGcc                                        >  2:3054786/1‑149 (MQ=255)
                                                                                                                aaTACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGt                                      >  2:916282/1‑149 (MQ=255)
                                                                                                                aaTACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGt                                      <  2:240971/149‑1 (MQ=255)
                                                                                                                  tACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAg                                    >  1:3584800/1‑149 (MQ=255)
                                                                                                                  tACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAg                                    >  2:4128894/1‑149 (MQ=255)
                                                                                                                   aCCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGt                                   <  2:4144584/149‑1 (MQ=255)
                                                                                                                     cAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTgg                                 >  2:3564423/1‑149 (MQ=255)
                                                                                                                         gtgACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACt                             >  1:3005750/1‑149 (MQ=255)
                                                                                                                         gtgACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACt                             >  2:647154/1‑149 (MQ=255)
                                                                                                                             cTTGACCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTAGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAac                         >  2:202094/1‑149 (MQ=255)
                                                                                                                              ttGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAacg                        >  1:2295558/1‑149 (MQ=255)
                                                                                                                                gTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAacgac                      <  2:1812824/149‑1 (MQ=255)
                                                                                                                                 tCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAacgacg                     >  2:2729910/1‑149 (MQ=255)
                                                                                                                                   cGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACgat                   >  1:3502420/1‑149 (MQ=255)
                                                                                                                                    ggCAGACCCAGGAAAATTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACgatg                  >  2:769124/1‑149 (MQ=255)
                                                                                                                                     gCAGACCCAGGAAAATTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACgatga                 >  1:2089561/1‑149 (MQ=255)
                                                                                                                                     gCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACgatga                 >  1:2554212/1‑149 (MQ=255)
                                                                                                                                     gCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACgatga                 >  2:2551830/1‑149 (MQ=255)
                                                                                                                                        gACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACgatgatga              >  2:1411266/1‑149 (MQ=255)
                                                                                                                                          cccAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGatc            >  2:4257917/1‑149 (MQ=255)
                                                                                                                                           ccAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGatca           >  1:4907487/1‑149 (MQ=255)
                                                                                                                                            cAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTCGCCGTAGTGGAACTAAACGACGATGATGatcat          <  2:2990170/149‑1 (MQ=255)
                                                                                                                                             aGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGatcatc         >  2:1783598/1‑149 (MQ=255)
                                                                                                                                             aGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGatcatc         >  1:1640162/1‑149 (MQ=255)
                                                                                                                                               gAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTATCCGTAGTTGAACTAAACGACGATGATGatcatcat       >  2:2155533/1‑149 (MQ=255)
                                                                                                                                               gAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGatcatcat       >  2:1135564/1‑149 (MQ=255)
                                                                                                                                                aaaaCTCAGTCCGTAAGAGCTTAAACGCTTACCCTCTTAATGCTCTCCTTTAGTCTAAAAAAACCTGATGCCTTGGCTTTCCACAGGCGCTGCTGCTTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGatcatcatc      <  1:2886478/149‑1 (MQ=255)
                                                                                                                                                aaaaCTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGatcatcatc      >  1:5043453/1‑149 (MQ=255)
                                                                                                                                                   aCTCTGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACGTGATGCCTTGGCTTTCCACCGGCGCTGCGGCGTTAATCCCGGCCTTAGCCGTTGTGGAACTACACGACGATGATGatcatcatcatc   <  1:4520532/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATTTATTATGCAAAGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCG  >  CP009273/1459642‑1459938

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: