Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,566,347 C→A 9.3% pseudogene (2612/2937 nt) yjiV → pseudogene; conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,566,3470CA 9.3% 838.7 / 40.0 232pseudogene (2612/2937 nt)yjiVpseudogene; conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (112/98);  new base A (7/15);  total (119/113)
Fisher's exact test for biased strand distribution p-value = 7.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 232 total aligned reads displayed.

AAAACCCGGTTGCGCTGGCGTTGTTGGGATCAATTCTGAGACCGCTAAAAACAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCACCATCGTAAACTGACGGTGACTTTTGAGGATGGTCAGGTATTGAAGATCCGCTTC  >  CP009273/4566200‑4566495
                                                                                                                                                   |                                                                                                                                                    
aaaaCCCGGTTGCGCTGGCGTTGTTGGGATCAATTCTGAGACCGCTAAAAACAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGAgg                                                                                                                                                    <  1:2980067/150‑1 (MQ=255)
  aaCCCGGTTGCGCTGGCGTTGTTGGGATCAATTCTGAGACCGCTAAAAACAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGAt                                                                                                                                                  >  1:449241/1‑150 (MQ=255)
  aaCCCGGTTGCGCTGGCGTTGTTGGGATCAATTCTGAGACCGCTAAAAACAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGAGGAt                                                                                                                                                  <  1:3767263/150‑1 (MQ=255)
    cccGGTTGCGCTGGCGTTGTTGGGATCAATTCTGAGACCGCTAAAAACAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGAttt                                                                                                                                                >  1:3311315/1‑150 (MQ=255)
    cccGGTTGCGCTGGCGTTGTTGGGATCAATTCTGAGACCGCTAAAAACAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGAttt                                                                                                                                                >  1:3312185/1‑150 (MQ=255)
      cGGTTGCGCTGTCGTTGTTGGGATCAATTCTGAGACCGCTAAAAACAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTc                                                                                                                                              <  2:4414709/150‑1 (MQ=255)
      cGGTTGCGCTGGCGTTGTTGGGATCAATTCTGAGACCGCTAAACACAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTc                                                                                                                                              <  2:4414955/150‑1 (MQ=255)
      cGGTTGCGCTGGCGTTGTTGGGATCAATTCTGAGACCGCTAAAAACAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTc                                                                                                                                              <  1:2492896/150‑1 (MQ=255)
      cGGTTGCGCTGGCGTTGTTGGGATCAATTCTGAGACCGCTAAAAACAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTc                                                                                                                                              <  1:2692694/150‑1 (MQ=255)
       ggTTGCGCTGGCGTTGTTGGGATCAATTCTGAGACCGCTAAAAACAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCa                                                                                                                                             <  1:3509057/150‑1 (MQ=255)
         ttGCGCTGGCGTTGTTGGGATCAATTCTGAGACCGCTAAAAACAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAgg                                                                                                                                           >  2:3411284/1‑150 (MQ=255)
            cgcTGGCGTTGTTGGGATCAATTCTGAGACCGCTAAAAACAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGAtt                                                                                                                                        >  1:2133566/1‑150 (MQ=255)
            cgcTGGCGTTGTTGGGATCAATTCTGAGACCGCTAAAAACAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGAtt                                                                                                                                        >  1:4338791/1‑150 (MQ=255)
             gcTGGCGTTGTTGGGATCAATTCTGAGACCGCTAAAAACAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGAttt                                                                                                                                       >  1:867667/1‑150 (MQ=255)
               tGGCGTTGTTGGGATCAATTCTGAGACCGCTAAAAACAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGAGGATTTTCAGGATTTCg                                                                                                                                     <  1:1739468/150‑1 (MQ=255)
                ggCGTTGTTGGGATCAATTCTGAGACCGCTAAAAACAAAGCTGACCGATGGCGCTGAAGTGACTCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGc                                                                                                                                    >  1:2696849/1‑150 (MQ=255)
                ggCGTTGTTGGGATCAATTCTGAGACCGCTAAAAACAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGc                                                                                                                                    >  1:2696728/1‑150 (MQ=255)
                ggCGTTGTTGGGATCAATTCTGAGACCGCTAAAAACAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGc                                                                                                                                    >  1:2696815/1‑150 (MQ=255)
                ggCGTTGTTGGGATCAATTCTGAGACCGCTAAAAACAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGc                                                                                                                                    >  1:407525/1‑150 (MQ=255)
                ggCGTTGTTGGGATCAATTCTGAGACCGCTAAAAACAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGc                                                                                                                                    >  1:2313522/1‑150 (MQ=255)
                            aTCAATTCTGAGACCGCTAAAAACAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGtt                                                                                                                        <  1:712667/150‑1 (MQ=255)
                             tCAATTCTGAGACCGCTAAAAACAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGAGGATTTTCAGGATTTCGCTGATCAGTGGttt                                                                                                                       <  1:3572108/150‑1 (MQ=255)
                              cAATTCTGAGACCGCTAAAAACAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGAGGATTTTCAGGATTTCGCTGATCAGTGGTTtg                                                                                                                      >  1:513236/1‑150 (MQ=255)
                               aaTTCTGAGACCGCTAAAAACAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTtgc                                                                                                                     >  2:784488/1‑150 (MQ=255)
                               aaTTCTGAGACCGCTAAAAACAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTtgc                                                                                                                     >  1:548414/1‑150 (MQ=255)
                                aTTCTGAGACCGCTAAAAACAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTtgct                                                                                                                    >  1:736342/1‑150 (MQ=255)
                                aTTCTGAGACCGCTAAAAACAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTtgct                                                                                                                    >  1:2896296/1‑150 (MQ=255)
                                 ttCAGAGACCGCTAAAAACAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTtgctg                                                                                                                   <  2:2614801/150‑1 (MQ=255)
                                  tCTGAGACCGCTAAAAACAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTtgctgc                                                                                                                  >  1:453667/1‑150 (MQ=255)
                                     gagaCCGCTAAAAACAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTtgctgctgc                                                                                                               <  2:2354379/150‑1 (MQ=255)
                                      agaCCGCTAAAAACAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTgcg                                                                                                              >  1:152324/1‑150 (MQ=255)
                                        aCCGCTAAAAACAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCCTGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCt                                                                                                            <  2:2082092/150‑1 (MQ=255)
                                        aCCGCTAAAAACAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCt                                                                                                            <  2:1824369/150‑1 (MQ=255)
                                        aCCGCTAAAAACAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCt                                                                                                            <  2:1824514/150‑1 (MQ=255)
                                        aCCGCTAAAAACAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCt                                                                                                            <  1:1069549/150‑1 (MQ=255)
                                                   cAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCgg                                                                                                 <  1:724219/150‑1 (MQ=255)
                                                   cAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCgg                                                                                                 >  1:2000132/1‑150 (MQ=255)
                                                    aaaGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGt                                                                                                <  2:2442969/150‑1 (MQ=255)
                                                    aaaGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGt                                                                                                <  1:2393935/150‑1 (MQ=255)
                                                        cTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGaa                                                                                            <  2:3152224/150‑1 (MQ=255)
                                                        cTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGaa                                                                                            <  2:2068146/150‑1 (MQ=255)
                                                                 ggCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGgtg                                                                                   >  1:1104235/1‑150 (MQ=255)
                                                                  gcgcTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGgtgt                                                                                  <  1:118688/150‑1 (MQ=255)
                                                                   cgcTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGtt                                                                                 <  2:1320078/150‑1 (MQ=255)
                                                                   cgcTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGtt                                                                                 <  1:2857425/150‑1 (MQ=255)
                                                                    gcTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTc                                                                                <  1:3459389/150‑1 (MQ=255)
                                                                    gcTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTc                                                                                <  1:3459809/150‑1 (MQ=255)
                                                                       gAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGAc                                                                             <  2:2282360/150‑1 (MQ=255)
                                                                        aaGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACt                                                                            >  2:4179515/1‑150 (MQ=255)
                                                                         aGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTc                                                                           >  1:960541/1‑150 (MQ=255)
                                                                            gACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTcgcc                                                                        <  2:3563935/150‑1 (MQ=255)
                                                                                cTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGAGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCcgcgc                                                                    <  2:255097/150‑1 (MQ=255)
                                                                                 tGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCcgcgcg                                                                   >  1:2476340/1‑150 (MQ=255)
                                                                                  ggATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGa                                                                  >  1:2660334/1‑150 (MQ=255)
                                                                                     tACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGAGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGatat                                                               <  1:4281977/150‑1 (MQ=255)
                                                                                      aCGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATc                                                              >  1:3130508/1‑150 (MQ=255)
                                                                                      aCGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATc                                                              >  1:3985445/1‑150 (MQ=255)
                                                                                      aCGATATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTCATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATc                                                              >  2:2897737/1‑150 (MQ=255)
                                                                                          tATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGGTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTc                                                          >  1:4222201/1‑150 (MQ=255)
                                                                                           aTTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCa                                                         <  2:675640/150‑1 (MQ=255)
                                                                                             ttAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCAcc                                                       >  2:2823883/1‑150 (MQ=255)
                                                                                             ttAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCAcc                                                       >  2:2824122/1‑150 (MQ=255)
                                                                                             ttAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCAcc                                                       <  2:4119307/150‑1 (MQ=255)
                                                                                             ttAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGAGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCAcc                                                       >  2:880714/1‑150 (MQ=255)
                                                                                                 agataaagatCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCACCATCg                                                   <  1:3406809/150‑1 (MQ=255)
                                                                                                    taaagatCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCACCATCGTaa                                                <  2:1730800/150‑1 (MQ=255)
                                                                                                      aagatCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCACCATCGTAAAc                                              >  1:2454836/1‑150 (MQ=255)
                                                                                                       agatCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCACCATCGTAAACt                                             <  1:34943/150‑1 (MQ=255)
                                                                                                        gatCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCACCATCGTAAACTg                                            >  2:2292605/1‑150 (MQ=255)
                                                                                                        gatCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAAATGACGGTGTTCGACTCGCCGCGCGATATCCCTCACCATCGTAAACTg                                            >  1:1909014/1‑150 (MQ=255)
                                                                                                           cGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCACCATCGTAAACTGACg                                         <  1:448576/150‑1 (MQ=255)
                                                                                                           cGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCACCATCGTAAACTGACg                                         <  2:1681604/150‑1 (MQ=255)
                                                                                                           cGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCACCATCGTAAACTGACg                                         >  2:1509576/1‑150 (MQ=255)
                                                                                                              cccGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCACCATCGTAAACTGACGGTg                                      <  2:285519/150‑1 (MQ=255)
                                                                                                                 ggTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCACCATCGTAAACTGACGGTGACt                                   >  2:1163716/1‑150 (MQ=255)
                                                                                                                   tAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATGTTCAGGATTTCGGTGAGCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCTCCGCGCGCTATCCCTCACCATCGTAAACTGACGGTGACttt                                 >  2:2684739/1‑150 (MQ=255)
                                                                                                                    aaCCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCACCATCGTAAACTGACGGTGACtttt                                >  2:3632554/1‑150 (MQ=255)
                                                                                                                    aaCCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCACCATCGTAAACTGACGGTGACtttt                                >  1:1301401/1‑150 (MQ=255)
                                                                                                                    aaCCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCACCATCGTAAACTGACGGTGACtttt                                <  1:1953142/150‑1 (MQ=255)
                                                                                                                       cGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCACCATCGTAAACTGACGGTGACTTTTGAg                             >  1:189616/1‑150 (MQ=255)
                                                                                                                        ggCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCACCATCGTAAACTGACGGTGACTTTTGAgg                            <  1:976703/150‑1 (MQ=255)
                                                                                                                        ggCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCACCATCGTAAACTGACGGTGACTTTTGAgg                            <  2:55131/150‑1 (MQ=255)
                                                                                                                        ggCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCACCATCGTAAACTGACGGTGACTTTTGAgg                            >  1:4473465/1‑150 (MQ=255)
                                                                                                                           cTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCACCATCGTAAACTGACGGTGACTTTTGAGGATg                         >  2:4374760/1‑150 (MQ=255)
                                                                                                                           cTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCACCATCGTAAACTGACGGTGACTTTTGAGGATg                         <  2:1910088/150‑1 (MQ=255)
                                                                                                                               ccATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCGGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCACCATCGTAAACTGACGGTGACTTTTGAGGATGGTCa                     <  2:3960318/150‑1 (MQ=255)
                                                                                                                                 aTGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCACCATCGTAAACTGACGGTGACTTTTGAGGATGGTCAgg                   >  2:404205/1‑150 (MQ=255)
                                                                                                                                 aTGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCACCATCGTAAACTGACGGTGACTTTTGAGGATGGTCAgg                   >  2:1048059/1‑150 (MQ=255)
                                                                                                                                    aCTGGATGAGCATAGAGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGACCTGACGGTGTTCGACTCGCCGCGCGATATCCCTCACCATCGTAAACTGACGGTGACTTTTGAGGATGGTCAGGTAt                <  1:2381829/150‑1 (MQ=255)
                                                                                                                                     cTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCACCATCGTAAACTGACGGTGACTTTTGAGGATGGTCAGGTAtt               <  2:1378886/150‑1 (MQ=255)
                                                                                                                                     cTGGATGAGCAGAGAGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGAGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCACTCACAATCGTAAACTGACGGTGACTTTTGAGGATGGTCAGGTAtt               <  1:2843703/150‑1 (MQ=255)
                                                                                                                                      tGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCACCATCGTAAACTGACGGTGACTTTTGAGGATGGTCAGGTATTg              >  2:487154/1‑150 (MQ=255)
                                                                                                                                          tGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCACCATCGTAAACTGACGGTGACTTTTGAGGATGGTCAGGTATTGAAGa          <  2:736298/150‑1 (MQ=255)
                                                                                                                                          tGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCACCATCGTAAACTGACGGTGACTTTTGAGGATGGTCAGGTATTGAAGa          >  2:1976635/1‑150 (MQ=255)
                                                                                                                                           gAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCACCATCGTAAACTGACGGTGACTTTTGAGGATGGTCAGGTATTGAAGAt         <  1:2729702/150‑1 (MQ=255)
                                                                                                                                             gCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTTGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCACCATCGTAAACTGACGGTGACTTTTGAGGATGGTCAGGTATTGAAGATcc       <  1:1448517/150‑1 (MQ=255)
                                                                                                                                              cATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCACCATCGTAAACTGACGGTGACTTTTGAGGATGGTCAGGTATTGAAGATCCg      >  2:3342548/1‑150 (MQ=255)
                                                                                                                                              cATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCACCATCGTAAACTGACGGTGACTTTTGAGGATGGTCAGGTATTGAAGATCCg      <  2:453667/150‑1 (MQ=255)
                                                                                                                                                tcgCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCACCATCGTAAACTGACGGTGACTTTTGAGGATGGTCAGGTATTGAAGATCCGCt    <  2:3242328/148‑1 (MQ=255)
                                                                                                                                                tAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCACCATCGTAAACTGACGGTGACTTTTGAGGATGGTCAGGTATTGAAGATCCGCt    <  1:2277807/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AAAACCCGGTTGCGCTGGCGTTGTTGGGATCAATTCTGAGACCGCTAAAAACAAAGCTGACCGATGGCGCTGAAGTGACGCTGGATACGCTATTTAAAGATAAAGATCGCCCCGGTAACCGGCCTTTCCATGACTGGATGAGCATAGCGGATTTTCAGGATTTCGCTGATCAGTGGTTTGCTGCTGCGCTGGGCAGGCCGGTAGAACTGACGGTGTTCGACTCGCCGCGCGATATCCCTCACCATCGTAAACTGACGGTGACTTTTGAGGATGGTCAGGTATTGAAGATCCGCTTC  >  CP009273/4566200‑4566495

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: