Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,182,646 G→T 7.9% I40I (ATC→ATA thiS ← immediate sulfur donor in thiazole formation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,182,6460GT 7.9% 1095.3 / 39.8 296I40I (ATC→ATAthiSimmediate sulfur donor in thiazole formation
Reads supporting (aligned to +/- strand):  ref base G (149/123);  new base T (16/8);  total (165/131)
Fisher's exact test for biased strand distribution p-value = 2.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 298 total aligned reads displayed.

AGAAGCGAATTTCCCTGTGCCGGTAAACAGATGTGAATCAAACGTTTTGTCCGCAATACGTAACATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGTTAAACAGGATCTGCATTGCTTCCTCCGCATACCGGGCAACCACT  >  CP009273/4182498‑4182793
                                                                                                                                                    |                                                                                                                                                   
agaagCGAATTTCCCTGTGCCGGTAAACAGATGTGAATCAAACGTTTTGTCCGCAATACGTAACATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACg                                                                                                                                                     <  1:4264064/149‑1 (MQ=255)
   agCGAATTTCCCTGTGCCGGTAAACAGATGTGAATCAAACGTTTTGTCCGCAATACGTAACATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGAtt                                                                                                                                                  <  1:2113594/149‑1 (MQ=255)
        aTTTCCCTGTGCCGGTAAACAGATGTGAATCAAACGTTTTGTCCGCAATACGTAACATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGAGTTGCTg                                                                                                                                             >  1:3203979/1‑149 (MQ=255)
         tttCCCTGTGCCGGTAAACAGATGTGAATCAAACGTTTTGTCCGCAATACGTAACATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGa                                                                                                                                            >  1:497655/1‑149 (MQ=255)
         tttCCCTGTGCCGGTAAACAGATGTGAATCAAACGTTTTGTCCGCAATACGTAACATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGa                                                                                                                                            >  1:2235933/1‑149 (MQ=255)
           tCCCTGTGCCGGTAAACAGATGTGAATCAAACGTTTTGTCCGCAATACGTAACATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGAtt                                                                                                                                          <  1:2592195/149‑1 (MQ=255)
               tgtgCCGGTAAACAGATGTGAATCAAACGTTTTGTCCGCAATACGTAACATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATc                                                                                                                                      <  1:608879/149‑1 (MQ=255)
                   ccGGTAAACAGATGTGAATCAAACGTTTTGTCCGCAATACGTAACATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCa                                                                                                                                  <  1:176553/149‑1 (MQ=255)
                    cGGTAAACAGATGTGAATCAAACGTTTTGTCCGCAATACGTAACATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCag                                                                                                                                 <  2:88952/149‑1 (MQ=255)
                      gTAAACAGATGTGAATCAAACGTTTTGTCCGCAATACGTAACATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCagag                                                                                                                               >  1:2818076/1‑149 (MQ=255)
                       tAAACAGATGTGAATCAAACGTTTTGTCCGCAATACGTAACATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGc                                                                                                                              >  2:4023091/1‑149 (MQ=255)
                            aGATGTGAATCAAACGTTTTGTCCGCAATACGTAACATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCgcgc                                                                                                                         >  1:2694839/1‑149 (MQ=255)
                                gtgAATCAAACGTTTTTTCCGCACTACGTAACATGTCAACCCCCTGCAATATCCTGAAAAAGCAGGATCTGGTCGCCATCCTGCAAGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTTATTAATCGCCAGAGCCGCGCCCGc                                                                                                                     <  1:541685/149‑1 (MQ=255)
                                  gAATCAAACGTTTTGTCCGCAATACGTAACATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACTATTTGCTGATTAATCGCCAGAGCCGCGCCCGCtt                                                                                                                   <  1:4402674/149‑1 (MQ=255)
                                  gAATCAAACGTTTTGTCCGCAATACGTAACATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCtt                                                                                                                   >  2:3416910/1‑149 (MQ=255)
                                   aaTCAAACGTTTTGTCCGCAATACGTAACATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACTATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTg                                                                                                                  >  2:509033/1‑149 (MQ=255)
                                    aTCAAACGTTTTGTCCGCAATACGTAACATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACTATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGt                                                                                                                 <  2:4185787/149‑1 (MQ=255)
                                    aTCAAACGTTTTGTCCGCAATACGTAACATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGt                                                                                                                 >  2:737213/1‑149 (MQ=255)
                                    aTCAAACGTTTTGTCCGCAATACGTAACATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGt                                                                                                                 >  2:3343172/1‑149 (MQ=255)
                                    aTCAAACGTTTTGTCCGCAATACGTAACATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGAGTAATCGCCAGAGCCGCGCCCGCTTGt                                                                                                                 >  2:1844575/1‑149 (MQ=255)
                                     tCAAACGTTTTGTCCGCAATACGTAACATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTc                                                                                                                <  2:690436/149‑1 (MQ=255)
                                        aaCGTTTTGTCCGCAATACGTAACATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGtt                                                                                                             >  2:3895788/1‑149 (MQ=255)
                                           gTTTTGTCCGCAATACGTAACATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGt                                                                                                          >  1:2866172/1‑149 (MQ=255)
                                                 tCCGCAATACGTAACATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACTATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGtt                                                                                                    >  2:4730207/1‑149 (MQ=255)
                                                      aaTACGTAACATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTcc                                                                                               >  2:454148/1‑149 (MQ=255)
                                                       aTACGTAACATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACTATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTAGGTCGAGTTGCTCCa                                                                                              >  2:4192942/1‑149 (MQ=255)
                                                       aTACGTAACATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCa                                                                                              >  2:2351031/1‑149 (MQ=255)
                                                        tACGTAACATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCag                                                                                             <  2:1477367/149‑1 (MQ=255)
                                                        tACGTAACATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCag                                                                                             >  2:3562518/1‑149 (MQ=255)
                                                        tACGTAACATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCag                                                                                             >  2:630025/1‑149 (MQ=255)
                                                           gTAACATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGGTGGTCGAGTTGCTCCagtag                                                                                          <  1:3319841/149‑1 (MQ=255)
                                                            tAACATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCagtagt                                                                                         <  1:3587238/149‑1 (MQ=255)
                                                             aaCATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGtt                                                                                        <  2:1414442/149‑1 (MQ=255)
                                                              aCATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTc                                                                                       <  1:1483910/149‑1 (MQ=255)
                                                              aCATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTc                                                                                       >  1:684706/1‑149 (MQ=255)
                                                              aCATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTAGTCGTTGGTCGAGTTGCTCCAGTAGTTc                                                                                       >  2:1663795/1‑149 (MQ=255)
                                                               cATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCg                                                                                      >  1:626071/1‑149 (MQ=255)
                                                                aTTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGt                                                                                     <  1:3766606/149‑1 (MQ=255)
                                                                aTTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGt                                                                                     >  1:4541032/1‑149 (MQ=255)
                                                                 tttCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTg                                                                                    >  1:748980/1‑149 (MQ=255)
                                                                  ttCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGa                                                                                   <  1:1327993/149‑1 (MQ=255)
                                                                  ttCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGa                                                                                   >  2:2400161/1‑149 (MQ=255)
                                                                    cAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAAc                                                                                 <  1:1040158/149‑1 (MQ=255)
                                                                     aaCCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACa                                                                                <  2:610301/149‑1 (MQ=255)
                                                                       cccccTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGt                                                                              >  2:560449/1‑149 (MQ=255)
                                                                       cccccTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGt                                                                              >  2:1383584/1‑149 (MQ=255)
                                                                         cccTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTCATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGttt                                                                            >  1:1304532/1‑149 (MQ=255)
                                                                         cccTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGttt                                                                            <  1:2525044/149‑1 (MQ=255)
                                                                             gCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGccc                                                                        <  2:256623/149‑1 (MQ=255)
                                                                                aTACCCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCcggc                                                                     <  2:2752175/149‑1 (MQ=255)
                                                                                aTAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCcggc                                                                     >  2:1964833/1‑149 (MQ=255)
                                                                                aTAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCcggc                                                                     >  2:361971/1‑149 (MQ=255)
                                                                                      tGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCa                                                               >  1:4107738/1‑149 (MQ=255)
                                                                                        aaaaaGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACt                                                             >  1:1407636/1‑149 (MQ=255)
                                                                                         aaaaGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACTATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTg                                                            >  1:468577/1‑149 (MQ=255)
                                                                                            aGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCAt                                                         <  2:2029181/149‑1 (MQ=255)
                                                                                             gCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATc                                                        >  2:2074093/1‑149 (MQ=255)
                                                                                               aGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGc                                                      <  1:825003/149‑1 (MQ=255)
                                                                                                 gATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCTGCGGCGCACTGCATCGCtt                                                    <  2:1136428/149‑1 (MQ=255)
                                                                                                    cTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGAt                                                 <  2:493861/149‑1 (MQ=255)
                                                                                                     tGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATc                                                <  2:3560746/149‑1 (MQ=255)
                                                                                                      ggTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACTATTGGCTGATGAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTCGGTCGTGAACAGAGTGCCCGGCGGCGCACTGCATCGCTAGATCg                                               >  2:2644266/1‑149 (MQ=255)
                                                                                                       gTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGt                                              <  2:2297963/149‑1 (MQ=255)
                                                                                                       gTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGt                                              >  2:4449385/1‑149 (MQ=255)
                                                                                                            cATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGTTAAAc                                         <  2:2235933/149‑1 (MQ=255)
                                                                                                              tCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACTATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGTTAAACAg                                       <  1:3236232/149‑1 (MQ=255)
                                                                                                              tCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGTTAAACAg                                       >  1:1713817/1‑149 (MQ=255)
                                                                                                              tCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGTTAAACAg                                       >  1:2902388/1‑149 (MQ=255)
                                                                                                              tCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGTTAAACAg                                       >  1:300027/1‑149 (MQ=255)
                                                                                                              tCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGTTAAACAg                                       >  1:4731985/1‑149 (MQ=255)
                                                                                                               ccTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGTTAAACAgg                                      >  1:3052187/1‑149 (MQ=255)
                                                                                                                cTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGTTAAACAGGa                                     <  1:1868539/149‑1 (MQ=255)
                                                                                                                cTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGTTAAACAGGa                                     <  1:2530506/149‑1 (MQ=255)
                                                                                                                  gCACGATATGGTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGTTAAACAGGATc                                   <  2:2984549/149‑1 (MQ=255)
                                                                                                                   cACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGTTAAACAGGATCt                                  <  1:2455311/149‑1 (MQ=255)
                                                                                                                     cGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGTTAAACAGGATCTGc                                <  2:1238056/149‑1 (MQ=255)
                                                                                                                      gATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGTTAAACAGGATCTGCa                               <  1:3898040/149‑1 (MQ=255)
                                                                                                                         atGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGTTAAACAGGATCTGCATTg                            <  2:4071949/149‑1 (MQ=255)
                                                                                                                          tGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGTTAAACAGGATCTGCATTGc                           >  2:834236/1‑149 (MQ=255)
                                                                                                                          tGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGAGCGTTAAACAGGATCTGCATTGc                           >  2:1069384/1‑149 (MQ=255)
                                                                                                                           gTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGTTAAACAGGATCTGCATTGCt                          <  1:1095405/149‑1 (MQ=255)
                                                                                                                           gTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGTTAAACAGGATCTGCATTGCt                          <  1:3451214/149‑1 (MQ=255)
                                                                                                                             tGCGCCCACTGCTCACGCGGGACTATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCTGTTGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGTTAAACAGGATCTGCATTGCTtc                        <  1:2332976/149‑1 (MQ=255)
                                                                                                                             tGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGTTAAACAGGATCTGCATTGCTtc                        <  2:4058577/149‑1 (MQ=255)
                                                                                                                               cgcCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGTTAAACAGGATCTGCATTGCTtcct                      <  2:3605088/149‑1 (MQ=255)
                                                                                                                                  ccACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGTTAAACAGGATCTGCATTGCTTCCTCCg                   >  1:851526/1‑149 (MQ=255)
                                                                                                                                   cACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGTTAAACAGGATCTGCATTGCTTCCTCCGc                  >  2:1584959/1‑149 (MQ=255)
                                                                                                                                   cACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGTTAAACAGGATCTGCATTGCTTCCTCCGc                  >  1:1657936/1‑149 (MQ=255)
                                                                                                                                    aCTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGTTAAACAGGATCTGCATTGCTTCCTCCGCa                 >  2:252988/1‑149 (MQ=255)
                                                                                                                                    aCTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGTTAAACAGGATCTGCATTGCTTCCTCCGCa                 <  2:3744841/149‑1 (MQ=255)
                                                                                                                                      tGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGTTAAACAGGATCTGCATTGCTTCCTCCGCATa               <  1:2483762/149‑1 (MQ=255)
                                                                                                                                       gCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGTTAAACAGGATCTGCATTGCTTCCTCCGCATAc              <  2:3932806/149‑1 (MQ=255)
                                                                                                                                       gCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGTTAAACAGGATCTGCATTGCTTCCTCCGCATAc              >  2:857850/1‑149 (MQ=255)
                                                                                                                                         tCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGTTAAACAGGATCTGCATTGCTTCCTCCGCATACCg            >  2:1688768/1‑149 (MQ=255)
                                                                                                                                         tCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGTTAAACAGGATCTGCATTGCTTCCTCCGCATACCg            <  1:4444675/149‑1 (MQ=255)
                                                                                                                                         tCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGTTAAACAGGATCTGCATTGCTTCCTCCGCATACCg            >  1:4086728/1‑149 (MQ=255)
                                                                                                                                           aCGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGTTAAACAGGATCTGCATTGCGTCCTCCGCATACCggg          >  2:2908176/1‑149 (MQ=255)
                                                                                                                                              cgGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGTTAAACAGGATCTGCATTGCTTCCTCCGCATACCGGGCaa       <  2:2468796/149‑1 (MQ=255)
                                                                                                                                                  aCGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGTTAAACAGGATCTGCATTGCTTCCTCCGCATACCGGGCAACCAc   >  2:4351482/1‑149 (MQ=255)
                                                                                                                                                  aCGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGTTAAACAGGATCTGCATTGCTTCCTCCGCATACCGGGCAACCAc   <  1:3343172/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGAAGCGAATTTCCCTGTGCCGGTAAACAGATGTGAATCAAACGTTTTGTCCGCAATACGTAACATTTCAACCCCCTGCAATAACCTGAAAAAGCAGGATCTGGTCGCCATCCTGCACGATATGTTGCGCCCACTGCTCACGCGGGACGATTTGCTGATTAATCGCCAGAGCCGCGCCCGCTTGTCGTTGGTCGAGTTGCTCCAGTAGTTCGTGAACAGTTTGCCCGGCGGCGCACTGCATCGCTTGATCGTTAAACAGGATCTGCATTGCTTCCTCCGCATACCGGGCAACCACT  >  CP009273/4182498‑4182793

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: