Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 690,099 G→T 5.5% P105T (CCG→ACG)  miaB ← tRNA‑i(6)A37 methylthiotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273690,0990GT 5.5% 1376.5 / 34.3 344P105T (CCG→ACG) miaBtRNA‑i(6)A37 methylthiotransferase
Reads supporting (aligned to +/- strand):  ref base G (160/165);  new base T (7/12);  total (167/177)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 348 total aligned reads displayed.

CATGATGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAG  >  CP009273/689950‑690248
                                                                                                                                                     |                                                                                                                                                     
 atgatgGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGc                                                                                                                                                      >  2:3322948/1‑150 (MQ=255)
   gatgGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGccc                                                                                                                                                    >  1:932460/1‑150 (MQ=255)
      ggAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCaaa                                                                                                                                                 <  2:2723386/150‑1 (MQ=255)
          aCAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAaata                                                                                                                                             >  1:1402036/1‑150 (MQ=255)
           cAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAaataa                                                                                                                                            <  2:3194302/150‑1 (MQ=255)
               cgcgGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATc                                                                                                                                        >  1:4547236/1‑150 (MQ=255)
               cgcgGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATc                                                                                                                                        <  2:2916958/150‑1 (MQ=255)
               cgcgGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATc                                                                                                                                        <  1:2476580/150‑1 (MQ=255)
                 cgGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGa                                                                                                                                      <  1:4860104/150‑1 (MQ=255)
                     cGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATa                                                                                                                                  >  1:4741977/1‑150 (MQ=255)
                       gCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGTCCCAAAAATAATATCGACATAGt                                                                                                                                <  2:2979666/150‑1 (MQ=255)
                         ccTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTgg                                                                                                                              <  2:113920/150‑1 (MQ=255)
                         ccTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTgg                                                                                                                              <  2:1306796/150‑1 (MQ=255)
                         ccTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTgg                                                                                                                              >  2:903199/1‑150 (MQ=255)
                              ggCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGgcgcg                                                                                                                         <  1:2675565/150‑1 (MQ=255)
                               gcgcgcGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGgcgcgc                                                                                                                        <  1:3122964/150‑1 (MQ=255)
                                      gTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGaa                                                                                                                 <  2:3632964/150‑1 (MQ=255)
                                      gTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGaa                                                                                                                 >  2:2110280/1‑150 (MQ=255)
                                       ttCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAAt                                                                                                                >  1:2775162/1‑150 (MQ=255)
                                        tCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAAtg                                                                                                               >  2:1942356/1‑150 (MQ=255)
                                          cGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGACtgtg                                                                                                             >  2:3385599/1‑150 (MQ=255)
                                          cGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAAtgtg                                                                                                             >  2:3281425/1‑150 (MQ=255)
                                            gTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGACTGTGCt                                                                                                           >  2:4619849/1‑150 (MQ=255)
                                              aGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCg                                                                                                         >  1:282266/1‑150 (MQ=255)
                                                 cGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCct                                                                                                      <  1:2922053/150‑1 (MQ=255)
                                                  ggTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCctt                                                                                                     >  2:4737635/1‑150 (MQ=255)
                                                  ggTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCctt                                                                                                     >  1:651716/1‑150 (MQ=255)
                                                      aaaCTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTg                                                                                                 >  1:120235/1‑150 (MQ=255)
                                                          ttCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGAt                                                                                             <  1:579397/150‑1 (MQ=255)
                                                           gctcGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGTCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATg                                                                                            <  1:3106338/149‑1 (MQ=255)
                                                            ctcGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGc                                                                                           <  2:4833437/150‑1 (MQ=255)
                                                             tcGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGcc                                                                                          >  1:3034952/1‑150 (MQ=255)
                                                              cGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCa                                                                                         >  2:4991136/1‑150 (MQ=255)
                                                              cGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCa                                                                                         <  1:217112/150‑1 (MQ=255)
                                                               gATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCAc                                                                                        >  2:2344785/1‑150 (MQ=255)
                                                                aTTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACg                                                                                       >  2:2671349/1‑150 (MQ=255)
                                                                aTTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACg                                                                                       >  1:3223174/1‑150 (MQ=255)
                                                                aTTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACg                                                                                       >  1:771607/1‑150 (MQ=255)
                                                                aTTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACg                                                                                       >  1:2457422/1‑150 (MQ=255)
                                                                 tttCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGc                                                                                      >  2:3695654/1‑150 (MQ=255)
                                                                  ttCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCa                                                                                     <  2:395234/150‑1 (MQ=255)
                                                                   tCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAg                                                                                    <  1:1027997/150‑1 (MQ=255)
                                                                   tCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAg                                                                                    >  2:2234494/1‑150 (MQ=255)
                                                                     cGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGcc                                                                                  <  1:458100/150‑1 (MQ=255)
                                                                      gggAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCa                                                                                 >  1:144690/1‑150 (MQ=255)
                                                                       ggAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCAc                                                                                <  2:5032665/150‑1 (MQ=255)
                                                                        gAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATCTCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCAcc                                                                               <  1:3846005/150‑1 (MQ=255)
                                                                        gAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCAcc                                                                               <  1:4466616/150‑1 (MQ=255)
                                                                        gAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCAcc                                                                               <  1:3103817/150‑1 (MQ=255)
                                                                           gCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGAc                                                                            >  1:2192760/1‑150 (MQ=255)
                                                                             tGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGc                                                                          >  2:4549097/1‑150 (MQ=255)
                                                                                tatCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGa                                                                       <  1:1082743/150‑1 (MQ=255)
                                                                                  tCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATa                                                                     <  2:3005228/150‑1 (MQ=255)
                                                                                   cTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATaa                                                                    >  2:4583099/1‑150 (MQ=255)
                                                                                     aCAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATc                                                                  <  1:2543734/150‑1 (MQ=255)
                                                                                     aCAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATc                                                                  <  2:822034/150‑1 (MQ=255)
                                                                                       aaCCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAg                                                                <  1:5014779/150‑1 (MQ=255)
                                                                                          cGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTc                                                             >  2:2809301/1‑150 (MQ=255)
                                                                                          cGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTc                                                             <  1:3617818/150‑1 (MQ=255)
                                                                                           gggCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCt                                                            <  1:2607906/150‑1 (MQ=255)
                                                                                            ggCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTg                                                           <  1:426516/150‑1 (MQ=255)
                                                                                                gCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGAtt                                                       <  1:3718922/150‑1 (MQ=255)
                                                                                                 cGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGAttc                                                      <  1:3486926/150‑1 (MQ=255)
                                                                                                      gCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTtctc                                                 >  2:3167198/1‑150 (MQ=255)
                                                                                                      gCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTtctc                                                 >  2:3077420/1‑150 (MQ=255)
                                                                                                       ccACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTctct                                                <  2:3627266/150‑1 (MQ=255)
                                                                                                          cgcacgGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCtttt                                             >  1:4175975/1‑150 (MQ=255)
                                                                                                          cgcacgGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCtttt                                             <  1:2535702/150‑1 (MQ=255)
                                                                                                           gcacgGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTa                                            <  1:1591641/150‑1 (MQ=255)
                                                                                                             acgGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAAc                                          <  1:2202803/150‑1 (MQ=255)
                                                                                                             acgGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAAc                                          >  1:1651140/1‑150 (MQ=255)
                                                                                                               ggAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAg                                        <  2:561158/150‑1 (MQ=255)
                                                                                                               ggAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAg                                        >  2:1218368/1‑150 (MQ=255)
                                                                                                               ggAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAg                                        <  2:3649571/150‑1 (MQ=255)
                                                                                                                gAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGt                                       >  1:2297899/1‑150 (MQ=255)
                                                                                                                    tGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTcc                                   <  2:1395121/150‑1 (MQ=255)
                                                                                                                     gATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCa                                  >  1:3879360/1‑150 (MQ=255)
                                                                                                                     gATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCa                                  >  2:4443353/1‑150 (MQ=255)
                                                                                                                      atcatcTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAg                                 >  2:3774673/1‑150 (MQ=255)
                                                                                                                       tcatcTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGc                                <  1:2014562/150‑1 (MQ=255)
                                                                                                                       tcatcTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGc                                >  2:643465/1‑150 (MQ=255)
                                                                                                                        catcTCCGGCAGACGGTGCAGCGTTTGCGTCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCg                               >  2:2951194/1‑150 (MQ=255)
                                                                                                                         atcTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGa                              >  1:1871811/1‑150 (MQ=255)
                                                                                                                         atcTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGa                              <  2:4162650/150‑1 (MQ=255)
                                                                                                                           ctcCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGAcc                            <  1:4793591/150‑1 (MQ=255)
                                                                                                                                  aGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGa                     >  1:2416063/1‑150 (MQ=255)
                                                                                                                                    aCGGTGCAGCGTTTGCGTCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATg                   >  2:775964/1‑150 (MQ=255)
                                                                                                                                      ggTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGa                 >  1:4347297/1‑150 (MQ=255)
                                                                                                                                       gTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGaa                >  2:679075/1‑150 (MQ=255)
                                                                                                                                        tGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAg               >  1:2063076/1‑150 (MQ=255)
                                                                                                                                        tGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAg               <  1:2766281/150‑1 (MQ=255)
                                                                                                                                         gCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGa              >  2:2907213/1‑150 (MQ=255)
                                                                                                                                          cAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGAc             >  1:3785352/1‑150 (MQ=255)
                                                                                                                                            gCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACtt           <  2:871703/150‑1 (MQ=255)
                                                                                                                                               tttGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACttttt        >  2:3735373/1‑150 (MQ=255)
                                                                                                                                                ttGCGGCCCAAAAATATTTTCGAGATAGTGGGCGCTGTGGCGTTTGTGCTCGCCTTCTTGCGATGCCACGCTGCCACCGAGGCCGTTTATCAGGTCTGTATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACttttta       <  2:938235/150‑2 (MQ=255)
                                                                                                                                                 tGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTcc      >  2:2432609/1‑150 (MQ=255)
                                                                                                                                                  gCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCt     >  2:1605117/1‑150 (MQ=255)
                                                                                                                                                   cGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTg    >  2:1089241/1‑150 (MQ=255)
                                                                                                                                                    ggCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGa   <  2:2581585/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CATGATGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAG  >  CP009273/689950‑690248

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: