Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,077,698 A→T 10.7% intergenic (+34/+9) sbcB → / ← yeeD exodeoxyribonuclease I; exonuclease I/putative TusA family sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,077,6980AT10.7% 743.4 / 66.4 259intergenic (+34/+9)sbcB/yeeDexodeoxyribonuclease I; exonuclease I/putative TusA family sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base A (134/97);  new base T (11/17);  total (145/114)
Fisher's exact test for biased strand distribution p-value = 7.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

Only 100 of 266 total aligned reads displayed.

TTTTTGCAGGGTTATGCTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCAATCAAACTCAATCACTAACTCATC  >  CP009273/2077551‑2077844
                                                                                                                                                   |                                                                                                                                                  
tttttGCAGGGTTATGCTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACa                                                                                                                                                    <  1:1263595/148‑1 (MQ=255)
tttttGCAGGGTTATGCTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACa                                                                                                                                                    >  2:3298500/1‑148 (MQ=255)
tttttGCAGGGTTATGCTGATGAATTGCAGATGCTGGTACAAAAATATGCCGATGACAAAGAGAAAGTGGCGCTCTTAAAAGCACTTTGGCAGTACGCGGAAGAGTTTGACTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACa                                                                                                                                                    <  1:2250890/148‑1 (MQ=255)
           ttATGCTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTct                                                                                                                                         <  2:1603506/148‑2 (MQ=255)
            tATGCTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGACAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAg                                                                                                                                        <  1:3541930/148‑1 (MQ=255)
              tGCTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGc                                                                                                                                      >  2:2348488/1‑148 (MQ=255)
               gCTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACACAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCt                                                                                                                                     <  2:2753625/148‑1 (MQ=255)
                cTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCtt                                                                                                                                    <  2:1334603/148‑1 (MQ=255)
                 tgatgaATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACTATGTTGAGTCAGGCttt                                                                                                                                   <  1:203663/148‑1 (MQ=255)
                  gatgaATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCtttt                                                                                                                                  >  1:1786170/1‑148 (MQ=255)
                   atgaATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACTATGTTGAGTCAGGCttttt                                                                                                                                 >  1:980904/1‑148 (MQ=255)
                   atgaATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCttttt                                                                                                                                 >  2:486356/1‑148 (MQ=255)
                   atgaATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCttttt                                                                                                                                 <  2:3222402/148‑1 (MQ=255)
                     gaATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGa                                                                                                                               >  2:443040/1‑148 (MQ=255)
                       aTTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAAc                                                                                                                             >  1:1295292/1‑148 (MQ=255)
                       aTTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAAc                                                                                                                             >  1:2691259/1‑148 (MQ=255)
                        ttGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACg                                                                                                                            >  2:1453295/1‑148 (MQ=255)
                           cAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTg                                                                                                                         <  1:1601215/148‑1 (MQ=255)
                             gATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGAt                                                                                                                       >  2:3505096/1‑148 (MQ=255)
                              aTGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATg                                                                                                                      >  2:2659215/1‑148 (MQ=255)
                               tGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTGTTGAACGGTGATGc                                                                                                                     >  1:446502/1‑148 (MQ=255)
                                 cTGGTACAACAATCTGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTc                                                                                                                   <  2:1615453/148‑1 (MQ=255)
                                     tACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACg                                                                                                               <  1:1025484/148‑1 (MQ=255)
                                     tACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACg                                                                                                               <  2:3643214/148‑1 (MQ=255)
                                          caaTATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCa                                                                                                          >  2:740968/1‑148 (MQ=255)
                                           aaTATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCAt                                                                                                         <  1:2406506/148‑1 (MQ=255)
                                            atatGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATc                                                                                                        >  2:399441/1‑148 (MQ=255)
                                                 ccGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCaa                                                                                                   >  2:165080/1‑148 (MQ=255)
                                                  cGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACTATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAAt                                                                                                  <  1:2457972/148‑1 (MQ=255)
                                                  cGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAAt                                                                                                  <  1:2550368/148‑1 (MQ=255)
                                                  cGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAAt                                                                                                  <  1:115132/148‑1 (MQ=255)
                                                     tGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATctg                                                                                               <  2:2371820/148‑1 (MQ=255)
                                                       aCAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATctgct                                                                                             <  1:2711059/148‑1 (MQ=255)
                                                         aaaGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACTATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGa                                                                                           <  2:216290/148‑1 (MQ=255)
                                                          aaGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACTATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGAt                                                                                          >  2:2983856/1‑148 (MQ=255)
                                                          aaGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGAt                                                                                          >  1:2102070/1‑148 (MQ=255)
                                                          aaGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGAt                                                                                          >  2:1673109/1‑148 (MQ=255)
                                                           agagAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACTATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATa                                                                                         <  1:869915/148‑1 (MQ=255)
                                                           agagAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACTATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATa                                                                                         <  1:2780843/148‑1 (MQ=255)
                                                            gagaAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCt                                                                                                                    >  1:891793/1‑120 (MQ=255)
                                                            gagaAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATaa                                                                                        >  1:2605975/1‑148 (MQ=255)
                                                            gagaAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATaa                                                                                        >  2:2044375/1‑148 (MQ=255)
                                                             agaAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAAt                                                                                       >  2:536341/1‑148 (MQ=255)
                                                              gaAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGatcatc                                                                                      >  2:1121387/1‑148 (MQ=255)
                                                              gaAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATc                                                                                      <  1:3298500/148‑1 (MQ=255)
                                                               aaaGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCg                                                                                     >  1:1522676/1‑148 (MQ=255)
                                                                   tGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGa                                                                                 >  2:1972188/1‑148 (MQ=255)
                                                                   tGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGa                                                                                 <  2:1339270/148‑1 (MQ=255)
                                                                        cTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCa                                                                            <  1:3444052/148‑1 (MQ=255)
                                                                          gTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATg                                                                          >  2:35117/1‑148 (MQ=255)
                                                                          gTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATg                                                                          <  1:601306/148‑1 (MQ=255)
                                                                            tAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTc                                                                        >  1:2255939/1‑148 (MQ=255)
                                                                              aaaGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCct                                                                      >  1:1420809/1‑148 (MQ=255)
                                                                              aaaGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCct                                                                      >  1:2624257/1‑148 (MQ=255)
                                                                              aaaGCACGTTGGCAGTACGCGGAAGAGATTGTCTAAATTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCct                                                                      <  2:2443379/148‑1 (MQ=255)
                                                                               aaGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCctt                                                                     <  2:1850523/148‑1 (MQ=255)
                                                                                aGCACTTTGGCAGTACGCGGAAGAGATTGTCTATTTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCcttc                                                                    <  2:2222178/148‑1 (MQ=255)
                                                                                aGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCcttc                                                                    >  2:1349083/1‑148 (MQ=255)
                                                                                aGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCcttc                                                                    >  2:2546199/1‑148 (MQ=255)
                                                                                aGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCcttc                                                                    >  2:3251655/1‑148 (MQ=255)
                                                                                aGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCcttc                                                                    >  1:854241/1‑148 (MQ=255)
                                                                                 gCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACTATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCcttct                                                                   <  1:2880164/148‑1 (MQ=255)
                                                                                 gCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCcttct                                                                   >  2:3393523/1‑148 (MQ=255)
                                                                                  cACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACTATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCcttctt                                                                  >  2:3759381/1‑148 (MQ=255)
                                                                                  cACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCcttctt                                                                  >  2:2600302/1‑148 (MQ=255)
                                                                                  cACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCcttctt                                                                  >  1:362803/1‑148 (MQ=255)
                                                                                   aCTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACTATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCttcttc                                                                 >  1:371916/1‑148 (MQ=255)
                                                                                    cTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTtcttct                                                                >  2:861071/1‑148 (MQ=255)
                                                                                             tACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCAc                                                       >  1:2965819/1‑148 (MQ=255)
                                                                                             tACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCAc                                                       >  1:516239/1‑148 (MQ=255)
                                                                                                    aaGAGATTGTCTACTTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGTGa                                                <  1:3137095/148‑1 (MQ=255)
                                                                                                    aaGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGa                                                <  1:1408745/148‑1 (MQ=255)
                                                                                                    aaGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGa                                                <  1:1939650/148‑1 (MQ=255)
                                                                                                     agagATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGAt                                               <  1:2865022/148‑1 (MQ=255)
                                                                                                        gaTTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGc                                            <  2:2836061/148‑1 (MQ=255)
                                                                                                          ttGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCtt                                          <  2:2857064/148‑1 (MQ=255)
                                                                                                             tCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAg                                       >  1:1063157/1‑148 (MQ=255)
                                                                                                             tCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAg                                       >  1:269738/1‑148 (MQ=255)
                                                                                                             tCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAg                                       >  1:3301515/1‑148 (MQ=255)
                                                                                                             tCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAg                                       <  1:3516621/148‑1 (MQ=255)
                                                                                                             tCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAg                                       >  1:3592598/1‑148 (MQ=255)
                                                                                                              cTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGt                                      >  2:2541378/1‑148 (MQ=255)
                                                                                                              cTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGt                                      >  1:3079376/1‑148 (MQ=255)
                                                                                                               tAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGtt                                     >  2:1921354/1‑148 (MQ=255)
                                                                                                                  tttGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGTCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGcc                                  <  1:1157182/148‑1 (MQ=255)
                                                                                                                   ttGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCt                                 >  1:2968383/1‑148 (MQ=255)
                                                                                                                       aGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGt                             >  2:3586051/1‑148 (MQ=255)
                                                                                                                         ccTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGc                           >  1:2272584/1‑148 (MQ=255)
                                                                                                                          cTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCa                          <  2:1591236/148‑1 (MQ=255)
                                                                                                                                 gcTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCAATCAAAc                   >  2:2630954/1‑148 (MQ=255)
                                                                                                                                 gcTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCAATCAAAc                   <  1:1349083/148‑1 (MQ=255)
                                                                                                                                      tACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCAATCAAACTCAAt              <  2:2102070/148‑1 (MQ=255)
                                                                                                                                        cggcggAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCAATCAAACTCAATCa            <  1:2357608/148‑1 (MQ=255)
                                                                                                                                        cggcggAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCAATCAAACTCAATCa            <  1:1579772/148‑1 (MQ=255)
                                                                                                                                            ggatgaCAATTTTGAGTCAGGCTTTTTGAACGGTGATGCTACACGCTGCATCGCCAATCTTCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCAATCAAACTATATCACTaa        <  2:920627/148‑1 (MQ=255)
                                                                                                                                             gaggaCAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCAATCAAACTCAATCACTAAc       >  1:1458343/1‑148 (MQ=255)
                                                                                                                                             gaggaCAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCAATCAAACTCAATCACTAAc       >  1:2984958/1‑148 (MQ=255)
                                                                                                                                               ggaCAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCAATCAAACTCAATCACTAACTc     >  2:78945/1‑148 (MQ=255)
                                                                                                                                                gaCTATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCAATCAAACTCAATCACTAACTCa    <  2:3525079/148‑1 (MQ=255)
                                                                                                                                                 aCAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCAATCAAACTCAATCACTAACTCAt   >  1:3190542/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTTTTGCAGGGTTATGCTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCAATCAAACTCAATCACTAACTCATC  >  CP009273/2077551‑2077844

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: