Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 430,523 G→A 10.3% A141T (GCA→ACA)  ribE → riboflavin synthase beta chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273430,5230GA10.3% 727.1 / 51.1 223A141T (GCA→ACA) ribEriboflavin synthase beta chain
Reads supporting (aligned to +/- strand):  ref base G (105/95);  new base A (7/16);  total (112/111)
Fisher's exact test for biased strand distribution p-value = 5.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 225 total aligned reads displayed.

TATGTCGCTGGTGGTGCAAGCAACGGCCTGGCGCATGTTGCCCAGGACAGCGAAATTCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTGAGTGTGCCGTCCAGGCGCTCTACTCCTGGCAGTTGTCCCAGAACG  >  CP009273/430376‑430668
                                                                                                                                                   |                                                                                                                                                 
 atGTCGCTGGTGGTGCAAGCAACGGCCTGGCGCATGTTGCCCAGGACAGCGAAATTCCGGTTGCTTTTGGGGTTCTGATCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGc                                                                                                                                                  >  2:3432880/1‑148 (MQ=255)
  tGTCGCTGGTGGTGCAAGCAACGGCCTGGCGCATGTTGCCCAGGACAGCGAAATTCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCa                                                                                                                                                 >  2:3284070/1‑148 (MQ=255)
      gCTGGTGGTGCAAGCAACGGCCTGGCGCATGTTGCCCAGGACAGCGAAATTCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGa                                                                                                                                             <  2:493598/148‑1 (MQ=255)
        tggtggTGCAAGCAACGGCCTGGCGCATGTTGCCCAGGACAGCGAAATTCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGAcc                                                                                                                                           >  1:2752205/1‑148 (MQ=255)
         ggtggtGCAAGCAACGGCCTGGCGCATGTTGCCCAGGACAGCGAAATTCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACcg                                                                                                                                          >  1:3095940/1‑148 (MQ=255)
         ggtggtGCAAGCAACGGCCTGGCGCATGTTGCCCAGGACAGCGAAATTCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACcg                                                                                                                                          >  1:1955417/1‑148 (MQ=255)
           tggtgCAAGCAACGGCCTGGCGCATGTTGCCCAGGACAGCGAAATTCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACcgcg                                                                                                                                        <  1:2147243/148‑1 (MQ=255)
            ggtgCAAGCAACGGCCTGGCGCATGTTGCCCAGGACAGCGAAATTCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCgcgc                                                                                                                                       <  1:2837137/148‑1 (MQ=255)
             gtgCAAGCAACGGCCTGGCGCATGTTGCCCAGGACAGCGAAATTCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTACACTGACCGCGCt                                                                                                                                      >  2:397152/1‑148 (MQ=255)
                    caaCGGCCTGGCGCATGTTGCCCAGGACAGCGAAATTCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATg                                                                                                                               >  1:1017838/1‑148 (MQ=255)
                    caaCGGCCTGGCGCATGTTGCCCAGGACAGCGAAATTCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTACACTGACCGCGCTTGAAATg                                                                                                                               <  2:728251/148‑1 (MQ=255)
                    caaCGGCCTGGCGCATGTTGCCCAGGACAGCGAAATTCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTACACTGACCGCGCTTGAAATg                                                                                                                               <  2:2776699/148‑1 (MQ=255)
                      aCGGCCTGGCGCATGTTGCCCAGGACAGCGAAATTCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGAt                                                                                                                             >  2:1996254/1‑148 (MQ=255)
                      aCGGCCTGGCGCATGTTGCCCAGGACAGCGAAATTCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGAt                                                                                                                             <  2:1870523/148‑1 (MQ=255)
                       cGGCCTGGCGCATGTTGCCCAGGACAGCGAAATTCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGAtt                                                                                                                            >  2:939797/1‑148 (MQ=255)
                       cGGCCTGGCGCATGTTGCCCAGGACAGCGAAATTCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGAtt                                                                                                                            >  2:2253573/1‑148 (MQ=255)
                        ggCCTGGCGCATGGTGCCCAGGACAGCGAAATTCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTa                                                                                                                           <  2:2113471/148‑1 (MQ=255)
                            tGGCGCATGTTGCCCAGGACAGCGAAATTCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGt                                                                                                                       >  1:957361/1‑148 (MQ=255)
                            tGGCGCATGTTGCCCAGGACAGCGAAATTCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGt                                                                                                                       >  1:2479095/1‑148 (MQ=255)
                                 cATGTTGCCCAGGACCGCGAAATTCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGa                                                                                                                  >  2:37651/1‑148 (MQ=255)
                                   tGTTGCCCAGGACAGCGAAATTCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTACACTGACCGCGCTTGAAATGATTAATGTATTGaaa                                                                                                                <  2:3201142/148‑1 (MQ=255)
                                     ttGCCCAGGACAGCGAAATTCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTACACTGACCGCGCTTGAAATGATTAATGTATTGAAAGc                                                                                                              <  1:627931/148‑1 (MQ=255)
                                        cccAGGACAGCGAAATTCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCAt                                                                                                           <  2:2919036/148‑1 (MQ=255)
                                           aGGACAGCGAAATTCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCaa                                                                                                        <  2:872212/148‑1 (MQ=255)
                                             gACAGCGAAATTCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAgg                                                                                                      >  1:3449214/1‑148 (MQ=255)
                                                  cGAAATTCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGa                                                                                                 <  2:1223064/148‑1 (MQ=255)
                                                  cGAAATTCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGa                                                                                                 >  2:3448807/1‑148 (MQ=255)
                                                  cGAAATTCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAACGTGCAGAAGCTGCACTGACCGCGCTTGAAATCATTAATGTATTGAAAGCCATCAAGGCCTGa                                                                                                 >  2:2060774/1‑148 (MQ=255)
                                                   gAAATTCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGaa                                                                                                <  2:3495242/148‑1 (MQ=255)
                                                   gAAATTCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGaa                                                                                                <  2:2271995/148‑1 (MQ=255)
                                                   gAAATTCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGaa                                                                                                <  2:2359162/148‑1 (MQ=255)
                                                    aaaTTCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGaaa                                                                                               >  2:2550487/1‑148 (MQ=255)
                                                    aaaTTCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGaaa                                                                                               >  1:2555632/1‑148 (MQ=255)
                                                     aaTTCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAAt                                                                                              >  2:187902/1‑148 (MQ=255)
                                                      aTTCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAAtt                                                                                             <  2:1367522/148‑1 (MQ=255)
                                                        tCCGGTTGCTTTTTGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCCCGCTTGAAATGATTAATTTATTTAAAGCCAACAAGGCCTGAAATTAg                                                                                           >  1:1466008/1‑148 (MQ=255)
                                                        tCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAg                                                                                           >  2:273169/1‑148 (MQ=255)
                                                        tCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAg                                                                                           <  1:1996254/148‑1 (MQ=255)
                                                        tCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAACTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAg                                                                                           <  1:3496917/148‑1 (MQ=255)
                                                         ccGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGt                                                                                          >  2:3061322/1‑148 (MQ=255)
                                                          cGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTa                                                                                         <  2:461697/148‑1 (MQ=255)
                                                             ttGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAgg                                                                                      >  1:895257/1‑148 (MQ=255)
                                                              tGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAggg                                                                                     <  1:2786391/148‑1 (MQ=255)
                                                                  tttGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGaaa                                                                                 <  1:1456054/148‑1 (MQ=255)
                                                                      gggTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCg                                                                             >  2:2037676/1‑148 (MQ=255)
                                                                           cTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGaaa                                                                        <  1:711932/148‑1 (MQ=255)
                                                                            tGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAAc                                                                       <  2:3160898/148‑1 (MQ=255)
                                                                              aCCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACct                                                                     >  1:2074065/1‑148 (MQ=255)
                                                                               ccACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACctg                                                                    >  1:1006227/1‑148 (MQ=255)
                                                                               ccACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACctg                                                                    <  1:674764/148‑1 (MQ=255)
                                                                                  cTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACctgctg                                                                 <  2:764981/148‑1 (MQ=255)
                                                                                  cTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACctgctg                                                                 <  2:879037/148‑1 (MQ=255)
                                                                                  cTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACctgctg                                                                 >  2:530784/1‑148 (MQ=255)
                                                                                  cTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACctgctg                                                                 <  2:879100/148‑1 (MQ=255)
                                                                                      aaGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCtcg                                                             <  2:668287/148‑1 (MQ=255)
                                                                                      aaGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTACACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCtcg                                                             <  2:2687746/148‑1 (MQ=255)
                                                                                       aGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCtcgt                                                            <  1:2686028/148‑1 (MQ=255)
                                                                                        gCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCtcgtc                                                           >  2:3510751/1‑148 (MQ=255)
                                                                                        gCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAAGTAGAAAGGGGAAATCCGTGAAACCTGCTGCtcgtc                                                           >  1:978536/1‑148 (MQ=255)
                                                                                           ttGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTcgcc                                                        >  1:1486574/1‑148 (MQ=255)
                                                                                            tGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTcgccg                                                       >  2:1280040/1‑148 (MQ=255)
                                                                                            tGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTcgccg                                                       >  2:2095871/1‑148 (MQ=255)
                                                                                              aaCAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTACACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCcgcg                                                     <  2:307221/148‑1 (MQ=255)
                                                                                                 aaGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTc                                                  >  2:1715019/1‑148 (MQ=255)
                                                                                                    cGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTg                                               >  1:2428259/1‑148 (MQ=255)
                                                                                                    cGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTg                                               >  1:814209/1‑148 (MQ=255)
                                                                                                    cGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTg                                               >  1:3151114/1‑148 (MQ=255)
                                                                                                       tCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTGAgt                                            <  2:19081/148‑1 (MQ=255)
                                                                                                          aaCGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTGAGtgtg                                         <  2:2202187/148‑1 (MQ=255)
                                                                                                           aCGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTGAGTGTGc                                        >  1:1424383/1‑148 (MQ=255)
                                                                                                              tGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTGAGTGTGCCGt                                     <  2:581165/148‑1 (MQ=255)
                                                                                                               gCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTGAGTGTGCCGTc                                    >  1:3253445/1‑148 (MQ=255)
                                                                                                               gCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTGAGTGTGCCGTc                                    <  2:2752205/148‑1 (MQ=255)
                                                                                                                  ggCACCAAAGCTGGCAACAAAGGTGCAGAAGCTACACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTGAGTGTGCCGTCCAg                                 <  2:1544726/148‑1 (MQ=255)
                                                                                                                      ccAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTGAGTGTGCCGTCCAGgcgc                             <  1:1230715/148‑1 (MQ=255)
                                                                                                                            cTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTGAGTGTGCCGTCCAGGCGCTCTACt                       <  1:2233558/148‑1 (MQ=255)
                                                                                                                             tGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTGAGTGTGCCGTCCAGGCGCTCTACTc                      <  2:2589717/148‑1 (MQ=255)
                                                                                                                             tGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTGAGTGTGCCGTCCAGGCGCTCTACTc                      <  1:1061240/148‑1 (MQ=255)
                                                                                                                               gCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTGAGTGTGCCGTCCAGGCGCTCTACTCCt                    <  2:1028026/148‑1 (MQ=255)
                                                                                                                               gCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTGAGTGTGCCGTCCAGGCGCTCTACTCCt                    <  2:1028008/148‑1 (MQ=255)
                                                                                                                               gCAACAAAGGTGCAGAAGCTACACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTGAGTGTGCCGTCCAGGCGCTCTACTCCt                    <  1:2439298/148‑1 (MQ=255)
                                                                                                                                  acaaAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTGAGTGTGCCGTCCAGGCGCTCTACTCCTGGc                 <  1:2809089/148‑1 (MQ=255)
                                                                                                                                  acaaAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTGAGTGTGCCGTCCAGGCGCTCTACTCCTGGc                 >  1:3198482/1‑148 (MQ=255)
                                                                                                                                   caaAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTGAGTGTGCCGTCCAGGCGCTCTACTCCTGGCa                >  2:1688211/1‑148 (MQ=255)
                                                                                                                                    aaaGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTGAGTGTGCCGTCCAGGCGCTCTACTCCTGGCAg               <  1:1275060/148‑1 (MQ=255)
                                                                                                                                      aGGTGCAGAAGCTGCACTGACCGCGCTTGAAATCATTAATGTATTGAAAGCCATCAATGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTGAGTGTGCCGTCCAGGCGCTCTACTCCTGGCAGtt             >  1:3034065/1‑148 (MQ=255)
                                                                                                                                      aGGTGCAGAAGCTACACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTGAGTGTGCCGTCCAGGCGCTCTACTCCTGGCAGtt             >  1:225756/1‑148 (MQ=255)
                                                                                                                                         tGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTGAGTGTGCCGTCCAGGCGCTCTACTCCTGGCAGTTGTc          >  1:1096524/1‑148 (MQ=255)
                                                                                                                                           cAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTGAGTGTGCCGTCCAGGCGCTCTACTCCTGGCAGTTGTccc        <  1:273169/148‑1 (MQ=255)
                                                                                                                                           cAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTGAGTGTGCCGTCCAGGCGCTCTACTCCTGGCAGTTGTccc        <  2:2636649/148‑1 (MQ=255)
                                                                                                                                             gAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTGAGTGTGCCGTCCAGGCGCTCTACTCCTGGCAGTTGTCCCAg      >  1:3495938/1‑148 (MQ=255)
                                                                                                                                              aaGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTGAGTGTGCCGTCCAGGCGCTCTACTCCTGGCAGTTGTCCCAGa     <  1:3114548/148‑1 (MQ=255)
                                                                                                                                              aaGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTGAGTGTGCCGTCCAGGCGCTCTACTCCTGGCAGTTGTCCCAGa     >  2:3031156/1‑148 (MQ=255)
                                                                                                                                              aaGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTGAGTGTGCCGTCCAGGCGCTCTACTCCTGGCAGTTGTCCCAGa     >  2:79486/1‑148 (MQ=255)
                                                                                                                                               aGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTGAGTGTGCCGTCCAGGCGCTCTACTCCTGGCAGTTGTCCCAGaa    >  1:1734482/1‑148 (MQ=255)
                                                                                                                                                gCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTGAGTGTGCCGTCCAGGCGCTCTACTCCTGGCAGTTGTCCCAGAAc   >  1:136994/1‑148 (MQ=255)
                                                                                                                                                gCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTGAGTGTGCCGTCCAGGCGCTCTACTCCTGGCAGTTGTCCCAGAAc   >  1:1029513/1‑148 (MQ=255)
                                                                                                                                                gCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTGAGTGTGCCGTCCAGGCGCTCTACTCCTGGCAGTTGTCCCAGAAc   >  1:1699231/1‑148 (MQ=255)
                                                                                                                                                 cTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTGAGTGTGCCGTCCAGGCGCTCTACTCCTGGCAGTTGTCCCAGAACg  <  2:130394/148‑1 (MQ=255)
                                                                                                                                                 cTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTGAGTGTGCCGTCCAGGCGCTCTACTCCTGGCAGTTGTCCCAGAACg  <  2:129458/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TATGTCGCTGGTGGTGCAAGCAACGGCCTGGCGCATGTTGCCCAGGACAGCGAAATTCCGGTTGCTTTTGGGGTTCTGACCACTGAAAGCATTGAACAAGCGATCGAACGTGCTGGCACCAAAGCTGGCAACAAAGGTGCAGAAGCTGCACTGACCGCGCTTGAAATGATTAATGTATTGAAAGCCATCAAGGCCTGAAATTAGTAAGGGGAAATCCGTGAAACCTGCTGCTCGTCGCCGCGCTCGTGAGTGTGCCGTCCAGGCGCTCTACTCCTGGCAGTTGTCCCAGAACG  >  CP009273/430376‑430668

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: