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,1690GT98.2% 701.5 / 1.9 219A26A (GCG→GCTaidBDNA alkylation damage repair protein; flavin‑containing DNA binding protein, weak isovaleryl CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (3/1);  new base T (118/97);  total (121/98)
Fisher's exact test for biased strand distribution p-value = 6.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 228 total aligned reads displayed.

CTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGATGCGCA  >  CP009273/4404024‑4404315
                                                                                                                                                 |                                                                                                                                                  
ctctTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGctct                                                                                                                                                  >  2:2778996/1‑148 (MQ=255)
ctctTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCAATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGctct                                                                                                                                                  <  1:1754237/148‑1 (MQ=255)
  ctTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGctctct                                                                                                                                                >  1:3643690/1‑148 (MQ=255)
    ttCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGc                                                                                                                                              >  1:2803626/1‑148 (MQ=255)
        gTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAg                                                                                                                                          >  1:3323206/1‑148 (MQ=255)
             gAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTa                                                                                                                                     >  2:640519/1‑148 (MQ=255)
               aTCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAAc                                                                                                                                   >  1:2463855/1‑148 (MQ=255)
                tCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAAcg                                                                                                                                  <  2:2645815/148‑1 (MQ=255)
                   aTGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGt                                                                                                                               <  2:199471/148‑1 (MQ=255)
                          tGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTg                                                                                                                        >  1:2784661/1‑148 (MQ=255)
                             ccATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTgctg                                                                                                                     >  2:2299390/1‑148 (MQ=255)
                                 acgaaTGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACGTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTgctggctg                                                                                                                 <  2:1792345/145‑1 (MQ=255)
                                 aCTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTgctggctg                                                                                                                 <  2:179938/148‑1 (MQ=255)
                                   tAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTggg                                                                                                               >  1:2615495/1‑148 (MQ=255)
                                     aTGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGAt                                                                                                             <  1:1206249/148‑1 (MQ=255)
                                      tGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATa                                                                                                            >  2:3285950/1‑148 (MQ=255)
                                      tGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATa                                                                                                            >  2:1631386/1‑148 (MQ=255)
                                       ggTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAg                                                                                                           >  2:3556853/1‑148 (MQ=255)
                                       ggTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAg                                                                                                           <  2:1694216/148‑1 (MQ=255)
                                          gACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGa                                                                                                        <  1:1518031/148‑1 (MQ=255)
                                            cTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGAtt                                                                                                      <  1:3505970/148‑1 (MQ=255)
                                              gCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGCTAGCGAtttt                                                                                                    >  2:631086/1‑148 (MQ=255)
                                              gCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGAtttt                                                                                                    <  1:440212/148‑1 (MQ=255)
                                              gCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGAtttt                                                                                                    >  2:2552159/1‑148 (MQ=255)
                                              gCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGAtttt                                                                                                    <  2:2941472/148‑1 (MQ=255)
                                              gCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGAtttt                                                                                                    >  2:2319009/1‑148 (MQ=255)
                                                cATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCt                                                                                                  >  2:576501/1‑148 (MQ=255)
                                                   tGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGc                                                                                               >  2:387716/1‑148 (MQ=255)
                                                     aTGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCa                                                                                             >  2:3645497/1‑148 (MQ=255)
                                                         agggagACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTAt                                                                                         <  1:1846937/148‑1 (MQ=255)
                                                         agggagACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTAt                                                                                         <  2:943468/148‑1 (MQ=255)
                                                            gagaCACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTgg                                                                                      <  1:2923913/148‑1 (MQ=255)
                                                            gagaCACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTgg                                                                                      <  1:3309594/148‑1 (MQ=255)
                                                            gagaCACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTgg                                                                                      <  1:1400442/148‑1 (MQ=255)
                                                             agaCACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGt                                                                                     <  1:3103214/148‑1 (MQ=255)
                                                                cacaGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTcag                                                                                  >  2:884471/1‑148 (MQ=255)
                                                                 acaGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTcagc                                                                                 <  2:2012561/148‑1 (MQ=255)
                                                                     tGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGtt                                                                             >  1:3020700/1‑148 (MQ=255)
                                                                         cTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGa                                                                         <  1:1123136/148‑1 (MQ=255)
                                                                             cAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACgg                                                                     <  2:1865620/148‑1 (MQ=255)
                                                                                aCTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTg                                                                  >  1:54989/1‑148 (MQ=255)
                                                                                aCTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTg                                                                  >  2:2715674/1‑148 (MQ=255)
                                                                                aCTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTg                                                                  >  1:970354/1‑148 (MQ=255)
                                                                                 cTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGa                                                                 <  1:2542608/148‑1 (MQ=255)
                                                                                 cTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGa                                                                 >  1:615900/1‑148 (MQ=255)
                                                                                    acacCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATc                                                              <  2:3391358/148‑1 (MQ=255)
                                                                                    acacCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATc                                                              <  2:160180/148‑1 (MQ=255)
                                                                                    acacCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATc                                                              <  2:1485480/148‑1 (MQ=255)
                                                                                     cacCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATcc                                                             >  1:3603182/1‑148 (MQ=255)
                                                                                      acCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATccc                                                            <  2:540762/148‑1 (MQ=255)
                                                                                       ccGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCt                                                           >  2:2730758/1‑148 (MQ=255)
                                                                                       ccGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCt                                                           >  2:78616/1‑148 (MQ=255)
                                                                                       ccGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCt                                                           >  2:78894/1‑148 (MQ=255)
                                                                                        cGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCtt                                                          >  2:1752050/1‑148 (MQ=255)
                                                                                           ttttAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGaa                                                       <  1:906339/148‑1 (MQ=255)
                                                                                           ttttAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGaa                                                       <  1:2513378/148‑1 (MQ=255)
                                                                                                aTCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTggg                                                  >  1:225381/1‑148 (MQ=255)
                                                                                                aTCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTggg                                                  >  2:3170934/1‑148 (MQ=255)
                                                                                                aTCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTggg                                                  >  2:87123/1‑148 (MQ=255)
                                                                                                     ccTATACCATTAAATAACAGCAACATATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCACTAGAACTGGggcggc                                             >  1:1642085/1‑148 (MQ=255)
                                                                                                      cTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCt                                            <  1:1586871/148‑1 (MQ=255)
                                                                                                       tataCCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTg                                           <  2:1692250/148‑1 (MQ=255)
                                                                                                        ataCCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTgg                                          >  1:548377/1‑148 (MQ=255)
                                                                                                        ataCCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTgg                                          <  2:1091830/148‑1 (MQ=255)
                                                                                                           ccATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGa                                       <  2:3275075/148‑1 (MQ=255)
                                                                                                            cATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGaa                                      <  1:2291701/148‑1 (MQ=255)
                                                                                                            cATTAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGaa                                      <  1:2688233/148‑1 (MQ=255)
                                                                                                              ttAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAAtg                                    >  2:3702638/1‑148 (MQ=255)
                                                                                                               tAAATAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAAtgt                                   >  1:3742272/1‑148 (MQ=255)
                                                                                                                 aaTAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGa                                 <  2:2047728/148‑1 (MQ=255)
                                                                                                                  aTAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGaa                                <  1:3270243/148‑1 (MQ=255)
                                                                                                                   tAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTCACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAAt                               >  1:438813/1‑148 (MQ=255)
                                                                                                                   tAACAGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAAt                               >  2:2628575/1‑148 (MQ=255)
                                                                                                                       aGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGc                           >  1:141724/1‑148 (MQ=255)
                                                                                                                       aGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGc                           <  2:696559/148‑1 (MQ=255)
                                                                                                                       aGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGc                           <  1:1008381/148‑1 (MQ=255)
                                                                                                                       aGCAACTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGc                           <  2:2615495/148‑1 (MQ=255)
                                                                                                                           aCTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGa                       >  2:3181484/1‑148 (MQ=255)
                                                                                                                           aCTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGa                       <  2:3321515/148‑1 (MQ=255)
                                                                                                                           aCTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGa                       <  2:3323206/148‑1 (MQ=255)
                                                                                                                           aCTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGa                       <  1:2900571/148‑1 (MQ=255)
                                                                                                                           aCTTATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGa                       >  1:1615941/1‑148 (MQ=255)
                                                                                                                           aCTTATACCTGACTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTCGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTTGGGCGGCTGGCGAATGTGAATCCGCCTGa                       >  2:1231214/1‑148 (MQ=255)
                                                                                                                             ttATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAat                     >  1:301030/1‑148 (MQ=255)
                                                                                                                             ttATACCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAat                     >  1:865177/1‑148 (MQ=255)
                                                                                                                                taCCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAattat                  <  1:1695494/148‑1 (MQ=255)
                                                                                                                                taCCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAattat                  <  2:2426195/148‑1 (MQ=255)
                                                                                                                                 aCCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAattatt                 >  1:3361129/1‑148 (MQ=255)
                                                                                                                                 aCCTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAattatt                 >  1:1311835/1‑148 (MQ=255)
                                                                                                                                  ccTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTg                >  2:114022/1‑148 (MQ=255)
                                                                                                                                  ccTGTCTGATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTg                <  1:852696/148‑1 (MQ=255)
                                                                                                                                         gATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACg         <  2:2905573/148‑1 (MQ=255)
                                                                                                                                         gATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACg         >  2:626519/1‑148 (MQ=255)
                                                                                                                                         gATGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACg         >  1:752159/1‑148 (MQ=255)
                                                                                                                                          aTGGCGCTCTCTGCGAAGCGGTAACGCGTGAATGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGa        <  1:2870776/148‑1 (MQ=255)
                                                                                                                                          aTGGCGCTCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGa        <  1:3792375/148‑1 (MQ=255)
                                                                                                                                           tGGCGCTCTCTGCGAAGCGGTCACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGAt       >  2:1494083/1‑148 (MQ=255)
                                                                                                                                              cgctcTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGATgcg    >  2:809011/1‑148 (MQ=255)
                                                                                                                                              cgctcTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGATgcg    >  1:590968/1‑148 (MQ=255)
                                                                                                                                                cgcTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGATGCGCa  >  1:2985955/1‑148 (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: