Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 262,114 C→A 100% D51Y (GAT→TAT)  yafY ← lipoprotein, inner membrane; overproduction stimulates degP expression; CP4‑6 prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273262,1140CA96.2% 670.8 / 14.5 221D51Y (GAT→TAT) yafYlipoprotein, inner membrane; overproduction stimulates degP expression; CP4‑6 prophage
Reads supporting (aligned to +/- strand):  ref base C (4/5);  new base A (100/112);  total (104/117)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 225 total aligned reads displayed.

TCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTA  >  CP009273/261970‑262256
                                                                                                                                                |                                                                                                                                              
tcaGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACatcatc                                                                                                                                             <  1:732835/148‑1 (MQ=255)
 caGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATcc                                                                                                                                            >  2:2536706/1‑148 (MQ=255)
   gTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCtg                                                                                                                                          <  1:3379420/148‑1 (MQ=255)
   gTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCtg                                                                                                                                          <  1:2014969/148‑1 (MQ=255)
   gTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCtg                                                                                                                                          >  1:2734767/1‑148 (MQ=255)
    tGGCTGGACGGTTTTGTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCtgg                                                                                                                                         <  1:329874/148‑1 (MQ=255)
     ggCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCtggt                                                                                                                                        <  2:3087852/148‑1 (MQ=255)
       cTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCtggtgg                                                                                                                                      <  2:2931115/148‑1 (MQ=255)
        tGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGc                                                                                                                                     <  1:1307105/148‑1 (MQ=255)
         ggACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCt                                                                                                                                    >  1:1043111/1‑148 (MQ=255)
             ggTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGAcc                                                                                                                                >  1:2796303/1‑148 (MQ=255)
             ggTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGAcc                                                                                                                                >  1:2796391/1‑148 (MQ=255)
             ggTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGAcc                                                                                                                                >  2:2106439/1‑148 (MQ=255)
                tttttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACcaca                                                                                                                             >  1:1581416/1‑148 (MQ=255)
                tttttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACcaca                                                                                                                             >  2:1000683/1‑148 (MQ=255)
                  tttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATc                                                                                                                           <  2:1969779/148‑1 (MQ=255)
                       tGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGa                                                                                                                      >  2:3172810/1‑148 (MQ=255)
                        gCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGaa                                                                                                                     <  2:759729/148‑1 (MQ=255)
                        gCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGaa                                                                                                                     <  1:78570/148‑1 (MQ=255)
                           tatCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAAcg                                                                                                                  <  1:2760206/148‑1 (MQ=255)
                           tatCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAAcg                                                                                                                  <  1:3190772/148‑1 (MQ=255)
                           tatCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAAcg                                                                                                                  <  2:1341578/148‑1 (MQ=255)
                           tatCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAAcg                                                                                                                  >  1:2097478/1‑148 (MQ=255)
                             tCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAAcgcg                                                                                                                >  2:1140538/1‑148 (MQ=255)
                                   aGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTgggg                                                                                                          <  2:3326086/148‑1 (MQ=255)
                                    ggTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGt                                                                                                         <  2:3241004/148‑1 (MQ=255)
                                      tagtagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGa                                                                                                       >  2:1944876/1‑148 (MQ=255)
                                      tagtagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGa                                                                                                       <  1:2164269/148‑1 (MQ=255)
                                        gtagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAgg                                                                                                     >  2:2563397/1‑148 (MQ=255)
                                         tagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGAAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACAACACTGACACGCGGTGGGGTGAGGt                                                                                                    >  2:3157440/1‑148 (MQ=255)
                                            taACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCt                                                                                                 >  1:937127/1‑148 (MQ=255)
                                                tGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTaaa                                                                                             >  2:3345709/1‑148 (MQ=255)
                                                 gCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGACACGCGGTGGGGTGAGGTCCTTAAAt                                                                                            >  1:1285427/1‑148 (MQ=255)
                                                  ccATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAAtt                                                                                           <  1:3556742/148‑1 (MQ=255)
                                                   cATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAAttt                                                                                          >  1:1185251/1‑148 (MQ=255)
                                                     tCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTAc                                                                                        >  1:2062587/1‑148 (MQ=255)
                                                        tCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGa                                                                                     <  2:847436/148‑1 (MQ=255)
                                                        tCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGa                                                                                     <  2:3517266/148‑1 (MQ=255)
                                                        tCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGa                                                                                     <  2:2633362/148‑1 (MQ=255)
                                                        tCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGa                                                                                     <  1:50727/148‑1 (MQ=255)
                                                        tCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGa                                                                                     <  1:1881050/148‑1 (MQ=255)
                                                         cgccgcCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGTTGAGGTCCTTAAATTTACTGAt                                                                                    <  2:2115683/148‑1 (MQ=255)
                                                         cGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGAt                                                                                    <  1:2407535/148‑1 (MQ=255)
                                                         cGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGAt                                                                                    <  2:2115371/148‑1 (MQ=255)
                                                         cGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGAt                                                                                    >  2:873927/1‑148 (MQ=255)
                                                             gCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGcc                                                                                >  1:798806/1‑148 (MQ=255)
                                                             gCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGcc                                                                                >  1:3246603/1‑148 (MQ=255)
                                                             gCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGCTGGGGTGCGGTCCTTAAATTTACTGATGGcc                                                                                >  2:342211/1‑148 (MQ=255)
                                                               cATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCtt                                                                              <  1:86283/148‑1 (MQ=255)
                                                                aTCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTc                                                                             <  1:387964/148‑1 (MQ=255)
                                                                         cTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGtcatca                                                                    >  1:1450870/1‑148 (MQ=255)
                                                                            gCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCACCCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACtt                                                                 >  2:170307/1‑148 (MQ=255)
                                                                            gCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACtt                                                                 >  1:282596/1‑148 (MQ=255)
                                                                             cGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTc                                                                <  2:1899331/148‑1 (MQ=255)
                                                                              gggTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCt                                                               <  1:473602/148‑1 (MQ=255)
                                                                                gTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGt                                                             >  1:248384/1‑148 (MQ=255)
                                                                                gTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGt                                                             <  2:597765/148‑1 (MQ=255)
                                                                                      aaaaCCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGc                                                       <  1:2701437/148‑1 (MQ=255)
                                                                                      aaaaCCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGc                                                       <  2:591737/148‑1 (MQ=255)
                                                                                       aaaCCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCa                                                      <  1:163499/148‑1 (MQ=255)
                                                                                         aCCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATg                                                    >  2:404348/1‑148 (MQ=255)
                                                                                            tGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGa                                                 <  2:1621423/148‑1 (MQ=255)
                                                                                                  gtgtGACAGGATAAGTAGTCAACCTTTGTCCATTCTTCGCTTAAATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGTGTTCCTTAAATTTACTGATGGCCTTCATGTCATCACTTCTGTCATCGCATGCGATGAGaa                                           <  2:459402/148‑1 (MQ=255)
                                                                                                  gTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGaa                                           <  2:793888/148‑1 (MQ=255)
                                                                                                  gTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGaa                                           <  2:460352/148‑1 (MQ=255)
                                                                                                      gACCGGATAAGTAGTCACCCTGTGACCATTCGTCGCTGACATCAGCCTGGTGGCTGACCACATCACGGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACag                                       <  2:2042060/148‑1 (MQ=255)
                                                                                                      gACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACag                                       <  1:819120/148‑1 (MQ=255)
                                                                                                      gACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACag                                       >  1:1370585/1‑148 (MQ=255)
                                                                                                      gACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACag                                       <  2:1829021/148‑1 (MQ=255)
                                                                                                       aCCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACaga                                      >  2:2484878/1‑148 (MQ=255)
                                                                                                       aCCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACaga                                      <  1:1140538/148‑1 (MQ=255)
                                                                                                       aCCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACaga                                      >  2:250919/1‑148 (MQ=255)
                                                                                                        ccGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACagag                                     <  2:2821269/148‑1 (MQ=255)
                                                                                                         cGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAgt                                    >  1:428592/1‑148 (MQ=255)
                                                                                                         cGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAgt                                    <  1:2582744/148‑1 (MQ=255)
                                                                                                          ggATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAgtg                                   <  2:2764459/148‑1 (MQ=255)
                                                                                                           gATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGAGGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAgtgg                                  <  1:715978/148‑1 (MQ=255)
                                                                                                            aTAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAgtggt                                 >  2:2924584/1‑148 (MQ=255)
                                                                                                            aTAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAgtggt                                 >  2:229682/1‑148 (MQ=255)
                                                                                                            aTAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAgtggt                                 <  1:1000683/148‑1 (MQ=255)
                                                                                                             tAAGTAGTCAACCTGTGACCATTCTTCGCTGACATACTCCTGGTGGCTGACCACATCACTGAAACGCGGTTGGGTGAGGTCCTTAAATTTTCTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAgtggtg                                <  1:2293933/148‑1 (MQ=255)
                                                                                                                  agtCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAAc                           <  1:2375697/148‑1 (MQ=255)
                                                                                                                   gtCAACCTGTGACCATTCTTTGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAAcc                          >  2:1464215/1‑148 (MQ=255)
                                                                                                                        ccTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAgcgc                     <  1:1970080/148‑1 (MQ=255)
                                                                                                                          tgtgACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCca                   >  1:2659978/1‑148 (MQ=255)
                                                                                                                          tgtgACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCca                   >  1:384136/1‑148 (MQ=255)
                                                                                                                            tgACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCcagc                 >  2:2488574/1‑148 (MQ=255)
                                                                                                                             gACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCcagca                >  2:2907821/1‑148 (MQ=255)
                                                                                                                               ccATTCTTCGCTGACATAATCCTGGTGGCTGTCCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCaga              <  1:2488574/148‑1 (MQ=255)
                                                                                                                                  ttcttcGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGc           >  2:648928/1‑148 (MQ=255)
                                                                                                                                  ttcttcGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGc           <  2:1285427/148‑1 (MQ=255)
                                                                                                                                  ttcttcGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGc           <  1:493729/148‑1 (MQ=255)
                                                                                                                                  ttcttcGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGc           <  1:2287785/148‑1 (MQ=255)
                                                                                                                                   tcttcGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCa          <  2:3297746/148‑1 (MQ=255)
                                                                                                                                   tcttcGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCa          >  1:670691/1‑148 (MQ=255)
                                                                                                                                   tcttcGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCa          >  1:2460840/1‑148 (MQ=255)
                                                                                                                                    cttcGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAg         <  2:556373/148‑1 (MQ=255)
                                                                                                                                     ttcGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAgt        <  1:2926638/148‑1 (MQ=255)
                                                                                                                                          tGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCTAGCAGAGGCAGTGtttt   <  2:1581337/148‑1 (MQ=255)
                                                                                                                                          tGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCCGTGtttt   >  2:1147367/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                              
TCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTA  >  CP009273/261970‑262256

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: