Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,913,979 C→A 7.4% A90A (GCG→GCT atpF ← F0 sector of membrane‑bound ATP synthase, subunit b

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,913,9790CA 7.4% 1133.4 / 38.7 301A90A (GCG→GCTatpFF0 sector of membrane‑bound ATP synthase, subunit b
Reads supporting (aligned to +/- strand):  ref base C (153/125);  new base A (9/14);  total (162/139)
Fisher's exact test for biased strand distribution p-value = 1.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 307 total aligned reads displayed.

CGGAACGTTCGATGATCTTCTCGGCGCCAGCAACAGCCAGGATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAG  >  CP009273/3913831‑3914127
                                                                                                                                                    |                                                                                                                                                    
 ggAACGTTCGATGATCTTCTCGGCGCCAGCAACAGCCAGGATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCg                                                                                                                                                     >  2:3770151/1‑149 (MQ=255)
   aaCGTTCGATGATCTTCTCGGCGCCAGCAACAGCCAGGATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTAGCt                                                                                                                                                   >  1:3004395/1‑149 (MQ=255)
       ttCGATGATCTTCTCGGCGCCAGCAACAGCCAGGATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGt                                                                                                                                               >  1:417354/1‑149 (MQ=255)
        tCGATGATCTTCTCGGCGCCAGCAACAGCCAGGATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTc                                                                                                                                              <  1:91092/149‑1 (MQ=255)
           atgatCTTCTCGGCGCCAGCAACAGCCAGGATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAg                                                                                                                                           >  2:1024242/1‑149 (MQ=255)
           atgatCTTCTCGGCGCCAGCAACAGCCAGGATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTAGCTTCGTCCAg                                                                                                                                           >  2:5748/1‑149 (MQ=255)
            tgatCTTCTCGGCGCCAGCAACAGCCAGGATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGa                                                                                                                                          >  1:4414446/1‑149 (MQ=255)
             gatCTTCTCGGCGCCAGCAACAGCCAGGATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTAGCTTCGTCCAGaa                                                                                                                                         <  1:3307257/149‑1 (MQ=255)
               tcttctCGGCGCCAGCAACAGCCAGGATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTAGCTTCGTCCAGAATc                                                                                                                                       <  1:108016/149‑1 (MQ=255)
                      ggCGCCAGCAACAGCCAGGATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAgc                                                                                                                                <  2:1356421/149‑1 (MQ=255)
                       gcgcCAGCAACAGCCAGGATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAgcg                                                                                                                               >  1:2658861/1‑149 (MQ=255)
                       gcgcCAGCAACAGCCAGGATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTGTCGCTTCGTCCAGAATCTGCGAgcg                                                                                                                               >  1:1233050/1‑149 (MQ=255)
                          ccAGCAACAGCCAGGATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAgcggcg                                                                                                                            >  1:4252911/1‑149 (MQ=255)
                            aGCAACAGCCAGGATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGtt                                                                                                                          <  2:2435748/149‑1 (MQ=255)
                            aGCAACAGCCAGGATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGtt                                                                                                                          <  1:476945/149‑1 (MQ=255)
                            aGCAACAGCCAGGATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGtt                                                                                                                          >  2:1315665/1‑149 (MQ=255)
                              cAACAGCCAGGATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTg                                                                                                                        <  2:1471908/149‑1 (MQ=255)
                               aaCAGCCAGGATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGt                                                                                                                       <  1:4222977/149‑1 (MQ=255)
                                aCAGCCAGGATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGtt                                                                                                                      <  2:164265/149‑1 (MQ=255)
                                  aGCCAGGATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCg                                                                                                                    >  2:1286248/1‑149 (MQ=255)
                                  aGCCAGGATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCg                                                                                                                    >  2:1679399/1‑149 (MQ=255)
                                   gCCAGGATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGc                                                                                                                   >  1:3226721/1‑149 (MQ=255)
                                   gCCAGGATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATATCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGc                                                                                                                   >  2:160448/1‑149 (MQ=255)
                                       ggATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGc                                                                                                               >  1:1135577/1‑149 (MQ=255)
                                        gATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCt                                                                                                              <  1:1897188/149‑1 (MQ=255)
                                        gATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCt                                                                                                              >  2:615198/1‑149 (MQ=255)
                                         aTAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTc                                                                                                             >  2:3015351/1‑149 (MQ=255)
                                         aTAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTc                                                                                                             <  2:488092/149‑1 (MQ=255)
                                               aCTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCgatgat                                                                                                       >  1:2396278/1‑149 (MQ=255)
                                               aCTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCgatgat                                                                                                       <  1:2219838/149‑1 (MQ=255)
                                               aCTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTAGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCgatgat                                                                                                       <  2:2276239/149‑1 (MQ=255)
                                                cTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGAtt                                                                                                      >  2:2699774/1‑149 (MQ=255)
                                                 ttGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTa                                                                                                     >  1:811683/1‑149 (MQ=255)
                                                 ttGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTa                                                                                                     >  1:2171478/1‑149 (MQ=255)
                                                 ttGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTa                                                                                                     >  1:3798233/1‑149 (MQ=255)
                                                 gtGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTa                                                                                                     <  2:4080601/148‑1 (MQ=255)
                                                  tGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTAc                                                                                                    <  1:1360928/149‑1 (MQ=255)
                                                     ttACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTg                                                                                                 >  1:2375694/1‑149 (MQ=255)
                                                     ttACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTAGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTg                                                                                                 <  1:3998237/149‑1 (MQ=255)
                                                        cGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGc                                                                                              <  1:1198534/149‑1 (MQ=255)
                                                           aGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTc                                                                                           <  2:4679901/149‑1 (MQ=255)
                                                           aGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTc                                                                                           <  2:3234633/149‑1 (MQ=255)
                                                                 tCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCcgcttt                                                                                     >  2:4242349/1‑149 (MQ=255)
                                                                 tCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCcgcttt                                                                                     <  1:3104649/149‑1 (MQ=255)
                                                                  cACGGGCACGTGTACGCTCGGCTTCAAGTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCcgctttc                                                                                    <  2:1177325/149‑1 (MQ=255)
                                                                   aCGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCcgctttcg                                                                                   >  2:256179/1‑149 (MQ=255)
                                                                     gggCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTATCTGGGCTTCcgctttcgct                                                                                 <  2:4644930/149‑1 (MQ=255)
                                                                      ggCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCcgctttcgctt                                                                                >  1:4487785/1‑149 (MQ=255)
                                                                      ggCACGTTTACGCTCGGCTACAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCAGGGCATCcgcattcgctt                                                                                >  1:3560671/1‑149 (MQ=255)
                                                                       gCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTAGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCcgctttcgcttt                                                                               >  2:1385012/1‑149 (MQ=255)
                                                                        cACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCtttt                                                                              >  1:4741345/1‑149 (MQ=255)
                                                                          cGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCtttttt                                                                            >  1:3080259/1‑149 (MQ=255)
                                                                           gTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTc                                                                           >  2:1545904/1‑149 (MQ=255)
                                                                            tttACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCa                                                                          <  1:684123/149‑1 (MQ=255)
                                                                               aCGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCt                                                                       <  2:4252911/149‑1 (MQ=255)
                                                                                   tCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTggtc                                                                   >  1:4655338/1‑149 (MQ=255)
                                                                                   tCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTggtc                                                                   >  2:3049841/1‑149 (MQ=255)
                                                                                      gCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTggtcggt                                                                <  1:1577328/149‑1 (MQ=255)
                                                                                      gCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTggtcggt                                                                <  1:3265039/149‑1 (MQ=255)
                                                                                       cTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTAGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTggtcggtc                                                               <  2:4459132/149‑1 (MQ=255)
                                                                                            aTTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCt                                                          >  2:3516782/1‑149 (MQ=255)
                                                                                             tgtCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTg                                                         <  2:3683633/147‑1 (MQ=255)
                                                                                              ttCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTgg                                                        <  2:3623156/149‑1 (MQ=255)
                                                                                                ccGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGcc                                                      <  1:1472350/149‑1 (MQ=255)
                                                                                                 cgcctgcgcctgGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCt                                                     <  1:736584/149‑1 (MQ=255)
                                                                                                   cctgcgcctgGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCttt                                                   <  2:2153103/149‑1 (MQ=255)
                                                                                                   cctgcgcctgGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTAGCCTGCTCGTTGATTACCTGGGCTTCCGGTTTCGCTGTTTTCATCTGGTCGGTCGCGCTGGcggtt                                                   >  1:2948599/1‑145 (MQ=255)
                                                                                                    ctgcgcctgGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTg                                                  <  1:3457336/149‑1 (MQ=255)
                                                                                                         cctgGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGcaagg                                             >  2:869230/1‑149 (MQ=255)
                                                                                                           tgGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGcaaggtc                                           <  1:3698668/149‑1 (MQ=255)
                                                                                                           tgGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTAGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCATCTGGTCGGTCGCGCTGGCCGTTTcaaggtc                                           >  1:2814241/1‑149 (MQ=255)
                                                                                                             ggCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGcaaggtcaa                                         >  1:4188849/1‑149 (MQ=255)
                                                                                                               ccACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTAACTGGGCTTCCGCTTTCTCTTTTTTAATCTGGTCGGTCGCGCTGGCCGTTtcaaggtcaagg                                       >  1:1234304/1‑149 (MQ=255)
                                                                                                                 aCGATTGTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCaaggtcaaggtc                                     <  1:4013312/149‑1 (MQ=255)
                                                                                                                   gATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTAAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCt                                   >  1:1136760/1‑149 (MQ=255)
                                                                                                                   gATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTAAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCt                                   >  1:3524970/1‑149 (MQ=255)
                                                                                                                    aTTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAATGTCCtt                                  >  1:261547/1‑149 (MQ=255)
                                                                                                                    aTTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCtt                                  >  2:4212287/1‑149 (MQ=255)
                                                                                                                    aTTTTAGTACGGTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCtt                                  <  2:3085542/149‑1 (MQ=255)
                                                                                                                         aGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTAtgtg                             >  1:708651/1‑149 (MQ=255)
                                                                                                                              gTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTAGCCTGCTCGATGATTAACTGGGCTTACGCTTTCGCTTTTTTCATCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTTCTCGt                        >  1:1669149/1‑149 (MQ=255)
                                                                                                                                tCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTc                      <  2:3949274/149‑1 (MQ=255)
                                                                                                                                 ccTGTTCTGCCTCAGCTTTAGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCt                     <  2:2933121/149‑1 (MQ=255)
                                                                                                                                    gTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTAAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCg                  >  2:4012877/1‑149 (MQ=255)
                                                                                                                                     ttCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCgg                 >  1:3984934/1‑149 (MQ=255)
                                                                                                                                     ttCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCgg                 >  2:1726583/1‑149 (MQ=255)
                                                                                                                                     ttCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCATTTTTGATCGTGTAGGTGGCGCAGGCCATTGCAAGGTAAAGGGCCGTATGTTGTCGTTATGCgg                 >  1:2135075/1‑149 (MQ=255)
                                                                                                                                      tCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGa                <  2:3552857/149‑1 (MQ=255)
                                                                                                                                       cTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGAGGCTTACCTGGGCTTCCGCTTTCGCTTGTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGaa               >  2:3300366/1‑149 (MQ=255)
                                                                                                                                        ggCCTCAGCGTGCGCTACGTCCAGAATCTGCGAGCGGCGTTTGGTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGCTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAg              <  2:1099715/148‑1 (MQ=255)
                                                                                                                                          ccTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCa            <  1:1849321/149‑1 (MQ=255)
                                                                                                                                          ccTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCa            <  1:2739155/149‑1 (MQ=255)
                                                                                                                                          ccTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTCTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCa            <  2:3350336/149‑1 (MQ=255)
                                                                                                                                           cTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCaa           >  2:3529684/1‑149 (MQ=255)
                                                                                                                                           cTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCaa           >  2:1827220/1‑149 (MQ=255)
                                                                                                                                            tCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAg          <  2:1957395/149‑1 (MQ=255)
                                                                                                                                              aGCTTTCGCTTCGTCCAGACTCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGc        <  1:1736976/149‑1 (MQ=255)
                                                                                                                                                  ttCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTc    >  2:4718119/1‑149 (MQ=255)
                                                                                                                                                    cGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAg  <  2:3398836/149‑1 (MQ=255)
                                                                                                                                                    cGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAg  >  1:3919281/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGGAACGTTCGATGATCTTCTCGGCGCCAGCAACAGCCAGGATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAG  >  CP009273/3913831‑3914127

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: