Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,960,036 G→A 100% R87H (CGC→CAC)  rho → transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,960,0360GA100.0% 809.7 / NA 245R87H (CGC→CAC) rhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (122/123);  total (122/123)
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 245 total aligned reads displayed.

TGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAAC  >  CP009273/3959892‑3960181
                                                                                                                                                |                                                                                                                                                 
tGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGt                                                                                                                                               >  2:185584/1‑149 (MQ=255)
 gAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGtt                                                                                                                                              >  1:3586221/1‑149 (MQ=255)
   agcagcACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTc                                                                                                                                            <  2:2571114/149‑1 (MQ=255)
   agcagcACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTc                                                                                                                                            >  1:3302853/1‑149 (MQ=255)
   agcagcACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTc                                                                                                                                            >  1:538701/1‑149 (MQ=255)
   agcagcACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGAGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTACCGCAGAAAGCTCCAACCTCGCCGGTCCAGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTc                                                                                                                                            >  1:2838567/1‑149 (MQ=255)
       gcaCGCAAAGAGTCCCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAAcc                                                                                                                                        >  2:117121/1‑149 (MQ=255)
             aaaGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCa                                                                                                                                  >  2:3250688/1‑149 (MQ=255)
             aaaGAGGGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCa                                                                                                                                  <  2:2228624/149‑1 (MQ=255)
              aaGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCAc                                                                                                                                 >  2:553233/1‑149 (MQ=255)
              aaGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCAc                                                                                                                                 >  2:3677076/1‑149 (MQ=255)
                 agTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTgg                                                                                                                              >  2:494990/1‑149 (MQ=255)
                  gTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGt                                                                                                                             <  1:3739985/149‑1 (MQ=255)
                   tGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTg                                                                                                                            <  2:420723/149‑1 (MQ=255)
                    ggCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGa                                                                                                                           >  1:1180202/1‑149 (MQ=255)
                       gAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATAc                                                                                                                        >  2:3156133/1‑149 (MQ=255)
                         aGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCa                                                                                                                      >  1:1432417/1‑149 (MQ=255)
                          gATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCAt                                                                                                                     <  1:1924954/149‑1 (MQ=255)
                          gATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCAt                                                                                                                     <  2:262539/149‑1 (MQ=255)
                           atatCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCAtc                                                                                                                    <  2:3305583/149‑1 (MQ=255)
                           atatCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCAtc                                                                                                                    >  2:1138289/1‑149 (MQ=255)
                            tatCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCAtct                                                                                                                   >  2:4030276/1‑149 (MQ=255)
                            tatCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCAtct                                                                                                                   <  1:2574861/149‑1 (MQ=255)
                            tatCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCAtct                                                                                                                   >  1:1776475/1‑149 (MQ=255)
                            tatCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCAtct                                                                                                                   >  1:2048593/1‑149 (MQ=255)
                              tCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATctct                                                                                                                 <  1:1689358/149‑1 (MQ=255)
                               cTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTg                                                                                                                >  2:2387569/1‑149 (MQ=255)
                                tttGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTgg                                                                                                               >  2:2942864/1‑149 (MQ=255)
                                 ttGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGt                                                                                                              <  1:2223803/149‑1 (MQ=255)
                                   ggTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTaa                                                                                                            >  2:2735649/1‑149 (MQ=255)
                                    gTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAg                                                                                                           >  2:3805507/1‑149 (MQ=255)
                                      gATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGAt                                                                                                         >  2:200796/1‑149 (MQ=255)
                                      gATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGAt                                                                                                         >  1:563754/1‑149 (MQ=255)
                                        tGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTc                                                                                                       <  1:188590/149‑1 (MQ=255)
                                        tGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTc                                                                                                       >  1:992878/1‑149 (MQ=255)
                                          gCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgc                                                                                                     >  1:2686319/1‑149 (MQ=255)
                                           cGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgcc                                                                                                    <  2:2170256/149‑1 (MQ=255)
                                             tACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgcccg                                                                                                  >  2:1572145/1‑149 (MQ=255)
                                               cTGGAGATATTGCAGGATGGATTTTGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGACTGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgcccgcc                                                                                                <  1:477986/149‑1 (MQ=255)
                                                   agaTATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGaaa                                                                                            <  1:3004357/149‑1 (MQ=255)
                                                    gaTATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAaag                                                                                           >  1:3264240/1‑149 (MQ=255)
                                                        ttGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAgg                                                                                       >  1:3688887/1‑149 (MQ=255)
                                                        ttGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAgg                                                                                       >  1:3619885/1‑149 (MQ=255)
                                                            aGGATGGATTTGGTTTCCTCCGGTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGaa                                                                                   <  2:2616834/149‑1 (MQ=255)
                                                              gatggatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACg                                                                                 >  2:3507223/1‑149 (MQ=255)
                                                                tggatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCt                                                                               >  2:1705094/1‑149 (MQ=255)
                                                                   atTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTAtt                                                                            <  2:2991487/149‑1 (MQ=255)
                                                                    tttGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTAttt                                                                           <  2:133252/149‑1 (MQ=255)
                                                                    tttGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTAttt                                                                           >  2:2813554/1‑149 (MQ=255)
                                                                    tttGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTAttt                                                                           <  2:660589/149‑1 (MQ=255)
                                                                       ggTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTgc                                                                        >  1:2153911/1‑149 (MQ=255)
                                                                       ggTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTgc                                                                        >  1:1248858/1‑149 (MQ=255)
                                                                       ggTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTgc                                                                        <  1:728510/149‑1 (MQ=255)
                                                                        gTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTgcg                                                                       <  1:1317827/149‑1 (MQ=255)
                                                                         tttCCCCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTgcgc                                                                      <  2:180250/149‑1 (MQ=255)
                                                                            cctccGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgctgc                                                                   <  1:3891193/149‑1 (MQ=255)
                                                                             ctccGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgctgct                                                                  >  1:3126515/1‑149 (MQ=255)
                                                                               ccGTTCCGCAGAAAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGAACTGGTGATACCATATCTGGTAAGAATCGCCCGCCGAAAGAAGGTGAACGCTATTTTCCGCTGCTGa                                                                >  2:927888/1‑149 (MQ=255)
                                                                                 gTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGaaa                                                              <  2:661975/149‑1 (MQ=255)
                                                                                 gTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGaaa                                                              <  2:2906289/149‑1 (MQ=255)
                                                                                 gTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGaaa                                                              >  1:291326/1‑149 (MQ=255)
                                                                                   tCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGt                                                            <  1:3677076/149‑1 (MQ=255)
                                                                                   tCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACAGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGt                                                            <  1:3746654/149‑1 (MQ=255)
                                                                                    ccGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGtt                                                           <  2:2876137/149‑1 (MQ=255)
                                                                                          aCAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGaa                                                     <  1:164718/149‑1 (MQ=255)
                                                                                          aCAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGaa                                                     <  1:511089/149‑1 (MQ=255)
                                                                                            aGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGt                                                   <  2:3946017/149‑1 (MQ=255)
                                                                                            aGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGt                                                   <  1:1504853/149‑1 (MQ=255)
                                                                                              cTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTa                                                 <  1:1036898/149‑1 (MQ=255)
                                                                                              cTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTa                                                 <  1:2881664/149‑1 (MQ=255)
                                                                                               tCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTaa                                                >  2:3522093/1‑149 (MQ=255)
                                                                                                  tACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACtt                                             >  2:3206991/1‑149 (MQ=255)
                                                                                                  tACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACtt                                             <  2:3886251/149‑1 (MQ=255)
                                                                                                    ccTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCg                                           >  2:239353/1‑149 (MQ=255)
                                                                                                      tCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCTTCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTTAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGAc                                         <  1:2337155/149‑1 (MQ=255)
                                                                                                      tCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGAc                                         <  2:1033542/149‑1 (MQ=255)
                                                                                                       cGCCGGTCCTGATGACATCTACGTGTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACa                                        <  2:782735/149‑1 (MQ=255)
                                                                                                         ccGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACaaa                                      <  1:1649019/149‑1 (MQ=255)
                                                                                                         ccGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAAAGCTATTTTGCGCTGCTGAAAGTTAACTAAGTTAACTTCGACaaa                                      >  1:3781595/1‑149 (MQ=255)
                                                                                                          cGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAAc                                     <  2:269239/149‑1 (MQ=255)
                                                                                                           ggTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAAcc                                    >  1:3397206/1‑149 (MQ=255)
                                                                                                            gTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCt                                   >  1:1164271/1‑149 (MQ=255)
                                                                                                             tCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTg                                  >  1:3864604/1‑149 (MQ=255)
                                                                                                             tCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTg                                  <  1:2624881/149‑1 (MQ=255)
                                                                                                                      cATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGccc                         <  2:563754/149‑1 (MQ=255)
                                                                                                                       aTCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCg                        >  2:1404295/1‑149 (MQ=255)
                                                                                                                         cTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGca                      >  1:1071185/1‑149 (MQ=255)
                                                                                                                         cTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGca                      >  1:1033176/1‑149 (MQ=255)
                                                                                                                           aCGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaac                    >  1:2318577/1‑149 (MQ=255)
                                                                                                                                    cTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCctc           >  1:234685/1‑149 (MQ=255)
                                                                                                                                     tAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCctct          >  1:2057299/1‑149 (MQ=255)
                                                                                                                                     tAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCctct          >  1:2711549/1‑149 (MQ=255)
                                                                                                                                     tAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCctct          <  2:2153911/149‑1 (MQ=255)
                                                                                                                                       gCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCttt        <  2:2690848/149‑1 (MQ=255)
                                                                                                                                       gCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCttt        <  2:2309541/149‑1 (MQ=255)
                                                                                                                                        ccAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTg       >  1:1944494/1‑149 (MQ=255)
                                                                                                                                        ccAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTg       >  1:3981046/1‑149 (MQ=255)
                                                                                                                                          aaaTCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTgag     >  1:297233/1‑149 (MQ=255)
                                                                                                                                           aaTCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTgaga    <  2:3617214/149‑1 (MQ=255)
                                                                                                                                             tCCACCGTTACAACCTCCGCACGGTTGATACCTTCTCTGGTATGATTCGCCCGCCTAATTAAGGTGAACGCTATTTTGCGCTGCTGAACGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAAc  <  2:1432417/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
TGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAAC  >  CP009273/3959892‑3960181

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: