Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,351,573 C→A 27.1% 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,5730CA27.1% 571.9 / 163.4 258R60I (AGA→ATA) cadCcadBA operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (99/89);  new base A (43/27);  total (142/116)
Fisher's exact test for biased strand distribution p-value = 2.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 269 total aligned reads displayed.

TCGCTGTACCAGATAACCGGCACCATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATGGACGGAGTAACAAGCC  >  CP009273/4351426‑4351720
                                                                                                                                                   |                                                                                                                                                   
 cGCTGTACCAGATAACCGGCACCATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTtct                                                                                                                                                    <  1:2138628/148‑1 (MQ=255)
     gTACCAGATAACCGGCACCATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTc                                                                                                                                                <  2:1043766/148‑1 (MQ=255)
      tACCAGATAACCGGCACCATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTATCTTcc                                                                                                                                               >  1:2746399/1‑148 (MQ=255)
            aTAACCGGCACCATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACAt                                                                                                                                         <  1:1221718/148‑1 (MQ=255)
            aTAACCGGCACCATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACAt                                                                                                                                         >  1:3272938/1‑148 (MQ=255)
             tAACCGGCACCATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACAtt                                                                                                                                        >  2:622621/1‑148 (MQ=255)
              aaCCGGCACCATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTATCTTCCAGACATTa                                                                                                                                       >  2:2058506/1‑148 (MQ=255)
               aCCGGCACCATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTAt                                                                                                                                      >  1:1767468/1‑148 (MQ=255)
                ccGGCTCCATTAATTTATAGCCGCGCTTTTGTACAGTAGCGATATAGAAAGGTCTATCTTAATCATTATCTTTTAATGACTTACTTAGTTCTGATATACTCTGCGTCACAACGTGATTGGTTACAATACTTCTCTTCCAGACATTATc                                                                                                                                     <  2:1107715/148‑1 (MQ=255)
                 cGGCACCATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTAtag                                                                                                                                    <  2:1125844/148‑3 (MQ=255)
                 cGGCACCATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCg                                                                                                                                    <  2:1767468/148‑1 (MQ=255)
                 cGGCACCATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTATCTTCCAGACATTATCg                                                                                                                                    <  2:1926047/148‑1 (MQ=255)
                    cACCATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTATCTTCCAGACATTATCGATa                                                                                                                                 >  2:519592/1‑148 (MQ=255)
                       cATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGt                                                                                                                              <  2:784302/148‑1 (MQ=255)
                       cATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGt                                                                                                                              <  2:3047184/148‑1 (MQ=255)
                       cATTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGt                                                                                                                              >  2:254443/1‑148 (MQ=255)
                        aTTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGtt                                                                                                                             >  2:1350746/1‑148 (MQ=255)
                        aTTAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGtt                                                                                                                             >  1:3541507/1‑148 (MQ=255)
                         ttAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTc                                                                                                                            <  2:881314/148‑1 (MQ=255)
                         ttAATTTATAGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTATCTTCCAGACATTATCGATAAGTTc                                                                                                                            <  2:3331071/148‑1 (MQ=255)
                               tataGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCt                                                                                                                      <  2:224617/148‑1 (MQ=255)
                               tataGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCt                                                                                                                      >  1:1568484/1‑148 (MQ=255)
                               tataGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTATCTTCCAGACATTATCGATAAGTTCATCCCt                                                                                                                      >  2:834213/1‑148 (MQ=255)
                                ataGCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTg                                                                                                                     >  1:117639/1‑148 (MQ=255)
                                   gCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTa                                                                                                                  >  2:3159110/1‑148 (MQ=255)
                                   gCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTa                                                                                                                  <  2:1202599/148‑1 (MQ=255)
                                   gCCGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTa                                                                                                                  <  2:2326550/148‑1 (MQ=255)
                                    ccGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTaa                                                                                                                 >  2:320165/1‑148 (MQ=255)
                                    ccGCGCTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTaa                                                                                                                 <  1:3463981/148‑1 (MQ=255)
                                      gcgcTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGt                                                                                                               <  1:1475528/148‑1 (MQ=255)
                                      gcgcTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGt                                                                                                               >  1:869089/1‑148 (MQ=255)
                                      gcgcTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTATCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGt                                                                                                               <  1:238303/148‑1 (MQ=255)
                                        gcTTTGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTAc                                                                                                             <  2:651082/148‑1 (MQ=255)
                                          tttGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACtt                                                                                                           >  2:1446604/1‑148 (MQ=255)
                                          tttGGTACAGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTATCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACtt                                                                                                           >  1:3373749/1‑148 (MQ=255)
                                                  aGTAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTCAGTACTTCGCCACtg                                                                                                   >  2:2417575/1‑148 (MQ=255)
                                                    tAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACtgtg                                                                                                 >  1:2493717/1‑148 (MQ=255)
                                                    tAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTATCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACtgtg                                                                                                 <  2:891801/148‑1 (MQ=255)
                                                    tAGCGATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTATCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACtgtg                                                                                                 >  1:2823002/1‑148 (MQ=255)
                                                        gATATAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGa                                                                                             <  2:2536812/148‑1 (MQ=255)
                                                         atatAGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAg                                                                                            >  1:767456/1‑148 (MQ=255)
                                                          tataGACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGc                                                                                           >  2:2139793/1‑148 (MQ=255)
                                                              gACAGGACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAg                                                                                       >  1:303218/1‑148 (MQ=255)
                                                                  ggACTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGaaaa                                                                                   >  1:3362900/1‑148 (MQ=255)
                                                                    aCTATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTATCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAAcc                                                                                 >  2:1193652/1‑148 (MQ=255)
                                                                      tATCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTATCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCag                                                                               >  2:263057/1‑148 (MQ=255)
                                                                       aTCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCaga                                                                              >  2:856952/1‑148 (MQ=255)
                                                                       aTCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTATCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCaga                                                                              >  1:990028/1‑148 (MQ=255)
                                                                       aTCTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTATCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCaga                                                                              >  2:3072310/1‑148 (MQ=255)
                                                                         cTTCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTATCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCagaag                                                                            <  1:70038/148‑1 (MQ=255)
                                                                          ttCATCATTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCagaaga                                                                           <  2:3305757/148‑1 (MQ=255)
                                                                             atcatTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCg                                                                        <  2:3551927/148‑1 (MQ=255)
                                                                               catTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTATCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGAt                                                                      <  2:2970223/148‑1 (MQ=255)
                                                                                atTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGAtt                                                                     <  2:3439610/148‑1 (MQ=255)
                                                                                atTATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATACCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGAtt                                                                     >  1:521527/1‑148 (MQ=255)
                                                                                  tATCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTaa                                                                   >  1:2223396/1‑148 (MQ=255)
                                                                                   aTCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAAt                                                                  >  1:1368334/1‑148 (MQ=255)
                                                                                   aTCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAAt                                                                  <  1:416284/148‑1 (MQ=255)
                                                                                   aTCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTATCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTTCTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAAt                                                                  <  2:2514546/148‑1 (MQ=255)
                                                                                   aTCTTTTAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTATCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAAt                                                                  >  2:342968/1‑148 (MQ=255)
                                                                                      ttttAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAAtctc                                                               >  1:1976088/1‑148 (MQ=255)
                                                                                        ttAATGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCgg                                                             >  2:1809242/1‑148 (MQ=255)
                                                                                          aaTGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTATCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCt                                                           <  2:187443/148‑1 (MQ=255)
                                                                                           aTGACTTACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTATCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTc                                                          >  2:1099368/1‑148 (MQ=255)
                                                                                                 tACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGt                                                    >  2:1410315/1‑148 (MQ=255)
                                                                                                 tACGTAGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTATCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGt                                                    >  1:1056564/1‑148 (MQ=255)
                                                                                                      aGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGtt                                               >  2:42571/1‑148 (MQ=255)
                                                                                                      aGTTCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTATCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGtt                                               <  1:921652/148‑1 (MQ=255)
                                                                                                        ttCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGa                                             <  2:2883771/148‑1 (MQ=255)
                                                                                                        ttCTGAGATACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGa                                             <  2:2190307/148‑1 (MQ=255)
                                                                                                             agaTACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGccc                                        >  1:124425/1‑148 (MQ=255)
                                                                                                              gaTACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCa                                       <  2:3240612/148‑1 (MQ=255)
                                                                                                                tACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCAtt                                     >  2:2059684/1‑148 (MQ=255)
                                                                                                                tACTCTGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCAtt                                     <  1:546101/148‑1 (MQ=255)
                                                                                                                 aCTCTGCGTCACAACGTGATTGGTGACAATACTTATCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTg                                    >  1:2966084/1‑148 (MQ=255)
                                                                                                                  ctctGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGc                                   >  1:1307784/1‑148 (MQ=255)
                                                                                                                    ctGCGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCgg                                 <  1:2873909/148‑1 (MQ=255)
                                                                                                                     tGCGTCACAACGTGATTGGTGACAATACTTATCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGc                                <  1:3544722/148‑1 (MQ=255)
                                                                                                                       cGTCACAACGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTa                              <  1:1119606/148‑1 (MQ=255)
                                                                                                                               cGTGATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGtt                      >  2:1544394/1‑148 (MQ=255)
                                                                                                                               cGTGATTGGTGACAATACTTATCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGtt                      >  2:748795/1‑148 (MQ=255)
                                                                                                                                  gATTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTAt                   >  1:2543552/1‑148 (MQ=255)
                                                                                                                                   aTTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATg                  <  1:1409843/148‑1 (MQ=255)
                                                                                                                                   aTTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATg                  <  1:2184317/148‑1 (MQ=255)
                                                                                                                                   aTTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATg                  <  1:1411400/148‑1 (MQ=255)
                                                                                                                                   aTTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATg                  <  1:1545103/148‑1 (MQ=255)
                                                                                                                                   aTTGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATg                  >  2:1801273/1‑148 (MQ=255)
                                                                                                                                    ttGGTGACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATgg                 <  1:2736791/148‑1 (MQ=255)
                                                                                                                                         gACAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATGGACGGa            >  2:1419002/1‑148 (MQ=255)
                                                                                                                                           cAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATGGACGGAGt          >  1:1732435/1‑148 (MQ=255)
                                                                                                                                           cAATACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATGGACGGAGt          >  1:2065551/1‑148 (MQ=255)
                                                                                                                                            aaTACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATGGACGGAGTa         <  2:3065312/148‑1 (MQ=255)
                                                                                                                                            aaTACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATGGACGGAGTa         <  2:1986504/148‑1 (MQ=255)
                                                                                                                                            aaTACTTATCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATGGACGGAGTa         <  1:58307/148‑1 (MQ=255)
                                                                                                                                             aTACTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATGGACGGAGTaa        >  2:2124647/1‑148 (MQ=255)
                                                                                                                                             aTACTTATCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATGGACGGAGTaa        >  1:190493/1‑148 (MQ=255)
                                                                                                                                               aCTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATGGACGGAGTAACa      <  1:1016687/148‑1 (MQ=255)
                                                                                                                                                cTTCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATGGACGGAGTAACaa     >  1:1081629/1‑148 (MQ=255)
                                                                                                                                                cTTATCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATGGACGGAGTAACaa     <  2:2746399/148‑1 (MQ=255)
                                                                                                                                                 ttCTCTTCCAGACATTATCGATAAGTTCATCCCTGCTAAGTACTTCGCCACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATGGACGGAGTAACAAg    <  1:632411/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: