Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,351,573 C→A 16.2% R60I (AGA→ATA)  cadC ← cadBA operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,351,5730CA16.2% 741.0 / 82.1 245R60I (AGA→ATA) cadCcadBA operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (121/84);  new base A (19/21);  total (140/105)
Fisher's exact test for biased strand distribution p-value = 2.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 254 total aligned reads displayed.

TCGCTGTACCAGATAACCGGCACCATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATGGACGGAGTAACAAGCC  >  CP009273/4351426‑4351720
                                                                                                                                                   |                                                                                                                                                   
  gCTGTACCAGATAACCGGCACCATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTtctc                                                                                                                                                   <  1:1055513/148‑1 (MQ=255)
   cTGTACCAGATAACCGGCACCATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTctct                                                                                                                                                  <  2:2672955/148‑1 (MQ=255)
    tGTACCAGATAACCGGCACCATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCtt                                                                                                                                                 <  2:3226570/148‑1 (MQ=255)
     gTACCAGATAACCGGCACCATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTc                                                                                                                                                >  2:2976971/1‑148 (MQ=255)
      tACCAGATAACCGGCACCATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTATCTTcc                                                                                                                                               <  2:767091/148‑1 (MQ=255)
         cAGATAACCGGCACCATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGa                                                                                                                                            >  1:1767430/1‑148 (MQ=255)
            aTAACCGGCACCATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACAt                                                                                                                                         >  1:1435977/1‑148 (MQ=255)
             tAACCGGCACCATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACAtt                                                                                                                                        <  1:3100140/148‑1 (MQ=255)
             tAACCGGCACCATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACAtt                                                                                                                                        <  1:1827537/148‑1 (MQ=255)
             tAACCGGCACCATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTATCTTCCAGACAtt                                                                                                                                        <  2:1461272/148‑1 (MQ=255)
              aaCCGGCACCATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTATCTTCCAGACATTa                                                                                                                                       <  1:2974145/148‑1 (MQ=255)
               aCCGGCACCATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTAt                                                                                                                                      <  1:1926753/148‑1 (MQ=255)
               aCCGGCACCATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTAt                                                                                                                                      <  1:2274816/148‑1 (MQ=255)
                ccGGCACCATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATc                                                                                                                                     >  1:917341/1‑148 (MQ=255)
                ccGGCACCATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATc                                                                                                                                     >  2:3174850/1‑148 (MQ=255)
                ccGGCACCATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATc                                                                                                                                     >  2:3567441/1‑148 (MQ=255)
                ccGGCACCATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATc                                                                                                                                     >  2:501420/1‑148 (MQ=255)
                 cGGCACCATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCg                                                                                                                                    <  1:1011162/148‑1 (MQ=255)
                      ccATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAg                                                                                                                               >  1:3321008/1‑148 (MQ=255)
                       cATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGt                                                                                                                              >  2:3804600/1‑148 (MQ=255)
                       cATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGt                                                                                                                              >  2:1481594/1‑148 (MQ=255)
                        aTTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGtt                                                                                                                             >  2:774731/1‑148 (MQ=255)
                         ttAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTc                                                                                                                            <  2:1382391/148‑1 (MQ=255)
                              ttATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATccc                                                                                                                       <  2:405181/148‑1 (MQ=255)
                                ataGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTg                                                                                                                     <  2:1052354/148‑1 (MQ=255)
                                ataGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTg                                                                                                                     <  2:3099484/148‑1 (MQ=255)
                                ataGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTg                                                                                                                     >  2:1571391/1‑148 (MQ=255)
                                ataGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTATCTTCCAGACATTATCGATAAGTTCATCCCTg                                                                                                                     >  1:1262450/1‑148 (MQ=255)
                                   gCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTa                                                                                                                  <  2:656212/148‑1 (MQ=255)
                                   gCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTa                                                                                                                  <  1:546711/148‑1 (MQ=255)
                                          tttGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACtt                                                                                                           <  1:1248203/148‑1 (MQ=255)
                                           ttGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTc                                                                                                          >  1:2678368/1‑148 (MQ=255)
                                           ttGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTc                                                                                                          >  1:3435349/1‑148 (MQ=255)
                                                aCAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCAc                                                                                                     >  1:386213/1‑148 (MQ=255)
                                                     aGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTgtgt                                                                                                <  1:2799722/148‑1 (MQ=255)
                                                       cGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTg                                                                                              <  1:3078050/148‑1 (MQ=255)
                                                           ataGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCa                                                                                          >  2:3349769/1‑148 (MQ=255)
                                                             aGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCaaa                                                                                        >  2:3767838/1‑148 (MQ=255)
                                                              gACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAg                                                                                       >  2:2343262/1‑148 (MQ=255)
                                                              gACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAg                                                                                       >  1:2557810/1‑148 (MQ=255)
                                                               aCAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGa                                                                                      >  2:2114134/1‑148 (MQ=255)
                                                                  ggACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGaaaa                                                                                   >  1:2589805/1‑148 (MQ=255)
                                                                   gACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAAc                                                                                  <  2:2555364/148‑1 (MQ=255)
                                                                     cTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCa                                                                                >  2:1857176/1‑148 (MQ=255)
                                                                     cTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCa                                                                                <  2:2243365/148‑1 (MQ=255)
                                                                      tATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTATCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCag                                                                               >  2:2199765/1‑148 (MQ=255)
                                                                        tCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCagaa                                                                             >  2:592226/1‑148 (MQ=255)
                                                                        tCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCagaa                                                                             >  2:813599/1‑148 (MQ=255)
                                                                         cTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTATCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCagaag                                                                            <  2:1015834/148‑1 (MQ=255)
                                                                          ttCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCagaaga                                                                           <  1:3344245/148‑1 (MQ=255)
                                                                          ttCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCagaaga                                                                           <  1:3374449/148‑1 (MQ=255)
                                                                          ttCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCagaaga                                                                           >  1:2332691/1‑148 (MQ=255)
                                                                           tcatcaTTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGAt                                                                          >  2:230512/1‑148 (MQ=255)
                                                                             atcatTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCg                                                                        <  1:2902615/148‑1 (MQ=255)
                                                                             atcatTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCg                                                                        >  1:3681437/1‑148 (MQ=255)
                                                                              tcatTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTATCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGa                                                                       <  2:531105/148‑1 (MQ=255)
                                                                               catTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGAt                                                                      >  2:1323370/1‑148 (MQ=255)
                                                                                 ttATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTa                                                                    >  2:173763/1‑148 (MQ=255)
                                                                                 ttATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTATCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTa                                                                    >  1:348378/1‑148 (MQ=255)
                                                                                   aTCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAAt                                                                  >  2:1868373/1‑148 (MQ=255)
                                                                                   aTCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAAt                                                                  <  2:1666030/148‑1 (MQ=255)
                                                                                    tCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAAtc                                                                 >  1:2338834/1‑148 (MQ=255)
                                                                                     cTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAAtct                                                                >  2:310587/1‑148 (MQ=255)
                                                                                     cTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAAtct                                                                >  2:2455363/1‑148 (MQ=255)
                                                                                     cTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAAtct                                                                <  2:1680298/148‑1 (MQ=255)
                                                                                        ttAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTATCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCgg                                                             <  1:2209423/148‑1 (MQ=255)
                                                                                            tGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCa                                                         <  2:2037936/148‑1 (MQ=255)
                                                                                            tGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCa                                                         >  1:436656/1‑148 (MQ=255)
                                                                                            tGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCa                                                         <  1:1736865/148‑1 (MQ=255)
                                                                                            tGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCa                                                         <  2:3333538/148‑1 (MQ=255)
                                                                                             gACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCaa                                                        <  1:3547353/148‑1 (MQ=255)
                                                                                             gACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCaa                                                        <  1:1981590/148‑1 (MQ=255)
                                                                                              aCTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAg                                                       <  1:2407889/148‑1 (MQ=255)
                                                                                               cTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAgg                                                      >  1:2564692/1‑148 (MQ=255)
                                                                                                ttACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAggg                                                     >  2:349264/1‑148 (MQ=255)
                                                                                                ttACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTATCTTCCAGACCTTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAggg                                                     >  1:1485599/1‑148 (MQ=255)
                                                                                                 tACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTATCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGt                                                    >  2:166291/1‑148 (MQ=255)
                                                                                                    gTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAg                                                 >  2:2895032/1‑148 (MQ=255)
                                                                                                        ttCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTATCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGa                                             <  2:2783796/148‑1 (MQ=255)
                                                                                                          cTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACg                                           <  1:2096524/148‑1 (MQ=255)
                                                                                                            gagaTACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGcc                                         <  1:853903/148‑1 (MQ=255)
                                                                                                               aTACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCAt                                      >  1:998859/1‑148 (MQ=255)
                                                                                                                tACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCAtt                                     >  1:77214/1‑148 (MQ=255)
                                                                                                                 aCTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTg                                    >  1:732143/1‑148 (MQ=255)
                                                                                                                      gCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCt                               >  1:2667008/1‑148 (MQ=255)
                                                                                                                       cGTCACAACGTGATTGGTGACAATACTTATCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTa                              <  1:2468142/148‑1 (MQ=255)
                                                                                                                         tCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAAt                            <  1:3732942/148‑1 (MQ=255)
                                                                                                                           acaACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAAttt                          <  2:3742511/148‑1 (MQ=255)
                                                                                                                              aCGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGt                       <  2:2400988/148‑1 (MQ=255)
                                                                                                                              aCGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGt                       >  2:1301679/1‑148 (MQ=255)
                                                                                                                                 tGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTa                    >  1:3769136/1‑148 (MQ=255)
                                                                                                                                 tGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTa                    >  1:1171029/1‑148 (MQ=255)
                                                                                                                                  gATTGGTGACAATACTTATCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTAt                   >  1:3272640/1‑148 (MQ=255)
                                                                                                                                    ttGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATgg                 >  2:1930801/1‑148 (MQ=255)
                                                                                                                                    ttGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATgg                 >  2:1400053/1‑148 (MQ=255)
                                                                                                                                         gACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATGGACGGa            <  1:1571391/148‑1 (MQ=255)
                                                                                                                                           cAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATGGACGGAGt          >  2:3642590/1‑148 (MQ=255)
                                                                                                                                             aTACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATGGACGGAGTaa        >  1:2149912/1‑148 (MQ=255)
                                                                                                                                              tACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATGGACGGAGTAAc       >  2:1600784/1‑148 (MQ=255)
                                                                                                                                                 ttCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATGGACGGAGTAACAAg    <  1:3180781/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TCGCTGTACCAGATAACCGGCACCATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATGGACGGAGTAACAAGCC  >  CP009273/4351426‑4351720

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