Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,383,318 G→T 8.1% V189V (GTC→GTA queG ← epoxyqueuosine reductase, cobalamine‑stimulated; queosine biosynthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,383,3180GT 8.1% 1166.9 / 58.0 309V189V (GTC→GTAqueGepoxyqueuosine reductase, cobalamine‑stimulated; queosine biosynthesis
Reads supporting (aligned to +/- strand):  ref base G (153/131);  new base T (19/6);  total (172/137)
Fisher's exact test for biased strand distribution p-value = 3.65e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 313 total aligned reads displayed.

TTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAG  >  CP009273/4383170‑4383466
                                                                                                                                                    |                                                                                                                                                    
  aaCGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGAc                                                                                                                                                    <  2:3038877/149‑1 (MQ=255)
     ggTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTgg                                                                                                                                                 >  1:509742/1‑149 (MQ=255)
      gTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCTACTGGt                                                                                                                                                <  1:2914538/149‑1 (MQ=255)
      gTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGt                                                                                                                                                >  2:2880843/1‑149 (MQ=255)
      gTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGt                                                                                                                                                <  2:180915/149‑1 (MQ=255)
        cGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGCTGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCACATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTg                                                                                                                                              >  2:4485641/1‑149 (MQ=255)
        cGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTg                                                                                                                                              >  1:4368768/1‑149 (MQ=255)
          cAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGAt                                                                                                                                            >  2:1254886/1‑149 (MQ=255)
             ctctTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCa                                                                                                                                         <  1:4049515/149‑1 (MQ=255)
              tctTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCAc                                                                                                                                        <  1:3011832/149‑1 (MQ=255)
                ttCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCTACTGGTTGATCCACgg                                                                                                                                      <  1:1195894/149‑1 (MQ=255)
                    gggATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGTGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAg                                                                                                                                  >  1:230167/1‑149 (MQ=255)
                    gggATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAg                                                                                                                                  >  1:2957669/1‑149 (MQ=255)
                      gATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCTACTGGTTGATCCACGGGCAGCg                                                                                                                                <  2:866270/149‑1 (MQ=255)
                      gATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCg                                                                                                                                >  2:1647305/1‑149 (MQ=255)
                       aTCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCgg                                                                                                                               >  2:3775058/1‑149 (MQ=255)
                       aTCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCgg                                                                                                                               >  2:1157722/1‑149 (MQ=255)
                         cGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGaa                                                                                                                             <  2:4329106/149‑1 (MQ=255)
                           ccccTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAata                                                                                                                           <  2:1707171/149‑1 (MQ=255)
                            cccTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAatat                                                                                                                          <  2:4993350/149‑1 (MQ=255)
                            cccTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAatat                                                                                                                          <  1:3626462/149‑1 (MQ=255)
                              cTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCACATCCTTCCTCGACTGGTTGATCCACGGGCAGAGGAATATCg                                                                                                                        >  1:3149401/1‑149 (MQ=255)
                               ttCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGa                                                                                                                       <  1:348069/149‑1 (MQ=255)
                                tCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGAc                                                                                                                      >  2:4901030/1‑149 (MQ=255)
                                tCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGAc                                                                                                                      >  2:2578856/1‑149 (MQ=255)
                                 cAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCTACTGGTTGATCCACGGGCAGCGGAATATCGAcc                                                                                                                     >  1:2554436/1‑149 (MQ=255)
                                  aaGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCa                                                                                                                    >  2:4755113/1‑149 (MQ=255)
                                   aGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAg                                                                                                                   <  1:1627004/149‑1 (MQ=255)
                                    gTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGc                                                                                                                  >  2:822191/1‑149 (MQ=255)
                                     ttCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCTACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCa                                                                                                                 >  2:179515/1‑149 (MQ=255)
                                      tCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCaa                                                                                                                <  2:4804897/149‑1 (MQ=255)
                                       cGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAAt                                                                                                               <  2:4096800/149‑1 (MQ=255)
                                       cGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAAt                                                                                                               <  2:1479957/149‑1 (MQ=255)
                                       cGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAAt                                                                                                               <  1:1726502/149‑1 (MQ=255)
                                        gATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAAtt                                                                                                              <  1:5219707/149‑1 (MQ=255)
                                           ggTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGc                                                                                                           <  1:889945/149‑1 (MQ=255)
                                            gTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGcc                                                                                                          <  1:4917570/149‑1 (MQ=255)
                                            gTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGcc                                                                                                          >  2:2460383/1‑149 (MQ=255)
                                              gagaTAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTa                                                                                                        >  1:4294082/1‑149 (MQ=255)
                                                gaTAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTaaa                                                                                                      <  1:1007245/149‑1 (MQ=255)
                                                    agagATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAaaga                                                                                                  >  1:3107336/1‑149 (MQ=255)
                                                     gagaTACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAaagaa                                                                                                 >  1:2355711/1‑149 (MQ=255)
                                                       gaTACAGCGGCGAGAATCGACGGTATTTGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAagaaga                                                                                               <  1:2813485/149‑1 (MQ=255)
                                                        aTACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAgaagaa                                                                                              <  2:3091938/149‑1 (MQ=255)
                                                         tACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAgaac                                                                                             <  2:2409279/149‑1 (MQ=255)
                                                         tACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAgaac                                                                                             >  2:1947876/1‑149 (MQ=255)
                                                            aGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAgaacgaa                                                                                          <  2:1024006/149‑1 (MQ=255)
                                                            aGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAgaacgaa                                                                                          >  2:2133461/1‑149 (MQ=255)
                                                              cggcgAGCTTCGACGGTATATGGCTCGACGTTGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCTGCGGAATATCGACCAGCATTTCGCCTAAAAAGAAGAACGAAcc                                                                                        <  1:471286/149‑1 (MQ=255)
                                                                    gCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCTACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTc                                                                                  >  2:4192843/1‑149 (MQ=255)
                                                                    gCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTc                                                                                  >  1:1564872/1‑149 (MQ=255)
                                                                      aTCGACGGTATCTGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTGGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTcgc                                                                                <  1:3240816/149‑1 (MQ=255)
                                                                      aTCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTcgc                                                                                <  1:617861/149‑1 (MQ=255)
                                                                      aTCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTcgc                                                                                <  1:4987235/149‑1 (MQ=255)
                                                                      aTCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTcgc                                                                                >  1:3802155/1‑149 (MQ=255)
                                                                       tCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTcgcg                                                                               >  2:2589953/1‑149 (MQ=255)
                                                                       tCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTcgcg                                                                               >  2:4711927/1‑149 (MQ=255)
                                                                          aCGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGAtt                                                                            >  2:5201516/1‑149 (MQ=255)
                                                                          aCGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGAtt                                                                            >  2:1667778/1‑149 (MQ=255)
                                                                            ggTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGa                                                                          >  2:736222/1‑149 (MQ=255)
                                                                             gTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAg                                                                         <  1:3339926/149‑1 (MQ=255)
                                                                                   ggCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAag                                                                   <  2:3689661/149‑1 (MQ=255)
                                                                                   ggCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAag                                                                   <  1:455805/149‑1 (MQ=255)
                                                                                    gCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAagt                                                                  <  2:3214868/149‑1 (MQ=255)
                                                                                           gaTGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCt                                                           <  1:4011377/149‑1 (MQ=255)
                                                                                                 aCCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTg                                                     >  2:703184/1‑149 (MQ=255)
                                                                                                  ccGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGt                                                    <  1:2476214/149‑1 (MQ=255)
                                                                                                  ccGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGt                                                    >  2:4778907/1‑149 (MQ=255)
                                                                                                    ggTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTcc                                                  <  1:2114996/149‑1 (MQ=255)
                                                                                                     gTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCa                                                 <  1:5166717/149‑1 (MQ=255)
                                                                                                     gTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCa                                                 <  1:1448543/149‑1 (MQ=255)
                                                                                                     gTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCa                                                 <  2:193014/149‑1 (MQ=255)
                                                                                                        gggCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGcc                                              <  1:3954801/149‑1 (MQ=255)
                                                                                                         ggCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCg                                             >  2:3446295/1‑149 (MQ=255)
                                                                                                          gCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGa                                            >  2:1407128/1‑149 (MQ=255)
                                                                                                             aaTCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGcc                                         >  2:85281/1‑149 (MQ=255)
                                                                                                               tCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCa                                       >  1:1597301/1‑149 (MQ=255)
                                                                                                                  tCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCt                                    >  1:4783902/1‑149 (MQ=255)
                                                                                                                  tCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCt                                    >  1:3766886/1‑149 (MQ=255)
                                                                                                                    aTACAGGCCACGCATTTGCCGCATCCTTCCTCTACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCttt                                  >  2:3707185/1‑149 (MQ=255)
                                                                                                                        aGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCa                              >  2:3933070/1‑149 (MQ=255)
                                                                                                                         ggCCACGCATTTGCCGCATCCTTCCTCTACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAg                             >  2:2193097/1‑149 (MQ=255)
                                                                                                                         ggCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAg                             >  1:4497775/1‑149 (MQ=255)
                                                                                                                           ccACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCt                           <  1:1435036/149‑1 (MQ=255)
                                                                                                                             aCGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTaa                         <  1:2175265/149‑1 (MQ=255)
                                                                                                                             aCGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTaa                         >  2:3830022/1‑149 (MQ=255)
                                                                                                                               gCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACg                       <  2:2329443/149‑1 (MQ=255)
                                                                                                                                 aTTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACggg                     <  2:2138749/149‑1 (MQ=255)
                                                                                                                                     gCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCt                 >  1:1876023/1‑149 (MQ=255)
                                                                                                                                      ccGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTc                <  2:2202862/149‑1 (MQ=255)
                                                                                                                                       cGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCg               >  2:1291766/1‑149 (MQ=255)
                                                                                                                                        gCATCCTTCCTCTACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGCTAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCga              >  1:1800226/1‑149 (MQ=255)
                                                                                                                                         cATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCgag             >  1:691572/1‑149 (MQ=255)
                                                                                                                                          aTCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCgaga            >  1:4438704/1‑149 (MQ=255)
                                                                                                                                               tcctCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAgg       >  2:622557/1‑149 (MQ=255)
                                                                                                                                                cctCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGgc      <  2:3498441/149‑1 (MQ=255)
                                                                                                                                                cctCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGgc      <  1:3576551/149‑1 (MQ=255)
                                                                                                                                                 ctCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGgcg     <  2:3984563/149‑1 (MQ=255)
                                                                                                                                                   cGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCa   <  1:3109006/149‑1 (MQ=255)
                                                                                                                                                    gACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAg  >  2:1327317/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAG  >  CP009273/4383170‑4383466

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: