Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,077,698 A→T 20.3% 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,6980AT20.3% 415.6 / 103.4 187intergenic (+34/+9)sbcB/yeeDexodeoxyribonuclease I; exonuclease I/putative TusA family sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base A (87/62);  new base T (15/23);  total (102/85)
Fisher's exact test for biased strand distribution p-value = 4.49e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

Only 100 of 191 total aligned reads displayed.

TTTTTGCAGGGTTATGCTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCAATCAAACTCAATCACTAACTCATCG  >  CP009273/2077551‑2077845
                                                                                                                                                   |                                                                                                                                                   
tttttGCAGGGTTATGCTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACa                                                                                                                                                     <  2:1415350/148‑1 (MQ=255)
tttttGCAGGGTTATGCTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTATAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACa                                                                                                                                                     <  2:2100/148‑1 (MQ=255)
  tttGCAGGGTTATGCTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAACAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACTAt                                                                                                                                                   <  2:94676/148‑1 (MQ=255)
   ttGCAGGGTTATGCTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATg                                                                                                                                                  <  2:511698/148‑1 (MQ=255)
   ttGCAGGGTTATGCTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATg                                                                                                                                                  <  1:2381988/148‑1 (MQ=255)
    tGCAGGGTTATGCTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGt                                                                                                                                                 <  1:991639/148‑1 (MQ=255)
    tGCAGGGTTATGCTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGt                                                                                                                                                 <  1:1974327/148‑1 (MQ=255)
             aTGCTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACTATGTTGAGTCAgg                                                                                                                                        <  1:1232081/148‑1 (MQ=255)
              tGCTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACTATGTTGAGTCAGGc                                                                                                                                       <  2:1231517/148‑1 (MQ=255)
              tGCTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGc                                                                                                                                       >  1:2272058/1‑148 (MQ=255)
                 tgatgaATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCttt                                                                                                                                    <  2:957778/148‑1 (MQ=255)
                  gatgaATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCtttt                                                                                                                                   <  1:105632/148‑1 (MQ=255)
                   atgaATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCttttt                                                                                                                                  >  1:3003957/1‑148 (MQ=255)
                   atgaATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCttttt                                                                                                                                  >  1:3521761/1‑148 (MQ=255)
                    ggaATTGCAGATGCTGGTACAACACTATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACTATGTTGAGTCAGGCTTTTTg                                                                                                                                 <  2:19846/147‑1 (MQ=255)
                       aTTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAAc                                                                                                                              >  2:3532080/1‑148 (MQ=255)
                         tGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACgg                                                                                                                            <  1:2391080/148‑1 (MQ=255)
                           cAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTg                                                                                                                          <  2:1888266/148‑1 (MQ=255)
                            aGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGa                                                                                                                         <  1:2987930/148‑1 (MQ=255)
                                gCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACTATGTTGAGTCAGGCTTTTTGAACGGTGATGCt                                                                                                                     >  1:2777149/1‑148 (MQ=255)
                                gCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCt                                                                                                                     >  2:657411/1‑148 (MQ=255)
                                gCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCt                                                                                                                     >  2:2449350/1‑148 (MQ=255)
                                 cTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACTATGTTGAGTCAGGCTTTTTGAACGGTGATGCTc                                                                                                                    <  2:1941987/148‑1 (MQ=255)
                                   ggTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGAACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGTGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCa                                                                                                                  >  1:3372139/1‑148 (MQ=255)
                                        aacaaTATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTg                                                                                                             >  1:842955/1‑148 (MQ=255)
                                         acaaTATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGc                                                                                                            >  1:824768/1‑148 (MQ=255)
                                             tatGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCg                                                                                                        >  2:415311/1‑148 (MQ=255)
                                              atGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGc                                                                                                       >  1:1492229/1‑148 (MQ=255)
                                                  cGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAAt                                                                                                   <  1:1252830/148‑1 (MQ=255)
                                                     tGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATctg                                                                                                <  2:410716/148‑1 (MQ=255)
                                                     tGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATctg                                                                                                >  2:749055/1‑148 (MQ=255)
                                                       aCAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACTATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATctgct                                                                                              >  2:2020288/1‑148 (MQ=255)
                                                       aCAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATctgct                                                                                              <  1:166494/148‑1 (MQ=255)
                                                        cAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTAAGGCGGATGACAATGTTGAGTCAGGCTTTTTGAAAGGTGATGCTCCACGCTGCATCGCCAATatgctt                                                                                             >  1:438285/1‑147 (MQ=255)
                                                         aaaGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGa                                                                                            >  2:2801752/1‑148 (MQ=255)
                                                           agagAAAGTGGCGCTGATAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACTATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATa                                                                                          <  2:2489711/148‑1 (MQ=255)
                                                            gagaAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACTATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATaa                                                                                         >  2:3058874/1‑148 (MQ=255)
                                                            gagaAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCGCGCTGCATCGCCAATCTGCTGATaa                                                                                         <  2:1921194/148‑1 (MQ=255)
                                                            gagaAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATaa                                                                                         >  1:967857/1‑148 (MQ=255)
                                                              gaAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATc                                                                                       >  2:2600894/1‑148 (MQ=255)
                                                                aaGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCgg                                                                                     <  2:3003957/148‑1 (MQ=255)
                                                                  gTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACTATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTg                                                                                   <  2:2781369/148‑1 (MQ=255)
                                                                       gcTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGc                                                                              <  1:162192/148‑1 (MQ=255)
                                                                        cTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCa                                                                             >  1:1427259/1‑148 (MQ=255)
                                                                        cTGATAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACTATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCa                                                                             <  1:1014302/148‑1 (MQ=255)
                                                                          gTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGACTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATg                                                                           >  2:1551223/1‑148 (MQ=255)
                                                                           ttAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACTATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGt                                                                          >  2:1674854/1‑148 (MQ=255)
                                                                           ttAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGt                                                                          >  1:1879449/1‑148 (MQ=255)
                                                                             aaaaTCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCAACGCTGCATCGCCAATCTGATGATCATCGGTGATGGCATGTcc                                                                        >  2:1844944/1‑148 (MQ=255)
                                                                              aaaGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCct                                                                       >  2:2163750/1‑148 (MQ=255)
                                                                              aaaGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCct                                                                       >  1:3364853/1‑148 (MQ=255)
                                                                               aaGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACTATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCctt                                                                      <  1:1458116/148‑1 (MQ=255)
                                                                               aaGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCctt                                                                      >  1:3083778/1‑148 (MQ=255)
                                                                               aaGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCctt                                                                      <  1:1594439/148‑1 (MQ=255)
                                                                                aGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCcttc                                                                     >  2:2266033/1‑148 (MQ=255)
                                                                                    cTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTtcttct                                                                 >  1:1375381/1‑148 (MQ=255)
                                                                                            gTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCa                                                         <  2:2492510/148‑1 (MQ=255)
                                                                                             tACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACTATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCAc                                                        >  2:2362676/1‑148 (MQ=255)
                                                                                              aCGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACt                                                       >  1:2749323/1‑148 (MQ=255)
                                                                                                   gtagagATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTggg                                                  <  2:1534420/146‑1 (MQ=255)
                                                                                                   gAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTggg                                                  >  1:2770263/1‑148 (MQ=255)
                                                                                                   gAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTggg                                                  >  1:763573/1‑148 (MQ=255)
                                                                                                   gAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTggg                                                  <  2:1458784/148‑1 (MQ=255)
                                                                                                    aaGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACTATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGa                                                 <  1:1674854/148‑1 (MQ=255)
                                                                                                     agagATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACTATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGAt                                                >  1:576266/1‑148 (MQ=255)
                                                                                                      gagaTTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACTATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGAtt                                               >  2:2857077/1‑148 (MQ=255)
                                                                                                        gaTTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGc                                             >  2:3329057/1‑148 (MQ=255)
                                                                                                         aTTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCt                                            >  1:2884624/1‑148 (MQ=255)
                                                                                                          ttGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCtt                                           >  2:606679/1‑148 (MQ=255)
                                                                                                           tGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTc                                          >  1:1548767/1‑148 (MQ=255)
                                                                                                             tCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACTATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAg                                        >  2:569565/1‑148 (MQ=255)
                                                                                                              cTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGt                                       >  2:348685/1‑148 (MQ=255)
                                                                                                                 aTTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGc                                    <  2:3099881/148‑1 (MQ=255)
                                                                                                                 aTTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGc                                    >  1:824334/1‑148 (MQ=255)
                                                                                                                  tttGAAGCCTCCGCCGCTGGTACGGCGGAGGACTATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGcc                                   <  1:1166336/148‑1 (MQ=255)
                                                                                                                  tttGAAGCCTCCGCCGCTGGTACGGCGGAGGACTATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGcc                                   <  1:246281/148‑1 (MQ=255)
                                                                                                                  tttGAAGCCTCCGCCGCTGGTACGGCGGAGGACTATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGcc                                   <  1:2301319/148‑1 (MQ=255)
                                                                                                                     gAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTgg                                <  1:59020/148‑1 (MQ=255)
                                                                                                                       aGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGt                              >  1:3213770/1‑148 (MQ=255)
                                                                                                                        gCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTg                             >  2:3238694/1‑148 (MQ=255)
                                                                                                                         ccTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGc                            >  2:1674185/1‑148 (MQ=255)
                                                                                                                          cTCCGCCGCTGGTACGGCGGAGGACTATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCa                           <  1:2020288/148‑1 (MQ=255)
                                                                                                                           tCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCaa                          >  2:1068031/1‑148 (MQ=255)
                                                                                                                            ccgccgCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCAAt                         <  2:1707865/148‑1 (MQ=255)
                                                                                                                                cgcTGGTACGTCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCACTCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCAATCaaa                     <  2:438285/148‑1 (MQ=255)
                                                                                                                                 gcTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCAATCAAAc                    >  2:2318414/1‑148 (MQ=255)
                                                                                                                                    ggTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCAATCAAACTCa                 <  2:899326/148‑1 (MQ=255)
                                                                                                                                         ggcggAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCAATCAAACTCAATCAc            <  1:2207228/148‑1 (MQ=255)
                                                                                                                                          gcggAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCAATCAAACTCAATCACt           <  1:1978968/148‑1 (MQ=255)
                                                                                                                                          gcggAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCAATCAAACTCAATCACt           >  1:2611068/1‑148 (MQ=255)
                                                                                                                                          gcggAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCAATCAAACTCAATCACt           >  1:2611341/1‑148 (MQ=255)
                                                                                                                                          gcggAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCAATAAAACTCAATCACt           >  2:329214/1‑148 (MQ=255)
                                                                                                                                          gcggAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCAATAAAACTCAATCACt           >  2:1113287/1‑148 (MQ=255)
                                                                                                                                            ggaggaCAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCAATCAAACTCAATCACTaa         <  1:2503334/148‑1 (MQ=255)
                                                                                                                                             gaggaCTATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCAATCAAACTCAATCACTAAc        >  1:2957373/1‑148 (MQ=255)
                                                                                                                                             gaggaCAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCAATCAAACTCAATCACTAAc        >  1:1403348/1‑148 (MQ=255)
                                                                                                                                              aggaCAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCAATCAAACTCAATCACTAACt       <  2:676580/148‑1 (MQ=255)
                                                                                                                                               ggaCAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGAGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCAATCAAACTCAATCACTAACTc      >  2:1161791/1‑148 (MQ=255)
                                                                                                                                                gaCAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCAATCAAACTCAATCACTAACTCa     <  2:543298/148‑1 (MQ=255)
                                                                                                                                                 aCAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCAATCAAACTCAATCACTAACTCAt    >  2:1629746/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TTTTTGCAGGGTTATGCTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCAATCAAACTCAATCACTAACTCATCG  >  CP009273/2077551‑2077845

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: