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,1690GT96.9% 600.7 / 9.1 195A26A (GCG→GCTaidBDNA alkylation damage repair protein; flavin‑containing DNA binding protein, weak isovaleryl CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (3/3);  new base T (100/89);  total (103/92)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 202 total aligned reads displayed.

CTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGATGCGCA  >  CP009273/4404024‑4404315
                                                                                                                                                 |                                                                                                                                                  
ctctTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGctct                                                                                                                                                  <  1:1569060/148‑1 (MQ=255)
 tctTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCtctc                                                                                                                                                 >  1:1034531/1‑148 (MQ=255)
      cTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGa                                                                                                                                            >  2:709518/1‑148 (MQ=255)
       tGTCATGAATCTATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGaa                                                                                                                                           >  2:624529/1‑148 (MQ=255)
         tCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGc                                                                                                                                         >  1:2001701/1‑148 (MQ=255)
             gAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTa                                                                                                                                     <  2:1167225/148‑1 (MQ=255)
               aTCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAAc                                                                                                                                   >  2:1957849/1‑148 (MQ=255)
                     ggCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTCACGCGTGa                                                                                                                             >  1:2889064/1‑148 (MQ=255)
                     ggCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGa                                                                                                                             >  1:976015/1‑148 (MQ=255)
                       cAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAg                                                                                                                           <  1:1489982/148‑1 (MQ=255)
                        aGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAgg                                                                                                                          <  2:105084/148‑1 (MQ=255)
                         gTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGt                                                                                                                         <  2:2514390/148‑1 (MQ=255)
                         gTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGt                                                                                                                         >  2:1744615/1‑148 (MQ=255)
                           gACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTgc                                                                                                                       <  2:707241/148‑1 (MQ=255)
                            aCCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTgct                                                                                                                      <  1:544206/148‑1 (MQ=255)
                              cATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTgctgg                                                                                                                    >  1:1278024/1‑148 (MQ=255)
                              cATAATAATGTTGAATGCCATTGATGGATGGAGACACCGTTCACTGGCAAACTCACACAGTTTTTAATCAACCTATAACATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTgcttg                                                                                                                    <  1:404980/148‑3 (MQ=255)
                               aTACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTgctggc                                                                                                                   >  2:1217749/1‑148 (MQ=255)
                                 aCTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTgctggctg                                                                                                                 <  1:463042/148‑1 (MQ=255)
                                 aCTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTgctggctg                                                                                                                 >  1:2191032/1‑148 (MQ=255)
                                 aCTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTgctggctg                                                                                                                 >  2:1333481/1‑148 (MQ=255)
                                    aaTGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGa                                                                                                              <  2:3545226/148‑1 (MQ=255)
                                    aaTGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGa                                                                                                              <  2:976015/148‑1 (MQ=255)
                                     aTGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGAt                                                                                                             >  1:393943/1‑148 (MQ=255)
                                       ggTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTCACGCGTGAAGGTGTTGGCTGGGATAg                                                                                                           >  1:431382/1‑148 (MQ=255)
                                       ggTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAg                                                                                                           <  2:1998953/148‑1 (MQ=255)
                                         tGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCg                                                                                                         <  1:1594968/148‑1 (MQ=255)
                                         tGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCg                                                                                                         <  1:3361960/148‑1 (MQ=255)
                                           aCTGCCATTGATGGAGGGAGACACAGTGCACTGACAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGAt                                                                                                       >  2:2756262/1‑148 (MQ=255)
                                            cTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGAtt                                                                                                      >  2:2327905/1‑148 (MQ=255)
                                             tGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGAttt                                                                                                     <  2:1306360/148‑1 (MQ=255)
                                                cATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCt                                                                                                  >  1:370233/1‑148 (MQ=255)
                                                cATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCt                                                                                                  <  2:173844/148‑1 (MQ=255)
                                                cATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCt                                                                                                  <  2:1856509/148‑1 (MQ=255)
                                                   tGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGc                                                                                               >  2:1858879/1‑148 (MQ=255)
                                                   tGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGc                                                                                               >  2:835492/1‑148 (MQ=255)
                                                   tGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGc                                                                                               >  1:2854647/1‑148 (MQ=255)
                                                   tGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGc                                                                                               >  2:53293/1‑148 (MQ=255)
                                                     aTGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCa                                                                                             <  2:780587/148‑1 (MQ=255)
                                                        gagggagACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTa                                                                                          <  1:2697721/148‑1 (MQ=255)
                                                         agggagACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTAt                                                                                         <  1:605812/148‑1 (MQ=255)
                                                                  caGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTcagca                                                                                >  1:1386627/1‑148 (MQ=255)
                                                                   aGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTcagcag                                                                               >  2:640506/1‑148 (MQ=255)
                                                                    gTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGt                                                                              <  1:2723838/148‑1 (MQ=255)
                                                                     tGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGtt                                                                             >  1:63935/1‑148 (MQ=255)
                                                                     tGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGtt                                                                             <  2:2369748/148‑1 (MQ=255)
                                                                      gCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTa                                                                            <  1:1654692/148‑1 (MQ=255)
                                                                       cACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAg                                                                           >  1:1372472/1‑148 (MQ=255)
                                                                       cACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAg                                                                           <  1:1097777/148‑1 (MQ=255)
                                                                             cAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACgg                                                                     >  1:3201473/1‑148 (MQ=255)
                                                                              aaaCTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGc                                                                    >  1:2207330/1‑148 (MQ=255)
                                                                              aaaCTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGc                                                                    >  2:2293038/1‑148 (MQ=255)
                                                                                aCTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTg                                                                  >  1:2650644/1‑148 (MQ=255)
                                                                                   cacaCCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAAt                                                               <  1:893087/148‑1 (MQ=255)
                                                                                     cacCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATcc                                                             >  2:3414709/1‑148 (MQ=255)
                                                                                          tttttAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGa                                                        <  2:2797839/148‑1 (MQ=255)
                                                                                          tttttAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGa                                                        >  2:1462449/1‑148 (MQ=255)
                                                                                           ttttAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGaa                                                       <  1:607805/148‑1 (MQ=255)
                                                                                           ttttAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGaa                                                       >  1:2309993/1‑148 (MQ=255)
                                                                                           ttttAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGaa                                                       >  1:2800632/1‑148 (MQ=255)
                                                                                             ttAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACt                                                     <  2:899320/148‑1 (MQ=255)
                                                                                              tAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTg                                                    <  2:3321640/148‑1 (MQ=255)
                                                                                                       tataCCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTg                                           >  2:2877689/1‑148 (MQ=255)
                                                                                                       tataCCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTg                                           >  2:2877592/1‑148 (MQ=255)
                                                                                                       tataCCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTg                                           >  1:3105522/1‑148 (MQ=255)
                                                                                                       tataCCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTg                                           >  1:3107681/1‑148 (MQ=255)
                                                                                                         taCCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGc                                         >  2:2636812/1‑148 (MQ=255)
                                                                                                           ccATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGa                                       <  1:627349/148‑1 (MQ=255)
                                                                                                            cATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGaa                                      >  1:1777839/1‑148 (MQ=255)
                                                                                                            cATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGaa                                      <  2:932148/148‑1 (MQ=255)
                                                                                                             aTTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAAt                                     >  1:2240108/1‑148 (MQ=255)
                                                                                                             aTTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGAGAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGact                                     >  1:2642600/1‑146 (MQ=255)
                                                                                                             aTTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAAt                                     >  2:1234094/1‑148 (MQ=255)
                                                                                                              ttAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAAtg                                    <  1:813513/148‑1 (MQ=255)
                                                                                                                 aaTAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGa                                 >  1:663311/1‑148 (MQ=255)
                                                                                                                  aTAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGaa                                >  1:1735253/1‑148 (MQ=255)
                                                                                                                    aaCAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATc                              >  2:2866302/1‑148 (MQ=255)
                                                                                                                     aCAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATcc                             >  2:2223144/1‑148 (MQ=255)
                                                                                                                     aCAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATcc                             >  2:700211/1‑148 (MQ=255)
                                                                                                                     aCAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATcc                             >  1:2845705/1‑148 (MQ=255)
                                                                                                                            cTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGaa                      <  2:1569648/148‑1 (MQ=255)
                                                                                                                            cTTATACCTGTCTGATGGCACTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGaa                      >  1:2105185/1‑148 (MQ=255)
                                                                                                                              tataCCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAatt                    <  2:912735/148‑1 (MQ=255)
                                                                                                                              tataCCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAatt                    <  1:3001950/148‑1 (MQ=255)
                                                                                                                              tataCCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAatt                    >  2:2429586/1‑148 (MQ=255)
                                                                                                                               ataCCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAatta                   >  2:1643839/1‑148 (MQ=255)
                                                                                                                                 aCCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAattatt                 >  1:448754/1‑148 (MQ=255)
                                                                                                                                  ccTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTg                >  1:3274015/1‑148 (MQ=255)
                                                                                                                                  ccTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTg                <  1:1359069/148‑1 (MQ=255)
                                                                                                                                  ccTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTg                >  1:2072007/1‑148 (MQ=255)
                                                                                                                                   cTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTgc               >  1:1228041/1‑148 (MQ=255)
                                                                                                                                    tGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTgcg              >  2:3405603/1‑148 (MQ=255)
                                                                                                                                    tGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTgcg              >  2:1651527/1‑148 (MQ=255)
                                                                                                                                       cTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTa           <  2:370233/148‑1 (MQ=255)
                                                                                                                                       cTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTa           <  1:1583541/148‑1 (MQ=255)
                                                                                                                                           tGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGACTTATTGCGCTACGAt       >  1:2495663/1‑148 (MQ=255)
                                                                                                                                            ggCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGATg      >  2:1192378/1‑148 (MQ=255)
                                                                                                                                             gcgctcTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGATgc     >  1:2955462/1‑148 (MQ=255)
                                                                                                                                             gcgctcTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGATgc     >  1:1176991/1‑148 (MQ=255)
                                                                                                                                                cgcTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGATGCGCa  <  1:1077873/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGATGCGCA  >  CP009273/4404024‑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: