Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,901,514 C→A 39.6% P504P (CCG→CCT)  pyrG ← CTP synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,901,5140CA39.6% 233.6 / 198.8 194P504P (CCG→CCT) pyrGCTP synthetase
Reads supporting (aligned to +/- strand):  ref base C (62/55);  new base A (37/40);  total (99/95)
Fisher's exact test for biased strand distribution p-value = 5.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 196 total aligned reads displayed.

TTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTG  >  CP009273/2901367‑2901660
                                                                                                                                                   |                                                                                                                                                  
ttGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTc                                                                                                                                                    >  2:1071063/1‑148 (MQ=255)
 tGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCg                                                                                                                                                   >  2:493294/1‑148 (MQ=255)
   ccGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGa                                                                                                                                                 >  2:3139603/1‑148 (MQ=255)
   ccGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGa                                                                                                                                                 <  1:2524017/148‑1 (MQ=255)
       tctAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACtt                                                                                                                                             <  1:1044271/148‑1 (MQ=255)
       tctAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACtt                                                                                                                                             <  2:2191123/148‑1 (MQ=255)
         tAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCg                                                                                                                                           <  2:2951486/148‑1 (MQ=255)
          aaCTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAACGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCga                                                                                                                                          <  2:1351433/148‑1 (MQ=255)
          aaCTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCga                                                                                                                                          >  2:2195240/1‑148 (MQ=255)
             ttttttACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCgatga                                                                                                                                       >  1:1424111/1‑148 (MQ=255)
             ttttttACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCgatga                                                                                                                                       >  2:2048996/1‑148 (MQ=255)
             ttttttACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCgatga                                                                                                                                       >  2:2122085/1‑148 (MQ=255)
                   acttacttCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGa                                                                                                                                 >  2:2428059/1‑148 (MQ=255)
                      tacttCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCa                                                                                                                              >  1:2658338/1‑148 (MQ=255)
                       acttCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCaa                                                                                                                             <  2:973002/148‑1 (MQ=255)
                        cttCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAAc                                                                                                                            >  2:1521489/1‑148 (MQ=255)
                              cTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGatca                                                                                                                      <  2:2994286/148‑1 (MQ=255)
                               tGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGatcat                                                                                                                     <  2:643700/148‑1 (MQ=255)
                                 aCGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATcc                                                                                                                   <  1:986589/148‑1 (MQ=255)
                                  cGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATccc                                                                                                                  >  1:549071/1‑148 (MQ=255)
                                    tttCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCg                                                                                                                >  1:1816018/1‑148 (MQ=255)
                                    tttCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCg                                                                                                                >  1:2796697/1‑148 (MQ=255)
                                     ttCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCgg                                                                                                               >  2:1046071/1‑148 (MQ=255)
                                     ttCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCgg                                                                                                               >  1:3216192/1‑148 (MQ=255)
                                         ggAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACg                                                                                                           >  1:105779/1‑148 (MQ=255)
                                          gAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACgg                                                                                                          <  1:281525/148‑1 (MQ=255)
                                            aCTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGcc                                                                                                        >  1:1168671/1‑148 (MQ=255)
                                            aCTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGcc                                                                                                        <  1:1355592/148‑1 (MQ=255)
                                                  tGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAAc                                                                                                  >  1:2635961/1‑148 (MQ=255)
                                                       gcTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCa                                                                                             >  1:2691376/1‑148 (MQ=255)
                                                       gcTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCa                                                                                             >  2:1471728/1‑148 (MQ=255)
                                                       gcTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCa                                                                                             >  2:1918749/1‑148 (MQ=255)
                                                          ttCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGAc                                                                                          <  1:373480/148‑1 (MQ=255)
                                                             acaAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTg                                                                                       <  2:261589/148‑1 (MQ=255)
                                                              caAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGc                                                                                      >  1:2145371/1‑148 (MQ=255)
                                                              caAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGc                                                                                      >  1:2145356/1‑148 (MQ=255)
                                                               aaaGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCa                                                                                     >  1:796027/1‑148 (MQ=255)
                                                               aaaGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCa                                                                                     >  1:795919/1‑148 (MQ=255)
                                                                      gCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCaa                                                                              <  1:3453947/148‑1 (MQ=255)
                                                                        aaaCAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATc                                                                            >  2:3381762/1‑148 (MQ=255)
                                                                          aCAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTg                                                                          >  2:3311092/1‑148 (MQ=255)
                                                                          aCAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTg                                                                          <  1:1508018/148‑1 (MQ=255)
                                                                          aCAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTg                                                                          >  1:875025/1‑148 (MQ=255)
                                                                          aCAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTg                                                                          >  1:1142386/1‑148 (MQ=255)
                                                                           cAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGt                                                                         <  2:2706530/148‑1 (MQ=255)
                                                                           cAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGt                                                                         <  1:2108366/148‑1 (MQ=255)
                                                                           cAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGt                                                                         >  2:2265411/1‑148 (MQ=255)
                                                                            aGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGtt                                                                        <  1:875445/148‑1 (MQ=255)
                                                                            aGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGtt                                                                        <  2:1760523/148‑1 (MQ=255)
                                                                              cGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTc                                                                      >  1:1475447/1‑148 (MQ=255)
                                                                                   gACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAg                                                                 <  1:725874/148‑1 (MQ=255)
                                                                                     ccATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCa                                                               >  2:2940827/1‑148 (MQ=255)
                                                                                     ccATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCa                                                               >  2:1165839/1‑148 (MQ=255)
                                                                                      cATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCAt                                                              >  1:2249585/1‑148 (MQ=255)
                                                                                         cACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCAtgtt                                                           >  2:2311782/1‑148 (MQ=255)
                                                                                         cACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCAtgtt                                                           <  1:2008161/148‑1 (MQ=255)
                                                                                         cACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCAtgtt                                                           <  2:3298556/148‑1 (MQ=255)
                                                                                          aCGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCAtgttg                                                          <  1:3054640/148‑1 (MQ=255)
                                                                                          aCGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCAtgttg                                                          <  1:1117875/148‑1 (MQ=255)
                                                                                           cGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCAtgttgt                                                         >  1:1210773/1‑148 (MQ=255)
                                                                                           cGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCAtgttgt                                                         >  2:2705908/1‑148 (MQ=255)
                                                                                              ggAGTAGAAGTAAACTCCGGAGGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGa                                                      <  1:1290508/148‑1 (MQ=255)
                                                                                               gAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGAc                                                     >  1:734541/1‑148 (MQ=255)
                                                                                                     aaGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTa                                               >  1:1917762/1‑148 (MQ=255)
                                                                                                     aaGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTa                                               <  1:667325/148‑1 (MQ=255)
                                                                                                     aaGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTa                                               <  2:767396/148‑1 (MQ=255)
                                                                                                      aGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGAGTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTaa                                              <  1:1828829/148‑1 (MQ=255)
                                                                                                          aaCTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGt                                          <  1:559733/148‑1 (MQ=255)
                                                                                                           aCTCCGGATGGAACTGGCAAGCCACTAACCACGGTTGATTAGGAACTTCTATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTg                                         <  1:1521489/148‑1 (MQ=255)
                                                                                                           aCTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTg                                         <  2:2683887/148‑1 (MQ=255)
                                                                                                           aCTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTg                                         <  2:461877/148‑1 (MQ=255)
                                                                                                           aCTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTg                                         >  2:3012491/1‑148 (MQ=255)
                                                                                                           aCTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTg                                         >  2:1335098/1‑148 (MQ=255)
                                                                                                             tCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGCGAc                                       <  1:2847494/148‑1 (MQ=255)
                                                                                                              ccGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACg                                      <  2:1824559/148‑1 (MQ=255)
                                                                                                                  aTGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGa                                  <  1:2940827/148‑1 (MQ=255)
                                                                                                                  aTGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGa                                  <  1:1165839/148‑1 (MQ=255)
                                                                                                                  aTGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGa                                  <  2:459646/148‑1 (MQ=255)
                                                                                                                  aTGGAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGa                                  <  2:2462815/148‑1 (MQ=255)
                                                                                                                    ggAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACg                                <  1:1824141/148‑1 (MQ=255)
                                                                                                                    ggAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACg                                >  1:2204762/1‑148 (MQ=255)
                                                                                                                    ggAACTGGCAAGCCACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACg                                <  2:4307/148‑1 (MQ=255)
                                                                                                                     gAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGc                               <  2:2145371/148‑1 (MQ=255)
                                                                                                                     gAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGc                               <  2:2145356/148‑1 (MQ=255)
                                                                                                                           gCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTcaaca                         >  1:2771470/1‑148 (MQ=255)
                                                                                                                            cAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTcaacaa                        <  2:680793/148‑1 (MQ=255)
                                                                                                                            cAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTcaacaa                        >  1:2662735/1‑148 (MQ=255)
                                                                                                                                ccACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGt                    >  1:2375942/1‑148 (MQ=255)
                                                                                                                                ccACGAACCACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGt                    >  1:1638138/1‑148 (MQ=255)
                                                                                                                                     aaCCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGgcg               <  2:1206859/148‑1 (MQ=255)
                                                                                                                                      aCCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGgcgc              <  1:828705/148‑1 (MQ=255)
                                                                                                                                        cACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCAt            <  1:3351063/148‑1 (MQ=255)
                                                                                                                                        cACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCAt            >  2:2536965/1‑148 (MQ=255)
                                                                                                                                        cACGGGTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCAt            <  1:1982428/148‑1 (MQ=255)
                                                                                                                                         aCGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCAtt           <  1:3259947/148‑1 (MQ=255)
                                                                                                                                          cGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTg          >  1:2359726/1‑148 (MQ=255)
                                                                                                                                           gggTGATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGt         <  2:886180/148‑1 (MQ=255)
                                                                                                                                              tGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACa      <  2:2033094/148‑1 (MQ=255)
                                                                                                                                                aTTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGc    <  2:1379824/148‑1 (MQ=255)
                                                                                                                                                  tCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTg  <  2:549071/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTG  >  CP009273/2901367‑2901660

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: ‑