Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,404,169 G→T 90.3% 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,1690GT90.3% 569.4 / 56.2 217A26A (GCG→GCT) aidBDNA alkylation damage repair protein; flavin‑containing DNA binding protein, weak isovaleryl CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (11/10);  new base T (116/79);  total (128/89)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 223 total aligned reads displayed.

TCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGATGC  >  CP009273/4404023‑4404312
                                                                                                                                                  |                                                                                                                                               
 ctctTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGctct                                                                                                                                               >  1:100113/1‑148 (MQ=255)
 ctctTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGctct                                                                                                                                               <  1:819345/148‑1 (MQ=255)
    tttCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTg                                                                                                                                            <  1:297346/148‑1 (MQ=255)
         gTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAg                                                                                                                                       >  1:514777/1‑148 (MQ=255)
            aTGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCgg                                                                                                                                    <  1:1860685/148‑1 (MQ=255)
              gAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTa                                                                                                                                  <  1:1318219/148‑1 (MQ=255)
                 tCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAAcg                                                                                                                               <  2:815931/148‑1 (MQ=255)
                 tCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAAcg                                                                                                                               <  1:3188182/148‑1 (MQ=255)
                        cAGTGACCATACTACTGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAg                                                                                                                        <  1:1332587/148‑1 (MQ=255)
                        cAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAg                                                                                                                        <  2:1046162/148‑1 (MQ=255)
                           tGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTg                                                                                                                     <  1:1691702/148‑1 (MQ=255)
                           tGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTCACGCGTGAAGGTg                                                                                                                     >  1:3137406/1‑148 (MQ=255)
                            gACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTgc                                                                                                                    <  2:2960380/148‑1 (MQ=255)
                            gACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTgc                                                                                                                    <  1:2781949/148‑1 (MQ=255)
                                aTACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCACCTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTgctggc                                                                                                                >  2:2283190/1‑148 (MQ=255)
                                 tACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTCACGCGTGAAGGTgctggct                                                                                                               >  2:2076053/1‑148 (MQ=255)
                                 tACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTgctggct                                                                                                               <  2:1446773/148‑1 (MQ=255)
                                  aCTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTgctggctg                                                                                                              <  1:435723/148‑1 (MQ=255)
                                   cTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGctggctgg                                                                                                             <  2:2619292/148‑1 (MQ=255)
                                    tAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTggg                                                                                                            <  2:153587/148‑1 (MQ=255)
                                      aTGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGAt                                                                                                          <  1:1105728/148‑1 (MQ=255)
                                       tGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATa                                                                                                         <  2:701767/148‑1 (MQ=255)
                                        ggTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAg                                                                                                        >  2:1024819/1‑148 (MQ=255)
                                        ggTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAg                                                                                                        <  1:2740321/148‑1 (MQ=255)
                                             cTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGAtt                                                                                                   <  2:2681457/148‑1 (MQ=255)
                                              tGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGAttt                                                                                                  <  2:1813695/148‑1 (MQ=255)
                                                 cATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCt                                                                                               <  1:1736654/148‑1 (MQ=255)
                                                  aTTGATGGAGGTAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACATATAACATTAAATAACAGCAACTTATACCTCTCTTATGGAGCTCTCTGCGAATCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATATTCTa                                                                                              >  2:3277417/1‑148 (MQ=255)
                                                  aTTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTa                                                                                              >  2:1217142/1‑148 (MQ=255)
                                                  aTTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTa                                                                                              >  2:1640024/1‑148 (MQ=255)
                                                          agggagACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTAt                                                                                      <  2:725381/148‑1 (MQ=255)
                                                           gggagACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTCtt                                                                                     >  2:691654/1‑148 (MQ=255)
                                                            ggagACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTg                                                                                    <  2:2740184/148‑1 (MQ=255)
                                                               gaCACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTc                                                                                 <  1:570522/148‑1 (MQ=255)
                                                               gaCACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTc                                                                                 >  2:2473685/1‑148 (MQ=255)
                                                                 cacaGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTcag                                                                               <  1:2650907/148‑1 (MQ=255)
                                                                   caGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTcagca                                                                             >  2:904235/1‑148 (MQ=255)
                                                                   caGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTcagca                                                                             >  2:755932/1‑148 (MQ=255)
                                                                    aGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTcagcag                                                                            >  1:1373326/1‑148 (MQ=255)
                                                                     gTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGt                                                                           >  1:480762/1‑148 (MQ=255)
                                                                     gTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGt                                                                           <  1:757166/148‑1 (MQ=255)
                                                                      tGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGtt                                                                          >  2:2817845/1‑148 (MQ=255)
                                                                             gCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGAGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACg                                                                   >  1:1488315/1‑148 (MQ=255)
                                                                                  cTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGa                                                              <  1:691368/148‑1 (MQ=255)
                                                                                  cTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGa                                                              >  1:787908/1‑148 (MQ=255)
                                                                                   tCACACCGTTTTTAATCGACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTCACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGaa                                                             >  2:1697214/1‑148 (MQ=255)
                                                                                      cacCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATcc                                                          >  2:534453/1‑148 (MQ=255)
                                                                                       acCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATccc                                                         <  2:2697746/148‑1 (MQ=255)
                                                                                         cGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCtt                                                       >  2:743792/1‑148 (MQ=255)
                                                                                         cGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCtt                                                       >  1:2358406/1‑148 (MQ=255)
                                                                                          gTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTCACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTg                                                      >  1:1340456/1‑148 (MQ=255)
                                                                                           tttttAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGa                                                     >  2:2804829/1‑148 (MQ=255)
                                                                                           tttttAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGa                                                     <  2:183162/148‑1 (MQ=255)
                                                                                           tttttAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGa                                                     >  2:1144071/1‑148 (MQ=255)
                                                                                            ttttAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGaa                                                    >  2:3269795/1‑148 (MQ=255)
                                                                                            ttttAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGaa                                                    >  2:3177190/1‑148 (MQ=255)
                                                                                             tttAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAAc                                                   <  2:2872869/148‑1 (MQ=255)
                                                                                                   cAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGggc                                             >  1:2209856/1‑148 (MQ=255)
                                                                                                   cAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGggc                                             >  2:380813/1‑148 (MQ=255)
                                                                                                   cAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGggc                                             >  2:2904940/1‑148 (MQ=255)
                                                                                                        tataCCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTg                                        >  2:1773863/1‑148 (MQ=255)
                                                                                                        tataCCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTg                                        >  1:1530643/1‑148 (MQ=255)
                                                                                                         ataCCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTgg                                       >  2:1095187/1‑148 (MQ=255)
                                                                                                         ataCCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTgg                                       >  2:1978023/1‑148 (MQ=255)
                                                                                                         ataCCATTAAATAACAGCAACTTATACATGTCTGATGGAGCACTCTGCTAAGCGGTAACGCGTGAAGGTGCTGGCTGGTATAGCGATTTTATAGCCAGTATTGGTAAGAAGTTAGGAACGGATGAATCCCTTGAACTGGGGCGGATgg                                       >  1:941352/1‑148 (MQ=255)
                                                                                                           aCCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCg                                     >  1:2675840/1‑148 (MQ=255)
                                                                                                             cATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGaa                                   >  1:347393/1‑148 (MQ=255)
                                                                                                              aTTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAAt                                  >  2:455739/1‑148 (MQ=255)
                                                                                                              aTTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAAt                                  >  2:314549/1‑148 (MQ=255)
                                                                                                              aTTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAAt                                  >  1:554390/1‑148 (MQ=255)
                                                                                                              aTTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAAt                                  <  1:2881804/148‑1 (MQ=255)
                                                                                                               ttAAATAACAGCAACTTCTACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAAtg                                 <  2:1921393/148‑1 (MQ=255)
                                                                                                                tAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAAtgt                                <  1:3009377/148‑1 (MQ=255)
                                                                                                                tAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAAtgt                                <  2:278695/148‑1 (MQ=255)
                                                                                                                  aaTAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGa                              >  1:162177/1‑148 (MQ=255)
                                                                                                                  aaTAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGa                              <  2:727224/148‑1 (MQ=255)
                                                                                                                   aTAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGaa                             >  1:2738143/1‑148 (MQ=255)
                                                                                                                   aTAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGaa                             <  2:3027836/148‑1 (MQ=255)
                                                                                                                   aTAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGaa                             <  2:2073061/148‑1 (MQ=255)
                                                                                                                   aTAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGaa                             >  1:2931702/1‑148 (MQ=255)
                                                                                                                   aTAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGaa                             <  2:2212662/148‑1 (MQ=255)
                                                                                                                    tAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAAt                            >  1:3269250/1‑148 (MQ=255)
                                                                                                                    tAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAAt                            <  2:235758/148‑1 (MQ=255)
                                                                                                                    tAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAAt                            <  2:2205602/148‑1 (MQ=255)
                                                                                                                         gCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGcc                       >  1:990120/1‑148 (MQ=255)
                                                                                                                         gCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGcc                       >  2:341881/1‑148 (MQ=255)
                                                                                                                         gCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGcc                       >  1:1748439/1‑148 (MQ=255)
                                                                                                                               tataCCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAatt                 >  1:1519651/1‑148 (MQ=255)
                                                                                                                                  aCCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAattatt              >  1:293897/1‑148 (MQ=255)
                                                                                                                                  aCCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAattatt              >  1:879184/1‑148 (MQ=255)
                                                                                                                                   ccTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGa                                                     >  1:1442599/1‑108 (MQ=255)
                                                                                                                                    cTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTgc            >  2:3021064/1‑148 (MQ=255)
                                                                                                                                     tGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTgcg           >  2:564885/1‑148 (MQ=255)
                                                                                                                                        cTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTa        >  2:2002039/1‑148 (MQ=255)
                                                                                                                                        cTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTa        <  2:2146180/148‑1 (MQ=255)
                                                                                                                                           aTGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGa     <  1:2167572/148‑1 (MQ=255)
                                                                                                                                           aTGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGa     >  1:1030607/1‑148 (MQ=255)
                                                                                                                                             ggCGCTCTCTGCGAAGCGGTCACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGATg   >  1:208498/1‑148 (MQ=255)
                                                                                                                                             ggCGCTCTCTGCGAAGCGGTCACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGATg   >  1:1240917/1‑148 (MQ=255)
                                                                                                                                             ggCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGATg   <  1:2398109/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
TCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGATGC  >  CP009273/4404023‑4404312

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: ‑