Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,980 A→G 6.7% intergenic (‑599/+32) trmF ← / ← trmE tRNA uridine 5‑carboxymethylaminomethyl modification protein/GTPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,9800AG 6.7% 219.9 / 4.6 77intergenic (‑599/+32)trmF/trmEtRNA uridine 5‑carboxymethylaminomethyl modification protein/GTPase
Reads supporting (aligned to +/- strand):  ref base A (33/38);  new base G (3/2);  total (37/40)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.44e-01

GCGAAAGGATGGGTACATCCACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGCGTAAGTCCTCAGCCAG  >  NC_002947/4840‑5122
                                                                                                                                            |                                                                                                                                              
gcgAAAGGATGGGTACATCCACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAgcggcg                                                                                                                                       >  1:1137890/1‑150 (MQ=255)
 cgAAAGGATGGGTACATCCACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCTTAGGGAgcggctt                                                                                                                                      <  1:1055601/150‑3 (MQ=255)
    aaGGATGGGTACATCCACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTg                                                                                                                                   <  2:1063305/150‑1 (MQ=255)
      ggATGGGTACATCCACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGc                                                                                                                                 >  1:133091/1‑150 (MQ=255)
           ggTACATACACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGc                                                                                                                            <  2:460624/150‑1 (MQ=255)
            gTACATCCACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCtt                                                                                                                              <  2:322753/147‑1 (MQ=255)
            gTACATCCACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCtt                                                                                                                              >  1:322753/1‑147 (MQ=255)
            gTACATCCACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCg                                                                                                                           >  2:604510/1‑150 (MQ=255)
             tACATCCACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGt                                                                                                                          <  2:133091/150‑1 (MQ=255)
                    acaGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCtgtg                                                                                                                   >  2:622335/1‑150 (MQ=255)
                    acaGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCtgtg                                                                                                                   <  2:118566/150‑1 (MQ=255)
                          aTCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACt                                                                                                            <  2:859549/150‑1 (MQ=255)
                             tGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCtgt                                                                                                                    <  2:1152360/140‑1 (MQ=255)
                             tGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCtgt                                                                                                                    >  1:1152360/1‑140 (MQ=255)
                                   ttGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATg                                                                                                    <  1:627118/150‑1 (MQ=255)
                                   ttGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATg                                                                                                    >  2:627118/1‑150 (MQ=255)
                                         gagCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAg                                                                                              >  1:821925/1‑150 (MQ=255)
                                          agCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAg                                                                                                                                               <  1:746875/100‑1 (MQ=255)
                                          agCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAg                                                                                                                                               >  2:746875/1‑100 (MQ=255)
                                            cTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTg                                                                                           <  2:821925/150‑1 (MQ=255)
                                              tcatcaCGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGACGCTCCCTTTGGAGCGACGCTTGGCCGTTCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGga                                                                                         >  1:389676/1‑150 (MQ=255)
                                                atcaCGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGgaga                                                                                       <  1:487267/150‑1 (MQ=255)
                                                 tcaCGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAg                                                                                                                                               >  1:765337/1‑93 (MQ=255)
                                                 tcaCGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAg                                                                                                                                               <  2:765337/93‑1 (MQ=255)
                                                   aCGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTgg                                                                                          <  2:146754/144‑1 (MQ=255)
                                                   aCGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTgg                                                                                          >  1:146754/1‑144 (MQ=255)
                                                   aCGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGaaaa                                                                                    <  1:622335/150‑1 (MQ=255)
                                                   aCGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGaaaa                                                                                    >  2:533995/1‑150 (MQ=255)
                                                   aCGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGaaaa                                                                                    >  2:1165673/1‑150 (MQ=255)
                                                     gggTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATc                                                                                  >  2:742921/1‑150 (MQ=255)
                                                     gggTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATc                                                                                  >  1:151155/1‑150 (MQ=255)
                                                       gTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGc                                                                                                         <  2:23953/125‑1 (MQ=255)
                                                       gTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGc                                                                                                         >  1:23953/1‑125 (MQ=255)
                                                       gTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCg                                                                                >  1:358851/1‑150 (MQ=255)
                                                           aCCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCa                                                                            >  1:402501/1‑150 (MQ=255)
                                                            cccGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCaa                                                                           >  1:301742/1‑150 (MQ=255)
                                                             ccGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAg                                                                          >  1:218751/1‑150 (MQ=255)
                                                                  cccACAGGTGCTTTGGTCAGCAGCCTCCATTCACCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGt                                                                     <  1:583996/150‑1 (MQ=255)
                                                                   caacaGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTc                                                                    <  2:472348/148‑1 (MQ=255)
                                                                       aGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATcc                                                                <  1:604510/150‑1 (MQ=255)
                                                                         gTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGa                                                              <  1:184021/150‑1 (MQ=255)
                                                                            attGGTTAAGCGGCCTCCATTCACCTCCTGCAACGAATTTCAGGCAAAAAAAAGCCTCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCt                                                           <  2:521847/149‑1 (MQ=255)
                                                                             ttGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTg                                                          >  1:306172/1‑150 (MQ=255)
                                                                             ttGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTg                                                          >  2:371157/1‑150 (MQ=255)
                                                                             tcgggTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTg                                                          >  2:1244778/3‑150 (MQ=255)
                                                                                   cAGCGGCCTGCATTCATCTTCTGCAACGAATTAAAGGAAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCt                                                    <  2:953799/150‑1 (MQ=255)
                                                                                       ggCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCtg                                                <  2:151155/150‑1 (MQ=255)
                                                                                        gCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCtgt                                               <  2:218751/150‑1 (MQ=255)
                                                                                          cTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGa                                             >  2:1245365/1‑150 (MQ=255)
                                                                                          cTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGa                                             <  1:371157/150‑1 (MQ=255)
                                                                                           tGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGc                                                                                                         <  2:786944/89‑1 (MQ=255)
                                                                                           tGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGc                                                                                                         >  1:786944/1‑89 (MQ=255)
                                                                                             cATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGAtct                                          >  1:515206/1‑150 (MQ=255)
                                                                                                 cATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTcccc                                      <  2:515206/150‑1 (MQ=255)
                                                                                                   tcttctGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAGTGCTCCTGTGATCTCCCCCa                                    <  1:533995/150‑1 (MQ=255)
                                                                                                   tcttctGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTATGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCa                                    >  2:627808/1‑150 (MQ=255)
                                                                                                     ttctGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAgc                                  >  2:91426/1‑150 (MQ=255)
                                                                                                              gAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCtgt                         >  1:931848/1‑150 (MQ=255)
                                                                                                                aTTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGc                       <  1:742921/150‑1 (MQ=255)
                                                                                                                  ttCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCtgtg                                                                                                                   >  2:80520/1‑56 (MQ=255)
                                                                                                                  ttCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCtgtg                                                                                                                   <  1:80520/56‑1 (MQ=255)
                                                                                                                   tCAGGGAAAAAAAAGCCGCTCCCTAAGTAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTg                    <  1:802803/150‑1 (MQ=255)
                                                                                                                     aGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGc                  <  2:358851/150‑1 (MQ=255)
                                                                                                                     aGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGc                  <  2:931848/150‑1 (MQ=255)
                                                                                                                         aaaaaaaaGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGCGTaa              <  1:553530/150‑1 (MQ=255)
                                                                                                                         aaaaaaaaGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGCGTaa              >  2:553530/1‑150 (MQ=255)
                                                                                                                               aaGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGCGTAAGTCCTc        <  2:389676/150‑1 (MQ=255)
                                                                                                                                aGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCt                           <  2:206603/130‑1 (MQ=255)
                                                                                                                                aGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCt                           >  1:206603/1‑130 (MQ=255)
                                                                                                                                aGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGCGTAAGTCCTCa       <  2:301742/150‑1 (MQ=255)
                                                                                                                                 gCCGCTCCTTAGGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTg                    >  1:1122220/1‑136 (MQ=255)
                                                                                                                                 gCCGCTCCTTAGGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTg                    <  2:1122220/136‑1 (MQ=255)
                                                                                                                                  ccGCTCCTTAGGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGCGTAAGTCCTCAGc     >  2:782901/1‑150 (MQ=255)
                                                                                                                                  ccGCTCCTTAGGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGCGTAAGTCCTCAGc     >  2:378191/1‑150 (MQ=255)
                                                                                                                                  ccGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGCGTAAGTCCTCAGc     >  1:372967/1‑150 (MQ=255)
                                                                                                                                    gCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGCGTAAGTCCTCAGCCa   <  1:300949/150‑1 (MQ=255)
                                                                                                                                     cTCCCTAAGGCGCGGCGCTCGGCCTTGCGTCGCTCTGCATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGCGTAAGTCCTCAGCCAg  <  2:372967/150‑1 (MQ=255)
                                                                                                                                           aaGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTc                                                   <  1:963424/95‑1 (MQ=255)
                                                                                                                                           aaGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTc                                                   >  2:963424/1‑95 (MQ=255)
                                                                                                                                            |                                                                                                                                              
GCGAAAGGATGGGTACATCCACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGCGTAAGTCCTCAGCCAG  >  NC_002947/4840‑5122

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