Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,919,922 G→T 100% D86Y (GAT→TAT)  ygdH → UPF0717 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,919,9220GT98.2% 892.4 / 5.4 281D86Y (GAT→TAT) ygdHUPF0717 family protein
Reads supporting (aligned to +/- strand):  ref base G (4/1);  new base T (130/146);  total (134/147)
Fisher's exact test for biased strand distribution p-value = 1.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 281 total aligned reads displayed.

TTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTT  >  CP009273/2919777‑2920066
                                                                                                                                                 |                                                                                                                                                
ttCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTAt                                                                                                                                                <  1:720204/148‑1 (MQ=255)
ttCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTAt                                                                                                                                                >  1:1468708/1‑148 (MQ=255)
 tCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATg                                                                                                                                               <  2:777670/148‑1 (MQ=255)
      tGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGa                                                                                                                                          >  1:1755177/1‑148 (MQ=255)
       gCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGaa                                                                                                                                         <  2:1421973/148‑1 (MQ=255)
        ccGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAat                                                                                                                                        >  1:2669652/1‑148 (MQ=255)
         cGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAatt                                                                                                                                       <  1:2525437/148‑1 (MQ=255)
         cGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAatt                                                                                                                                       <  1:2885034/148‑1 (MQ=255)
          gTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAatta                                                                                                                                      <  2:2915342/148‑1 (MQ=255)
           tACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAattat                                                                                                                                     <  2:1192781/148‑1 (MQ=255)
              tGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTcg                                                                                                                                  <  2:2684491/148‑1 (MQ=255)
                 aCTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCgcgc                                                                                                                               <  2:1911281/148‑1 (MQ=255)
                 aCTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCgcgc                                                                                                                               <  1:2018954/148‑1 (MQ=255)
                  cgccGGTAGTTTGACCGATAACAGCAAATAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCGTGCGCCGTGAACGCGGCGTAAATCTGGAACTGATTTATCCCCCGGAAGATGCTTTTGTCTATGGGCGAATTATTCGCGCt                                                                                                                              <  1:1994349/146‑1 (MQ=255)
                  cTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCt                                                                                                                              <  1:3557099/148‑1 (MQ=255)
                  cTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCt                                                                                                                              <  1:696191/148‑1 (MQ=255)
                   tCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCtt                                                                                                                             >  2:3102299/1‑148 (MQ=255)
                       gTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCa                                                                                                                         <  1:1874512/148‑1 (MQ=255)
                          gTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGc                                                                                                                      <  1:680837/148‑1 (MQ=255)
                            ttGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCa                                                                                                                    >  1:2581338/1‑148 (MQ=255)
                                ccGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACtt                                                                                                                <  1:1019337/148‑1 (MQ=255)
                                 cGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTg                                                                                                               >  1:3274175/1‑148 (MQ=255)
                                 cGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTg                                                                                                               >  1:1668524/1‑148 (MQ=255)
                                  gATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGt                                                                                                              <  1:1390720/148‑1 (MQ=255)
                                   aTAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGGATTATTCGCGCTTTGCAGGCCAACTTGtt                                                                                                             >  1:1761996/1‑148 (MQ=255)
                                   aTAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGtt                                                                                                             <  1:1808633/148‑1 (MQ=255)
                                   aTAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGtt                                                                                                             <  1:1623336/148‑1 (MQ=255)
                                    tAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTc                                                                                                            >  2:590410/1‑148 (MQ=255)
                                       cAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTcgcg                                                                                                         >  2:3027708/1‑148 (MQ=255)
                                         gCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGt                                                                                                       >  1:1149564/1‑148 (MQ=255)
                                          cAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTc                                                                                                      <  2:1938953/148‑1 (MQ=255)
                                          cAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTc                                                                                                      <  1:1245803/148‑1 (MQ=255)
                                           aaaGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTcc                                                                                                     >  2:3440380/1‑148 (MQ=255)
                                           aaaGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTcc                                                                                                     >  2:1679809/1‑148 (MQ=255)
                                            aaGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCt                                                                                                    >  2:2821870/1‑148 (MQ=255)
                                            aaGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCt                                                                                                    >  1:2043461/1‑148 (MQ=255)
                                             aGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTg                                                                                                   >  1:2103611/1‑148 (MQ=255)
                                             aGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTg                                                                                                   >  1:1457668/1‑148 (MQ=255)
                                              gAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGc                                                                                                  >  1:3329136/1‑148 (MQ=255)
                                               aaTTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCg                                                                                                 >  2:2599967/1‑148 (MQ=255)
                                               aaTTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCg                                                                                                 >  1:2630399/1‑148 (MQ=255)
                                                aTTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGt                                                                                                <  1:1713580/148‑1 (MQ=255)
                                                 ttGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTg                                                                                               >  2:1827987/1‑148 (MQ=255)
                                                      gTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACAtt                                                                                          >  1:3412697/1‑148 (MQ=255)
                                                         tcGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATtctc                                                                                       >  2:2013821/1‑148 (MQ=255)
                                                          cGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTctct                                                                                      <  1:487201/148‑1 (MQ=255)
                                                          cGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTctct                                                                                      <  2:2760880/148‑1 (MQ=255)
                                                          cGTTTTGAAAATTTCGATATGAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTctct                                                                                      <  2:797809/148‑1 (MQ=255)
                                                             tttGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCg                                                                                   >  1:3373764/1‑148 (MQ=255)
                                                             tttGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCg                                                                                   >  1:3373659/1‑148 (MQ=255)
                                                              ttGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGt                                                                                  >  1:2314120/1‑148 (MQ=255)
                                                              ttGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGt                                                                                  <  2:3354166/148‑1 (MQ=255)
                                                              ttGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGt                                                                                  <  1:1233117/148‑1 (MQ=255)
                                                               tGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGtt                                                                                 <  1:2095847/148‑1 (MQ=255)
                                                               tGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGtt                                                                                 <  2:2086386/148‑1 (MQ=255)
                                                                    aTTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACgg                                                                            >  1:1206677/1‑148 (MQ=255)
                                                                     tttCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACggg                                                                           >  2:2188723/1‑148 (MQ=255)
                                                                      ttCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGc                                                                          >  2:501516/1‑148 (MQ=255)
                                                                       tCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCa                                                                         >  1:1307873/1‑148 (MQ=255)
                                                                         gATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCaaa                                                                       <  1:1525632/148‑1 (MQ=255)
                                                                          atatTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAAt                                                                      <  2:1018527/148‑1 (MQ=255)
                                                                           tatTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATc                                                                     >  1:919722/1‑148 (MQ=255)
                                                                           tatTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATc                                                                     >  1:2932481/1‑148 (MQ=255)
                                                                            atTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATcc                                                                    <  1:422891/148‑1 (MQ=255)
                                                                             ttAACGTCGTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCa                                                                   <  2:2046659/148‑1 (MQ=255)
                                                                               aaCGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATa                                                                 <  1:2987647/148‑1 (MQ=255)
                                                                                 cGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAac                                                               >  1:1918684/1‑148 (MQ=255)
                                                                                   tCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAacac                                                             >  2:1348340/1‑148 (MQ=255)
                                                                                    cTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACAcc                                                            >  1:1112356/1‑148 (MQ=255)
                                                                                    cTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACAcc                                                            >  2:1470363/1‑148 (MQ=255)
                                                                                    cTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACAcc                                                            >  1:3546406/1‑148 (MQ=255)
                                                                                     ttGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACcg                                                           <  2:71624/148‑1 (MQ=255)
                                                                                      tGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACcgt                                                          >  1:1360725/1‑148 (MQ=255)
                                                                                      tGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACcgt                                                          <  2:390133/148‑1 (MQ=255)
                                                                                             tGAACGCGGCGTAAAGCTGGAAGTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGtttt                                                   >  1:2394549/1‑148 (MQ=255)
                                                                                                 cgcgGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCa                                               <  1:1687750/148‑1 (MQ=255)
                                                                                                   cggcgTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAAc                                             >  2:524524/1‑148 (MQ=255)
                                                                                                        tAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGaa                                        <  1:1510087/148‑1 (MQ=255)
                                                                                                          aaGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAAtc                                      <  2:2049041/148‑1 (MQ=255)
                                                                                                          aaGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAAtc                                      <  1:3026151/148‑1 (MQ=255)
                                                                                                                gAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACa                                >  1:1335928/1‑148 (MQ=255)
                                                                                                                gAACTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACa                                >  2:1288517/1‑148 (MQ=255)
                                                                                                                 aaCTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACaa                               <  1:3111890/148‑1 (MQ=255)
                                                                                                                   cTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGACTCTCGACAACt                             >  2:1962809/1‑148 (MQ=255)
                                                                                                                   cTGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACt                             >  2:1612148/1‑148 (MQ=255)
                                                                                                                    tGATTAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTc                            <  1:3042075/148‑1 (MQ=255)
                                                                                                                        tAATCCCCCGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTc                        <  2:3329136/148‑1 (MQ=255)
                                                                                                                          aTCCCCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTATTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCca                      >  2:873379/1‑148 (MQ=255)
                                                                                                                           gcccccGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCcac                     <  2:3198850/147‑1 (MQ=255)
                                                                                                                            cccccGGAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCcaca                    >  2:174384/1‑148 (MQ=255)
                                                                                                                            accccGGAAGAGGCTTTTGTCTATGGGCGAATTAATCGCGCTTTTCAGGCAAACTTGTTCGCGGTCCTGAGTGACATTCTCTTCGTTTACGGGCAATTACATAACACCGTTCGTTTTCCCAACCTGAATCTCTACAACTCCGTCcaca                    >  2:2466273/2‑148 (MQ=255)
                                                                                                                                   aaGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCTGTCCTGCGTGACATTCTCTTCGTTTAAGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTaa             <  2:2394560/148‑1 (MQ=255)
                                                                                                                                   aaGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTaa             >  2:1131166/1‑148 (MQ=255)
                                                                                                                                    agagGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACAtaactaac            <  2:2630399/148‑1 (MQ=255)
                                                                                                                                    agagGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACAtaactaac            <  1:2280340/148‑1 (MQ=255)
                                                                                                                                      agGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCt          <  1:501553/148‑1 (MQ=255)
                                                                                                                                        gCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTgg        >  1:2829161/1‑148 (MQ=255)
                                                                                                                                          ttttGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTc      <  2:1532624/148‑1 (MQ=255)
                                                                                                                                           tttGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCt     <  2:1209004/148‑1 (MQ=255)
                                                                                                                                           tttGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCt     <  2:1473234/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTT  >  CP009273/2919777‑2920066

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: