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% 674.4 / NA 205R87H (CGC→CAC) rhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (101/104);  total (101/104)
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 207 total aligned reads displayed.

TGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAA  >  CP009273/3959892‑3960180
                                                                                                                                                |                                                                                                                                                
tGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATccaccg                                                                                                                                               >  2:293625/1‑148 (MQ=255)
  aaGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGtt                                                                                                                                             >  1:1909417/1‑148 (MQ=255)
   agcagcACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGttt                                                                                                                                            <  2:1825212/148‑1 (MQ=255)
    gcagcaCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTc                                                                                                                                           >  1:1702272/1‑148 (MQ=255)
     cagcaCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCa                                                                                                                                          >  2:3505798/1‑148 (MQ=255)
       gcaCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAAc                                                                                                                                        <  2:1632251/148‑1 (MQ=255)
            cAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCg                                                                                                                                   >  2:1884226/1‑148 (MQ=255)
              aaGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCa                                                                                                                                 >  1:1551377/1‑148 (MQ=255)
                  gTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTgg                                                                                                                             <  1:2781827/148‑1 (MQ=255)
                   tGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGt                                                                                                                            >  1:2530877/1‑148 (MQ=255)
                   tGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCACTAGCCAAATCCACCGTTTCAACCTCCGCACTGGt                                                                                                                            <  1:2552073/148‑1 (MQ=255)
                    ggCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTg                                                                                                                           <  2:337692/148‑1 (MQ=255)
                      cGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGAt                                                                                                                         <  2:2391394/148‑1 (MQ=255)
                        aaGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATAc                                                                                                                       >  1:1966839/1‑148 (MQ=255)
                        aaGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATAc                                                                                                                       >  1:2536189/1‑148 (MQ=255)
                        aaGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATAc                                                                                                                       <  2:2527860/148‑1 (MQ=255)
                        aaGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAAACCACCGTATCAACCTCCGCACTCGTCATAc                                                                                                                       >  2:1956970/1‑148 (MQ=255)
                          gATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCa                                                                                                                     >  2:1662715/1‑148 (MQ=255)
                            tatCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCAtc                                                                                                                   >  2:1128953/1‑148 (MQ=255)
                             atCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCAtct                                                                                                                  <  1:1000549/148‑1 (MQ=255)
                             atCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCAtct                                                                                                                  >  1:1122533/1‑148 (MQ=255)
                               cTTTGGTGATGGCGTTCTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATctct                                                                                                                <  2:144253/148‑1 (MQ=255)
                               cTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATctct                                                                                                                <  2:2311941/148‑1 (MQ=255)
                                   ggTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTa                                                                                                            <  2:225372/148‑1 (MQ=255)
                                      gATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGa                                                                                                         <  2:2561570/148‑1 (MQ=255)
                                       aTGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGAt                                                                                                        <  1:2504178/148‑1 (MQ=255)
                                          gCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGAttag                                                                                                     <  2:2766105/148‑3 (MQ=255)
                                           cGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgc                                                                                                    >  1:43051/1‑148 (MQ=255)
                                             tACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgccc                                                                                                  >  2:280583/1‑148 (MQ=255)
                                                 ggAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCccgccg                                                                                              >  1:2844192/1‑148 (MQ=255)
                                                   agaTATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGaa                                                                                            <  2:1652936/148‑1 (MQ=255)
                                                   agaTATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGaa                                                                                            <  2:973088/148‑1 (MQ=255)
                                                    gaTATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGaaa                                                                                           <  1:2625935/148‑1 (MQ=255)
                                                        ttGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAaagaag                                                                                       >  2:2012565/1‑148 (MQ=255)
                                                        ttGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAaagaag                                                                                       <  2:1414997/148‑1 (MQ=255)
                                                         tGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAgg                                                                                      >  1:543425/1‑148 (MQ=255)
                                                          gCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGt                                                                                     <  2:2998739/148‑1 (MQ=255)
                                                          gCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGt                                                                                     <  2:1222086/148‑1 (MQ=255)
                                                            aGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGa                                                                                   <  1:874385/148‑1 (MQ=255)
                                                             ggatggatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGaa                                                                                  >  2:1301682/1‑148 (MQ=255)
                                                              gatggatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAAc                                                                                 >  2:2071694/1‑148 (MQ=255)
                                                                tggatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGc                                                                               <  1:1641007/148‑1 (MQ=255)
                                                                 ggatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCt                                                                              <  2:529198/148‑1 (MQ=255)
                                                                  gatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTa                                                                             <  1:2438983/148‑1 (MQ=255)
                                                                   atTTGGTTTCCTCCTTTACTCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTGTCCACTATCCAAATCCACCGTTTCAACCTCCGCTCTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTTAACGCttt                                                                            <  2:3069229/148‑3 (MQ=255)
                                                                       ggTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTg                                                                        <  2:749294/148‑1 (MQ=255)
                                                                       ggTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTg                                                                        >  1:933937/1‑148 (MQ=255)
                                                                            cctccGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgctg                                                                   <  2:3334834/148‑1 (MQ=255)
                                                                            cctccGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgctg                                                                   <  1:1468158/148‑1 (MQ=255)
                                                                             ctccGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgctgc                                                                  >  2:3433208/1‑148 (MQ=255)
                                                                                cGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGa                                                               >  2:2161652/1‑148 (MQ=255)
                                                                                   tCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAg                                                            >  2:758644/1‑148 (MQ=255)
                                                                                      gCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTa                                                         <  1:2269410/148‑1 (MQ=255)
                                                                                      gCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTa                                                         <  2:432468/148‑1 (MQ=255)
                                                                                      gCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTa                                                         <  2:2225802/148‑1 (MQ=255)
                                                                                       cAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTaa                                                        <  1:280583/148‑1 (MQ=255)
                                                                                         gACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTCACg                                                      >  2:2448624/1‑148 (MQ=255)
                                                                                         gACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACg                                                      >  2:2668930/1‑148 (MQ=255)
                                                                                          aCAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGa                                                     <  1:2653790/148‑1 (MQ=255)
                                                                                          aCAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGa                                                     >  1:3185514/1‑148 (MQ=255)
                                                                                             gCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGt                                                  <  1:2573546/148‑1 (MQ=255)
                                                                                              cTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGtt                                                 >  2:424206/1‑148 (MQ=255)
                                                                                              cTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGtt                                                 >  1:1635997/1‑148 (MQ=255)
                                                                                               tCCTACCTCGCCGGTCCTTATGACTTATACGTTTCCCCTAGCCAAATCCACGGTTTCAACCTCCGCACTGGTGATACCATCTCTGTTAAGATTCGCCCGCCGAAAGAAGGTGAACGTTATTTTGCGCTGCTGAAAGTTAACGATGTTa                                                <  2:1511/148‑1 (MQ=255)
                                                                                                ccTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTaa                                               <  1:1678184/148‑1 (MQ=255)
                                                                                                 cTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAAc                                              <  1:2549574/148‑1 (MQ=255)
                                                                                                  tACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACt                                             <  1:3140218/148‑1 (MQ=255)
                                                                                                  tACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACt                                             >  1:709157/1‑148 (MQ=255)
                                                                                                       cGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGAc                                        <  2:1997001/148‑1 (MQ=255)
                                                                                                         ccGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACaa                                      <  2:704064/148‑1 (MQ=255)
                                                                                                         ccGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACaa                                      <  2:210262/148‑1 (MQ=255)
                                                                                                         ccGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACaa                                      <  1:830370/148‑1 (MQ=255)
                                                                                                             tCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGACCGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCt                                  >  2:3230409/1‑148 (MQ=255)
                                                                                                                 gatgaCATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGaaa                              >  2:2573851/1‑148 (MQ=255)
                                                                                                                  atgaCATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCTCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGaaaa                             <  1:293625/148‑1 (MQ=255)
                                                                                                                   tgaCATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAAc                            >  1:3492934/1‑148 (MQ=255)
                                                                                                                   tgaCATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAAc                            >  1:1727536/1‑148 (MQ=255)
                                                                                                                     aCATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGc                          >  1:2143863/1‑148 (MQ=255)
                                                                                                                     aCATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGc                          >  1:1782695/1‑148 (MQ=255)
                                                                                                                        tCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCg                       >  2:678802/1‑148 (MQ=255)
                                                                                                                        tCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCg                       >  2:2478677/1‑148 (MQ=255)
                                                                                                                              tttCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaacaa                 >  1:3002169/1‑148 (MQ=255)
                                                                                                                               ttCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACaaa                >  1:3337655/1‑148 (MQ=255)
                                                                                                                                tCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACaaaa               >  1:1087424/1‑148 (MQ=255)
                                                                                                                                tCCCCTACCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGTTGAACGCTATTTTGCGCTGCTTAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACaaaa               <  2:1081055/148‑1 (MQ=255)
                                                                                                                                 ccccTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAAt              <  2:3524086/148‑1 (MQ=255)
                                                                                                                                 ccccTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAAt              <  2:2338425/148‑1 (MQ=255)
                                                                                                                                 ccccTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAAt              >  1:1854162/1‑148 (MQ=255)
                                                                                                                                 ccccTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAAt              <  2:2821605/148‑1 (MQ=255)
                                                                                                                                     tAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCctc          >  1:2411707/1‑148 (MQ=255)
                                                                                                                                      ccccAAATCCACCTTTTCAACCCCCGCAGTGGTGATACCATCTCTGGTAATTTTCGCCCGCCTAAAGACGGTGAAGGATATTTTGCGCTGCTTAAAGTTAACCAATTTAACTTCGAAAAACCTCAAAACGCCCGCAACAAAATCctct         <  1:2233340/146‑1 (MQ=255)
                                                                                                                                      aGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCctct         >  1:2324840/1‑148 (MQ=255)
                                                                                                                                       gCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCtt        >  1:1854781/1‑148 (MQ=255)
                                                                                                                                       gCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACACTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCtt        >  2:355779/1‑148 (MQ=255)
                                                                                                                                         cAAATCCACCGTTACAACCTCCGCACTGGTGATACAATCTCTGGTAAGATTCGCCCGCCGAATGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTAAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTg      <  2:619952/148‑1 (MQ=255)
                                                                                                                                          aaaTCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTga     >  2:2652627/1‑148 (MQ=255)
                                                                                                                                            aTCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTgaga   >  1:3147345/1‑148 (MQ=255)
                                                                                                                                            aTCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTgaga   >  2:879924/1‑148 (MQ=255)
                                                                                                                                             tCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGTTAAGAGTCGCCCGCCGAACGAAGGTGAACGCTATTTGGCGCTGCTTAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGaa  <  2:11566/148‑1 (MQ=255)
                                                                                                                                             tCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGaa  >  1:415361/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
TGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAA  >  CP009273/3959892‑3960180

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: