Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 689,079 C→A 14.2% E445* (GAA→TAA)  miaB ← tRNA‑i(6)A37 methylthiotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273689,0790CA14.2% 1009.7 / 101.8 313E445* (GAA→TAA) miaBtRNA‑i(6)A37 methylthiotransferase
Reads supporting (aligned to +/- strand):  ref base C (127/141);  new base A (21/24);  total (148/165)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 315 total aligned reads displayed.

TAAAGATGGGGATAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTACGCCCGGAAAGCTCCATGATGCTC  >  CP009273/688930‑689227
                                                                                                                                                     |                                                                                                                                                    
   aGATGGGGATAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATc                                                                                                                                                   <  1:1250722/150‑1 (MQ=255)
   aGATGGGGATAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATc                                                                                                                                                   <  2:2154616/150‑1 (MQ=255)
     aTGGGGATAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTAATCtt                                                                                                                                                 <  1:4640813/150‑1 (MQ=255)
         ggATAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGt                                                                                                                                             >  1:972203/1‑150 (MQ=255)
          gATAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTa                                                                                                                                            >  2:869074/1‑150 (MQ=255)
          gATAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTa                                                                                                                                            <  2:4189367/150‑1 (MQ=255)
           aTAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTAc                                                                                                                                           <  1:1029959/150‑1 (MQ=255)
           aTAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTAc                                                                                                                                           <  2:2359338/150‑1 (MQ=255)
           aTAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTAc                                                                                                                                           <  1:862450/150‑1 (MQ=255)
           aTAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTAc                                                                                                                                           <  1:3756508/150‑1 (MQ=255)
           aTAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTAc                                                                                                                                           <  1:3756419/150‑1 (MQ=255)
           aTAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTAc                                                                                                                                           <  1:3601354/150‑1 (MQ=255)
            tAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACg                                                                                                                                          <  1:1931515/150‑1 (MQ=255)
              aaGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGca                                                                                                                                        >  2:279641/1‑150 (MQ=255)
              aaGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGca                                                                                                                                        <  1:4462414/150‑1 (MQ=255)
               agagagAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGcac                                                                                                                                       >  2:3457281/1‑150 (MQ=255)
               agagagAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTAATCTTCAGTACGcac                                                                                                                                       <  2:2662381/150‑1 (MQ=255)
                gagagaAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGcacc                                                                                                                                      <  1:2471112/150‑1 (MQ=255)
                   agaAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACt                                                                                                                                   >  1:3651738/1‑150 (MQ=255)
                    gaAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACtt                                                                                                                                  <  2:3824945/150‑1 (MQ=255)
                     aaaaaaCAAGGCCCACCGGAACGGCAGGCCTGAGAATTCCGGCTGATAATAACCCACGCCACGGTCGTTTTCTGTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACttt                                                                                                                                 <  2:218169/150‑1 (MQ=255)
                      aaaaaCAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTAATCTTCAGTACGCACCACTTTa                                                                                                                                >  1:4894783/1‑150 (MQ=255)
                       aaaaCAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGTCCCATTTCATCTTCAGTACGCACCACTTTAc                                                                                                                               >  1:5135542/1‑150 (MQ=255)
                         aaCAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACcg                                                                                                                             >  2:4571266/1‑150 (MQ=255)
                         aaCAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTAATCTTCAGTACGCACCACTTTACcg                                                                                                                             <  1:2172935/150‑1 (MQ=255)
                             aGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTAATCTTCAGTACGCACCACTTTACCGCGCa                                                                                                                         <  2:3496318/150‑1 (MQ=255)
                              ggCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTAATCTTCAGTACGCACCACTTTACCGCGCag                                                                                                                        <  1:5085663/150‑1 (MQ=255)
                                cccACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCagag                                                                                                                      >  1:2391192/1‑150 (MQ=255)
                                cccACCGGAACGGCAGGCCTGAGAATAACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTATGCGGGTACGGGCAAACACTGATTCCGGTGATTCTGCCACGCGCAGCCCCATTTCATCAACAGTACGCACCACATTACCGCGCagag                                                                                                                      >  1:912822/1‑150 (MQ=255)
                                 ccACCGGCACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCagaga                                                                                                                     <  2:4186060/150‑1 (MQ=255)
                                  cACCGGAACGGCAGGCCTGAGAATTACGGCTGATAAAAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCagagag                                                                                                                    >  2:2233516/1‑150 (MQ=255)
                                   aCCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGt                                                                                                                   >  1:1697452/1‑150 (MQ=255)
                                   aCCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTAATCTTCAGTACGCACCACTTTACCGCGCAGAGAGt                                                                                                                   <  2:423480/150‑1 (MQ=255)
                                     cGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTc                                                                                                                 >  2:1366799/1‑150 (MQ=255)
                                     cGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTc                                                                                                                 <  2:176463/150‑1 (MQ=255)
                                      ggAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTAATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCg                                                                                                                >  1:105251/1‑150 (MQ=255)
                                       gAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCgg                                                                                                               <  2:3987084/150‑1 (MQ=255)
                                       gAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTAATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCgg                                                                                                               >  1:3019216/1‑150 (MQ=255)
                                          cGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTa                                                                                                            <  2:5108079/150‑1 (MQ=255)
                                           ggCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAg                                                                                                           <  2:1022801/150‑1 (MQ=255)
                                             cAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGAc                                                                                                         >  2:2707294/1‑150 (MQ=255)
                                             cAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGAc                                                                                                         >  2:2566268/1‑150 (MQ=255)
                                             cAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGAc                                                                                                         <  2:335257/150‑1 (MQ=255)
                                             cAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGAc                                                                                                         <  2:3897542/150‑1 (MQ=255)
                                             cAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTAATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGAc                                                                                                         <  1:3049415/150‑1 (MQ=255)
                                              aGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACg                                                                                                        <  1:4928121/150‑1 (MQ=255)
                                                    gagaATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTa                                                                                                  <  2:777250/150‑1 (MQ=255)
                                                      gaATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAAt                                                                                                <  1:101467/150‑1 (MQ=255)
                                                      gaATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAAt                                                                                                <  2:2880668/150‑1 (MQ=255)
                                                          tACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCt                                                                                            >  2:1721587/1‑150 (MQ=255)
                                                           aCGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTa                                                                                           >  1:1357190/1‑150 (MQ=255)
                                                             ggCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACa                                                                                         >  1:3420088/1‑150 (MQ=255)
                                                                tGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATcg                                                                                      >  2:2938738/1‑150 (MQ=255)
                                                                   taataaCCCACGCCTAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTTATTTCTACATcgacg                                                                                   <  1:2926128/150‑1 (MQ=255)
                                                                   taataaCCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATcgacg                                                                                   <  1:442684/150‑1 (MQ=255)
                                                                    aataaCCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATcgacga                                                                                  <  1:93499/150‑1 (MQ=255)
                                                                    aataaCCCACGCCAAGGTCGTTTTATTTTCGGGTACGGGCAAACACTGATTCCGGTGTTTAAGCCACGCGCAGACCCATTTCATATTAAGTACGCACCACTATAACGCGAAGAGATATCGGGTAGACGACGGTACATTCTACATcgacga                                                                                  >  1:5028331/1‑150 (MQ=255)
                                                                       aaCCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAAtt                                                                               >  1:2630006/1‑150 (MQ=255)
                                                                       aaCCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAAtt                                                                               <  1:4708712/150‑1 (MQ=255)
                                                                        aCCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAAttt                                                                              <  1:2939007/150‑1 (MQ=255)
                                                                          ccACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTAc                                                                            >  2:1570172/1‑150 (MQ=255)
                                                                             cGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTAAGTACGCACCACTTTAACGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCAACATCGACGAATTTAACGa                                                                         >  1:2707924/1‑150 (MQ=255)
                                                                                 aaGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCat                                                                     >  2:1202679/1‑150 (MQ=255)
                                                                                  aGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCata                                                                    >  1:1669347/1‑150 (MQ=255)
                                                                                   ggTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACAACGACGAATTTACCGAACatat                                                                   >  1:3216619/1‑150 (MQ=255)
                                                                                   ggTCGTTTTATTTTCGGGTACGGGCAATCACATATTACGGTGTTTATTCCACGCGCAGACCCAATTAATCTTCAGTACGCACCACATTACCGCCCAGAGAGTTCGGGTAGACGTCGGTAATTTATACAACGACGAAATTAACGATCatat                                                                   >  1:1504980/1‑150 (MQ=255)
                                                                                     tCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATcc                                                                 <  1:2972866/150‑1 (MQ=255)
                                                                                     tCGTTTTCTTTGCGGGTACGGGCAAACACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATcc                                                                 <  2:2026314/150‑1 (MQ=255)
                                                                                       gTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCgg                                                               >  1:1145878/1‑150 (MQ=255)
                                                                                       gTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCgg                                                               <  2:2264496/150‑1 (MQ=255)
                                                                                            cTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGc                                                          >  1:3946745/1‑150 (MQ=255)
                                                                                             tttGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACAACGACGAATTTACCGATCATATCCGGCGTGcc                                                         <  2:3379373/150‑1 (MQ=255)
                                                                                               tGCGGGTACGGGCACTCACAGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCtt                                                       <  2:2620409/150‑1 (MQ=255)
                                                                                               tGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCtt                                                       <  2:4648451/150‑1 (MQ=255)
                                                                                                gCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTc                                                      <  1:22642/150‑1 (MQ=255)
                                                                                                 cGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCg                                                     <  1:4202900/150‑1 (MQ=255)
                                                                                                   ggTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGaa                                                   >  1:4755865/1‑150 (MQ=255)
                                                                                                    gTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAg                                                  >  2:1850391/1‑150 (MQ=255)
                                                                                                        gggCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGa                                              >  1:1726802/1‑150 (MQ=255)
                                                                                                         ggCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGAc                                             <  2:1974076/150‑1 (MQ=255)
                                                                                                          gCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGAcc                                            <  2:477582/150‑1 (MQ=255)
                                                                                                           cAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCa                                           >  2:691631/1‑150 (MQ=255)
                                                                                                           cAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCa                                           <  2:4596566/150‑1 (MQ=255)
                                                                                                             aTCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTAATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCAcg                                         >  1:3520063/1‑150 (MQ=255)
                                                                                                              tCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCAcgc                                        >  1:4270021/1‑150 (MQ=255)
                                                                                                                  tGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGtt                                    >  1:668635/1‑150 (MQ=255)
                                                                                                                  tGATGCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGCGTTCGGGTAGACGTCGGTAATTTCTACATCGCCGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCCCGCGGtt                                    <  2:3167300/150‑1 (MQ=255)
                                                                                                                    aTTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTAt                                  <  1:4148388/150‑1 (MQ=255)
                                                                                                                      tCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTAttt                                <  1:2039172/150‑1 (MQ=255)
                                                                                                                             tttCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTAc                         >  2:3017592/1‑150 (MQ=255)
                                                                                                                              ttCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTACg                        <  2:4541331/150‑1 (MQ=255)
                                                                                                                              ttCAGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTACg                        <  2:4963271/150‑1 (MQ=255)
                                                                                                                                 tGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTACGccc                     <  1:3287739/150‑1 (MQ=255)
                                                                                                                                 tGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTACGccc                     <  1:711817/150‑1 (MQ=255)
                                                                                                                                     aCGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTACGCCCGGaa                 <  1:1316197/150‑1 (MQ=255)
                                                                                                                                     aCGCGCAGACCCATTTAATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTACGCCCGGaa                 <  2:788863/150‑1 (MQ=255)
                                                                                                                                        cgcAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTACGCCCGGAAAGc              <  2:677074/150‑1 (MQ=255)
                                                                                                                                         gcAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTACGCCCGGAAAGCt             >  2:4884417/1‑150 (MQ=255)
                                                                                                                                         gcAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTACGCCCGGAAAGCt             <  2:3589082/150‑1 (MQ=255)
                                                                                                                                                    tCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTACGCCCGGAAAGCTCCATGATGctc  >  1:1852574/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TAAAGATGGGGATAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTACGCCCGGAAAGCTCCATGATGCTC  >  CP009273/688930‑689227

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: ‑