Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,404,169 G→T 100% A26A (GCG→GCT aidB → DNA alkylation damage repair protein; flavin‑containing DNA binding protein, weak isovaleryl CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,404,1690GT97.6% 1057.3 / 8.9 331A26A (GCG→GCTaidBDNA alkylation damage repair protein; flavin‑containing DNA binding protein, weak isovaleryl CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (6/2);  new base T (162/161);  total (168/163)
Fisher's exact test for biased strand distribution p-value = 2.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 339 total aligned reads displayed.

AATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGATGCGCA  >  CP009273/4404021‑4404315
                                                                                                                                                    |                                                                                                                                                  
aaTCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGgcgcgc                                                                                                                                                   >  1:116477/1‑150 (MQ=255)
 aTCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGctct                                                                                                                                                  >  2:2482094/1‑150 (MQ=255)
  tctcTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCtctc                                                                                                                                                 <  1:3444599/150‑1 (MQ=255)
      tttCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTTCACTGGCAAACTCACACCGTTTTTAATCAACATATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCg                                                                                                                                             >  2:4750120/1‑150 (MQ=255)
      tttCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCg                                                                                                                                             >  2:4123620/1‑150 (MQ=255)
        tCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGaa                                                                                                                                           >  2:1814656/1‑150 (MQ=255)
        tCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGaa                                                                                                                                           >  1:4261532/1‑150 (MQ=255)
          tGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGc                                                                                                                                         <  2:4238845/150‑1 (MQ=255)
          tGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGc                                                                                                                                         >  1:1764490/1‑150 (MQ=255)
          tGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGc                                                                                                                                         <  2:1102009/150‑1 (MQ=255)
           gTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCg                                                                                                                                        <  1:4959232/150‑1 (MQ=255)
           gTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCg                                                                                                                                        <  1:1941068/150‑1 (MQ=255)
              aTGAATCCATGGCAGTGACCCTACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACGGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTa                                                                                                                                     <  2:3242974/150‑1 (MQ=255)
              aTGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACAATTAAATAACAGCAACTTATAACTGTCTGATGGCGCTCTCTGCGAAGCGGTa                                                                                                                                     >  1:932971/1‑150 (MQ=255)
               tGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTaa                                                                                                                                    >  2:1931864/1‑150 (MQ=255)
                     cATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTg                                                                                                                              <  2:1186804/150‑1 (MQ=255)
                      aTGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGa                                                                                                                             >  2:1294867/1‑150 (MQ=255)
                      aTGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGa                                                                                                                             <  2:1712174/150‑1 (MQ=255)
                      aTGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGa                                                                                                                             <  2:1124351/150‑1 (MQ=255)
                        ggCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAg                                                                                                                           <  1:5154525/150‑1 (MQ=255)
                        ggCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACGGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAg                                                                                                                           <  1:2357571/150‑1 (MQ=255)
                          cAGTGACCATACTAATGGTGACTGCCATTGGGGGTGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGt                                                                                                                         >  1:1737731/1‑150 (MQ=255)
                          cAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGt                                                                                                                         <  2:3410053/150‑1 (MQ=255)
                           aGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTg                                                                                                                        >  1:1018809/1‑150 (MQ=255)
                           aGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTg                                                                                                                        <  1:3573569/150‑1 (MQ=255)
                            gTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTgc                                                                                                                       <  1:2133059/150‑1 (MQ=255)
                            gTGACCATACTAAGGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTAGACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTgc                                                                                                                       <  2:951032/150‑1 (MQ=255)
                             tGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTgct                                                                                                                      <  1:4547148/150‑1 (MQ=255)
                              gACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTgctg                                                                                                                     <  1:4872586/150‑1 (MQ=255)
                              gACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTgctg                                                                                                                     <  2:4507930/150‑1 (MQ=255)
                                ccATACTTATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATCCCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTgctggc                                                                                                                   <  1:3104299/150‑1 (MQ=255)
                                 cATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTgctggct                                                                                                                  <  2:3115551/150‑1 (MQ=255)
                                 cATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTgctggct                                                                                                                  >  2:1387189/1‑150 (MQ=255)
                                  aTACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTgctggctg                                                                                                                 <  2:1266376/150‑1 (MQ=255)
                                   tACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGctggctgg                                                                                                                >  2:266343/1‑150 (MQ=255)
                                       aaTGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATa                                                                                                            <  1:2155945/150‑1 (MQ=255)
                                                  ccATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTa                                                                                                 >  2:2852213/1‑150 (MQ=255)
                                                    aTTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGc                                                                                               >  1:373557/1‑150 (MQ=255)
                                                    aTTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGc                                                                                               >  2:4814644/1‑150 (MQ=255)
                                                    aTTGATGGAGGGAGACACAGTGCACGGGCAAACTCACACCGTTTTTAATCAACCTATACAATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGAGAAGCGGTAACGCGTGAAGGTGCTGTCTGGGATAGCGATTTTCTAGc                                                                                               <  1:956299/150‑1 (MQ=255)
                                                       gATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAg                                                                                            >  2:1603222/1‑150 (MQ=255)
                                                        aTGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGt                                                                                           >  1:4540126/1‑150 (MQ=255)
                                                         tGGAGGGAGACACAGTGCACTGGCAACCTCACACCGTTTTTAATCAACCTATCCCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTa                                                                                          <  2:3054956/150‑1 (MQ=255)
                                                          ggagggagACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTAt                                                                                         <  1:203247/150‑1 (MQ=255)
                                                           gagggagACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTAtt                                                                                        >  1:2620447/1‑150 (MQ=255)
                                                             gggagACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTAGACCCTGAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTgg                                                                                      <  2:3919905/150‑1 (MQ=255)
                                                                agaCACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTca                                                                                   >  1:5034593/1‑150 (MQ=255)
                                                                 gaCACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTcag                                                                                  >  1:4802294/1‑150 (MQ=255)
                                                                 gaCACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTcag                                                                                  <  1:1104135/150‑1 (MQ=255)
                                                                   cacaGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCGGGGATCGCGATTTTCTAGCCAGTATTGGTcagca                                                                                >  2:1458708/1‑150 (MQ=255)
                                                                     caGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGt                                                                              >  2:2113465/1‑150 (MQ=255)
                                                                         gCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAgg                                                                          <  1:3170139/150‑1 (MQ=255)
                                                                          cACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGa                                                                         >  2:3204929/1‑150 (MQ=255)
                                                                          cACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGa                                                                         <  2:63452/150‑1 (MQ=255)
                                                                             tGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACg                                                                      >  1:1800796/1‑150 (MQ=255)
                                                                               gCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGc                                                                    >  1:1091411/1‑150 (MQ=255)
                                                                                cAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCt                                                                   >  2:1095377/1‑150 (MQ=255)
                                                                                 aaaCTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTg                                                                  >  1:3833216/1‑150 (MQ=255)
                                                                                 aaaCTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTg                                                                  >  1:3358856/1‑150 (MQ=255)
                                                                                 aaaCTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTg                                                                  <  2:2607089/150‑1 (MQ=255)
                                                                                  aaCTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGa                                                                 >  2:4435608/1‑150 (MQ=255)
                                                                                  aaCTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGa                                                                 >  1:671945/1‑150 (MQ=255)
                                                                                  aaCTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGa                                                                 >  1:4782108/1‑150 (MQ=255)
                                                                                        cacCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCt                                                           >  1:4408929/1‑150 (MQ=255)
                                                                                        cacCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCt                                                           >  1:3269492/1‑150 (MQ=255)
                                                                                               tttAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTg                                                    <  2:858473/150‑1 (MQ=255)
                                                                                                    tCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCTTGAACTGGggcgg                                              <  2:2353135/150‑1 (MQ=255)
                                                                                                     cAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGggcgg                                              >  1:379969/1‑150 (MQ=255)
                                                                                                       aCCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCt                                            >  1:2529792/1‑150 (MQ=255)
                                                                                                        ccTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTg                                           >  1:2225315/1‑150 (MQ=255)
                                                                                                          tataCCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTACTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGc                                         <  2:2200421/150‑1 (MQ=255)
                                                                                                            taCCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGa                                       <  2:4132466/150‑1 (MQ=255)
                                                                                                              ccATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGNACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAAt                                     <  1:2677336/150‑1 (MQ=255)
                                                                                                              ccATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAAt                                     <  2:1764490/150‑1 (MQ=255)
                                                                                                               cATTAAATAACAGCAACTTATACCTGTCTGACGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAAtg                                    >  1:4120406/1‑150 (MQ=255)
                                                                                                               cATTAAATAACAGCAACTTAGACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAAtg                                    <  2:4275523/150‑1 (MQ=255)
                                                                                                                  tAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGa                                 >  2:1844802/1‑150 (MQ=255)
                                                                                                                    aaTAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAAt                               >  1:2150723/1‑150 (MQ=255)
                                                                                                                    aaTAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAAt                               >  1:1286264/1‑150 (MQ=255)
                                                                                                                       aaCAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGACGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCg                            >  2:2166604/1‑150 (MQ=255)
                                                                                                                         cAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGcc                          >  1:196226/1‑150 (MQ=255)
                                                                                                                          aGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCAGCCt                         >  1:1867127/1‑150 (MQ=255)
                                                                                                                              aCTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAat                     <  2:3319749/150‑1 (MQ=255)
                                                                                                                                ttATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAatta                   >  1:4003177/1‑150 (MQ=255)
                                                                                                                                 tataCCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAattat                  >  1:1641059/1‑150 (MQ=255)
                                                                                                                                  ataCCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAattatt                 >  2:3594681/1‑150 (MQ=255)
                                                                                                                                  ataCCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAattatt                 >  1:3643094/1‑150 (MQ=255)
                                                                                                                                  ataCCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAattatt                 >  1:1493485/1‑150 (MQ=255)
                                                                                                                                  ataCCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAattatt                 >  2:4423991/1‑150 (MQ=255)
                                                                                                                                   taCCTGTCTGATGTCGCTCTCTGCGAAGCGGTAACTCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCTTATGTGAATCCGCATGAATTATTg                <  2:5023889/150‑1 (MQ=255)
                                                                                                                                   taCCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTg                >  1:4705142/1‑150 (MQ=255)
                                                                                                                                     ccTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTgcg              >  1:1530922/1‑150 (MQ=255)
                                                                                                                                     ccTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCGTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTgcg              >  1:5066621/1‑150 (MQ=255)
                                                                                                                                       tGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCt            <  1:258943/150‑1 (MQ=255)
                                                                                                                                       tGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCt            <  1:1557435/150‑1 (MQ=255)
                                                                                                                                        gTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTa           >  1:767284/1‑150 (MQ=255)
                                                                                                                                        ctCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGACAGCAGTTAGGAACGGCTGAATCCCTTGAACTCGGGCGGCTGGCGAATGTGAATCCACCTGAATGATTGCGCTa           >  2:2971031/2‑150 (MQ=255)
                                                                                                                                         tCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTAc          >  2:3958023/1‑150 (MQ=255)
                                                                                                                                         tCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTAc          >  1:3009920/1‑150 (MQ=255)
                                                                                                                                                 cgctcTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGATGCGCa  >  2:4472911/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGATGCGCA  >  CP009273/4404021‑4404315

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: