Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,481,791 C→A 12.6% P184P (CCG→CCT idnT ← L‑idonate and D‑gluconate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,481,7910CA12.6% 1020.3 / 75.9 289P184P (CCG→CCTidnTL‑idonate and D‑gluconate transporter
Reads supporting (aligned to +/- strand):  ref base C (127/125);  new base A (18/19);  total (145/144)
Fisher's exact test for biased strand distribution p-value = 8.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 293 total aligned reads displayed.

TGACCGGGATCACGGCAGCGAAAATACTGTTCCAGAAGGAGGGCATCTCCTCTTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGCCCTGGATGTGGCGGCAGAAAACAGTGGGTTACAGAGAGCGCTGCTACCATTGGTACGCCAACATACAACAGGGG  >  CP009273/4481643‑4481940
                                                                                                                                                    |                                                                                                                                                     
 gACCGGGATCACGGCAGCGAAAATACTGTTCCAGAAGGAGGGCATCTCCTCTTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCgg                                                                                                                                                     >  2:4523209/1‑150 (MQ=255)
    cGGGATCACGGCAGCGAAAATACTGTTCCAGAAGGAGGGCATCTCCTCTTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTAGGAAt                                                                                                                                                  >  1:682844/1‑150 (MQ=255)
     gggATCACGGCAGCGAAAATACTGTTCCAGAAGGAGGGCATCTCCTCTTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATg                                                                                                                                                 <  1:1100111/150‑1 (MQ=255)
     gggATCACGGCAGCGAAAATACTGTTCCAGAAGGAGGGCATCTCCTCTTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATg                                                                                                                                                 >  1:2552856/1‑150 (MQ=255)
       gATCACGGCAGCGAAAATACTGTTCCAGAAGGAGGGCATCTCCTCTTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGt                                                                                                                                               >  2:2182148/1‑150 (MQ=255)
        aTCACGGCAGCGAAAATACTGTTCCAGAAGGAGGGCATCTCCTCTTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTa                                                                                                                                              >  1:1588628/1‑150 (MQ=255)
         tCACGGCAGCGAAAATACTGTTCCAGAAGGAGGGCATCTCCTCTTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTaa                                                                                                                                             >  2:1598991/1‑150 (MQ=255)
          cACGGCAGCGAAAATACTGTTCCAGAAGGAGGGCATCTCCTCTTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAAt                                                                                                                                            <  1:2417704/150‑1 (MQ=255)
           aCGGCAGCGAAAATACTGTTCCAGAAGGAGGGCATCTCCTCTTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATg                                                                                                                                           <  1:281984/150‑1 (MQ=255)
             ggCAGCGAAAATACTGTTCCAGAAGGAGGGCATCTCCTCTTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGAt                                                                                                                                         >  1:4402967/1‑150 (MQ=255)
             ggCAGCGAAAATACTGTTCCAGAAGGAGGGCATCTCCTCTTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGAt                                                                                                                                         <  1:933047/150‑1 (MQ=255)
             ggCAGCGAAAATACTGTTCCAGAAGGAGGGCATCTCCTCTTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGAt                                                                                                                                         >  2:1880107/1‑150 (MQ=255)
             ggCAGCGAAAATACTGTTCCAGAAGGAGGGCATCTCCTCTTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTAGGAATGGTAATGAt                                                                                                                                         >  1:3648053/1‑150 (MQ=255)
              gCAGCGAAAATACTGTTCCAGAAGGAGGGCATCTCCTCTTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATa                                                                                                                                        <  1:4162722/150‑1 (MQ=255)
                aGCGAAAATACTGTTCCAGAAGGAGGGCATCTCCTCTTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTAGGAATGGTAATGATaaa                                                                                                                                      >  2:1685266/1‑150 (MQ=255)
                aGCGAAAATACTGTTCCAGAAGGAGGGCATCTCCTCTTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTAGGAATGGTAATGATaaa                                                                                                                                      >  1:2790220/1‑150 (MQ=255)
                 gCGAAAATACTGTTCCAGAAGGAGGGCATCTCCTCTTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAAt                                                                                                                                     <  1:2151674/150‑1 (MQ=255)
                  cGAAAATACTGTTCCAGAAGGAGGGCATCTCCTCTTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATc                                                                                                                                    <  2:4333743/150‑1 (MQ=255)
                  cGAAAATACTGTTCCAGAAGGAGGGCATCTCCTCTTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATc                                                                                                                                    >  2:294814/1‑150 (MQ=255)
                    aaaaTACTGTTCCAGAAGGAGGGCATCTCCTCTTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCa                                                                                                                                  <  1:735748/150‑1 (MQ=255)
                    aaaaTACTGTTCCAGAAGGAGGGCATCTCCTCTTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTAGGAATGGTAATGATAAATCCa                                                                                                                                  <  2:417771/150‑1 (MQ=255)
                      aaTACTGTTCCAGAAGGAGGGCATCTCCTCTTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATa                                                                                                                                >  1:4455152/1‑150 (MQ=255)
                      aaTACTGTTCCAGAAGGAGGGCATCTCCACTTCGCTAAACAGATGAGTATTAAAAAAGCCTTCCGGTGGTGCTTTCTCAAAGCGATTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATa                                                                                                                                >  1:254291/1‑150 (MQ=255)
                       atatTGTTCCAGAAGGAGGGTATCTCCTCTTCGCTAATCATATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGTGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCTTAc                                                                                                                               <  2:2030533/150‑1 (MQ=255)
                       aTACTGTTCCAGAAGGAGGGCATCTCCTCTTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATAc                                                                                                                               <  2:3721244/150‑1 (MQ=255)
                         aCTGTTCCAGAAGGAGGGCATCTCCTCTTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAg                                                                                                                             <  2:3683877/150‑1 (MQ=255)
                            gTTCCAGAAGGAGGGCATCTCCTCTTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTaa                                                                                                                          <  2:1598497/150‑1 (MQ=255)
                              tCCAGAAGGAGGGCATCTCCTCTTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAg                                                                                                                        <  1:1070335/150‑1 (MQ=255)
                                cAGAAGGAGGGCATCTCCTCTTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTc                                                                                                                      >  2:1380019/1‑150 (MQ=255)
                                cAGAAGGAGGGCATCTCCTCTTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTc                                                                                                                      >  2:1937665/1‑150 (MQ=255)
                                cAGAAGGAGGGCATCTCCTCTTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTAGGAATGGTAATGATAAATCCATACAGTAAAGTc                                                                                                                      <  1:1240027/150‑1 (MQ=255)
                                   aaGGAGGGCATCTCCTCTTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGtt                                                                                                                   <  1:1966559/150‑1 (MQ=255)
                                      gaggGCATCTCCTCTTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCg                                                                                                                >  2:2655539/1‑150 (MQ=255)
                                          gCATCTCCTCTTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGAt                                                                                                            >  1:3034605/1‑150 (MQ=255)
                                           cATCTCCTCTTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGAtt                                                                                                           >  1:1724238/1‑150 (MQ=255)
                                                cctcTTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCt                                                                                                      <  2:2204342/150‑1 (MQ=255)
                                                   ctTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCaa                                                                                                   >  1:4340343/1‑150 (MQ=255)
                                                   ctTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGAAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCTTTCCGAGATTAGACTCaa                                                                                                   >  1:3933715/1‑150 (MQ=255)
                                                     tCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAg                                                                                                 >  2:3762898/1‑150 (MQ=255)
                                                       gCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGAt                                                                                               >  1:4426691/1‑150 (MQ=255)
                                                       gCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTAGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGAt                                                                                               <  2:200312/150‑1 (MQ=255)
                                                        cTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATa                                                                                              >  2:3400278/1‑150 (MQ=255)
                                                        cTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTAGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATa                                                                                              >  1:87419/1‑150 (MQ=255)
                                                         tAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAg                                                                                             >  1:4185691/1‑150 (MQ=255)
                                                          aaaCAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGt                                                                                            >  2:862325/1‑150 (MQ=255)
                                                          aaaCAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGt                                                                                            >  2:862145/1‑150 (MQ=255)
                                                           aaCAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTc                                                                                           <  2:4366315/150‑1 (MQ=255)
                                                            aCAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTcg                                                                                          >  2:73432/1‑150 (MQ=255)
                                                            aCAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTcg                                                                                          >  2:1139431/1‑150 (MQ=255)
                                                            aCAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTcg                                                                                          <  1:1281382/150‑1 (MQ=255)
                                                            aCAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTAGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTcg                                                                                          <  1:1685266/150‑1 (MQ=255)
                                                             cAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTcgc                                                                                         >  1:3492367/1‑150 (MQ=255)
                                                             cAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTcgc                                                                                         >  1:3813435/1‑150 (MQ=255)
                                                               gATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGa                                                                                       >  1:367930/1‑150 (MQ=255)
                                                                atgatgATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGAt                                                                                      <  1:3881605/150‑1 (MQ=255)
                                                                aTGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGAt                                                                                      <  1:1767396/150‑1 (MQ=255)
                                                                   aGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGc                                                                                   <  1:4465224/150‑1 (MQ=255)
                                                                    ggATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCa                                                                                  <  1:4541743/150‑1 (MQ=255)
                                                                     gATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTAGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAg                                                                                 <  2:1749894/150‑1 (MQ=255)
                                                                      aTTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGt                                                                                <  2:413734/150‑1 (MQ=255)
                                                                        tAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAg                                                                              >  2:84995/1‑150 (MQ=255)
                                                                         aaaTAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAgg                                                                             >  2:2650662/1‑150 (MQ=255)
                                                                            tAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGccc                                                                          <  2:3629324/150‑1 (MQ=255)
                                                                             aaGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTAGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGCCCt                                                                         <  1:330720/150‑1 (MQ=255)
                                                                              aGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGCCCTg                                                                        >  2:1224771/1‑150 (MQ=255)
                                                                                ccTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGCCCTGGa                                                                      <  1:228368/150‑1 (MQ=255)
                                                                                     cGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGCCCTGGATGTgg                                                                 <  2:1103381/150‑1 (MQ=255)
                                                                                      ggtggtGCTTTCTCATAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGCCCTGGATGTggc                                                                <  2:2892062/150‑1 (MQ=255)
                                                                                      ggtggtGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGCCCTGGATGTggc                                                                >  2:496918/1‑150 (MQ=255)
                                                                                      ggtggtGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGCCCTGGATGTggc                                                                >  2:4271025/1‑150 (MQ=255)
                                                                                      ggtggtGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGCCCTGGATGTggc                                                                >  2:2667694/1‑150 (MQ=255)
                                                                                      ggtggtGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGCCCTGGATGTggc                                                                <  2:2155629/150‑1 (MQ=255)
                                                                                      ggtggtGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTAGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGCCCTGGATGTggc                                                                >  2:1157067/1‑150 (MQ=255)
                                                                                            gCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGTGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGCCCTGGATGTGGCGGCAGa                                                          >  1:1147817/1‑150 (MQ=255)
                                                                                              tttCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGCCCTGGATGTGGCGGCAGaaa                                                        >  1:1373758/1‑150 (MQ=255)
                                                                                                tctcAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGTTAATTATAAATCAATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGCCCTGGATGTGGCGGCAGAAAAc                                                      >  1:3592218/1‑150 (MQ=255)
                                                                                                     aaGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGCCCTGGATGTGGCGGCAGAAAACAGTgg                                                 >  1:2937781/1‑150 (MQ=255)
                                                                                                      aGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGCCCTGGATGTGGCGGCAGAAAACAGTggg                                                <  2:3942832/150‑1 (MQ=255)
                                                                                                          aGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGCCCTGGATGTGGCGGCAGAAAACAGTGGGTTAc                                            >  2:753808/1‑150 (MQ=255)
                                                                                                          aGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGCCCTGGATGTGGCGGCAGAAAACAGTGGGTTAc                                            >  2:2363954/1‑150 (MQ=255)
                                                                                                           gTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGCCCTGGATGTGGCGGCAGAAAACAGTGGGTTACa                                           >  2:3239369/1‑150 (MQ=255)
                                                                                                           gTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGCCCTGGATGTGGCGGCAGAAAACAGTGGGTTACa                                           <  1:1560081/150‑1 (MQ=255)
                                                                                                             tAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGCCCTGGATGTGGCGGCAGAAAACAGTGGGTTACaga                                         >  1:3490315/1‑150 (MQ=255)
                                                                                                                       gAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGCCCTGGATGTGGCGGCAGAAAACAGTGGGTTACAGAGAGCgctgct                               <  2:1523249/150‑1 (MQ=255)
                                                                                                                        aaaaCAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGCCCTGGATGTGGCGGCAGAAAACAGTGGGTTACAGAGAGCGCTGCTa                              >  1:3994439/1‑150 (MQ=255)
                                                                                                                        aaaaCAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGCCCTGGATGTGGCGGCAGAAAACAGTGGGTTACAGAGAGCGCTGCTa                              >  1:1314284/1‑150 (MQ=255)
                                                                                                                          aactggggTCCTGCTACAATAACTGTCGGAATGGTAATGATAAATCCTTACAGTTAAGTCGTTCAGATATTAGCATCAAAGATAGTCTCGATGGCATTAGGCCATGGATGTTGCGGCAGAAAACATTGGGTTACATAGTTCGCTTCTAcc                            <  2:3766644/144‑1 (MQ=255)
                                                                                                                          aaCAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGCCCTGGATGTGGCGGCAGAAAACAGTGGGTTACAGAGAGCGCTGCTAcc                            <  1:4592968/150‑1 (MQ=255)
                                                                                                                             aGCGGTCCTGCGACAATTACTGTCTGATTTTTAATGATAAATACATACAGTAATTTCGTTACGAGATTAGCCTCAACGATAGTCGCGATGGCAGTAGGCCCTGGATGTGGCGGCAGAAAACAGTGGGTTACAGAGAGCGCTGCTACCAtt                         <  2:254291/150‑1 (MQ=255)
                                                                                                                             aGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGCCCTGGATGTGGCGGCAGAAAACAGTGGGTTACAGAGAGCGCTGCTACCAtt                         <  2:4455146/150‑1 (MQ=255)
                                                                                                                             aGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGCCCTGGATGTGGCGGCAGAAAACAGTGGGTTACAGAGAGCGCTGCTACCAtt                         <  1:1860622/150‑1 (MQ=255)
                                                                                                                               cGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGCCCTGGATGTGGCGGCAGAAAACAGTGGGTTACAGAGAGCGCTGCTACCATTgg                       >  2:2633871/1‑150 (MQ=255)
                                                                                                                                 gTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGCCCTGGATGTGGCGGCAGAAAACAGTGGGTTACAGAGAGCGCTGCTACCATTGGTa                     <  2:3628179/150‑1 (MQ=255)
                                                                                                                                  tCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGCCCTGGATGTGGCGGCAGAAAACAGTGGGTTACAGAGAGCGCTGCTACCATTGGTAc                    <  2:3360950/150‑1 (MQ=255)
                                                                                                                                    cTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGCCCTGGATGTGGCGGCAGAAAACAGTGGGTTACAGAGAGCGCTGCTACCATTGGTACGc                  >  2:3906473/1‑150 (MQ=255)
                                                                                                                                          cAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGCCCTGGATGTGGCGGCAGAAAACAGTGGGTTACAGAGAGCGCTGCTACCATTGGTACGCCAACAt            <  1:143618/150‑1 (MQ=255)
                                                                                                                                             tAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGCCCTGGATGTGGCGGCAGAAAACAGTGGGTTACAGAGAGCGCTGCTACCATTGGTACGCCAACATaca         <  2:1053829/150‑1 (MQ=255)
                                                                                                                                              aaCTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGCCCTGGATGTGGCGGCAGAAAACAGTGGGTTACAGAGAGCGCTGCTACCATTGGTACGCCAACATacaa        <  1:1105668/150‑1 (MQ=255)
                                                                                                                                                cTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGCCCTGGATGTGGCGGCAGAAAACAGTGGGTTACAGAGAGCGCTGCTACCATTGGTACGCCAACATacaaca      <  1:3929822/150‑1 (MQ=255)
                                                                                                                                                    cGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGCCCTGGATGTGGCGGCAGAAAACAGTGGGTTACAGAGAGCGCTGCTACCATTGGTACGCCAACATACAACAgggg  <  1:147154/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGACCGGGATCACGGCAGCGAAAATACTGTTCCAGAAGGAGGGCATCTCCTCTTCGCTAAACAGATGAGGATTAAATAAGCCTTCCGGTGGTGCTTTCTCAAAGCGAGTTAGCAGTTTAGAAAACAGCGGTCCTGCGACAATAACTGTCGGAATGGTAATGATAAATCCATACAGTAAAGTCGTTCCGAGATTAGCCTCAAAGATAGTCGCGATGGCAGTAGGCCCTGGATGTGGCGGCAGAAAACAGTGGGTTACAGAGAGCGCTGCTACCATTGGTACGCCAACATACAACAGGGG  >  CP009273/4481643‑4481940

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: