Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,977 C→T 8.9% intergenic (‑596/+35) 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,9770CT 8.9% 190.8 / 7.3 69intergenic (‑596/+35)trmF/trmEtRNA uridine 5‑carboxymethylaminomethyl modification protein/GTPase
Reads supporting (aligned to +/- strand):  ref base C (30/32);  new base T (4/2);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGAGCCCGCGAAAGGATGGGTACATCCACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGCGTAAGTCCTCAGCCAG  >  NC_002947/4832‑5122
                                                                                                                                                  |                                                                                                                                                  
tCGAGCCCGCGAAAGGATGGGTACATCCACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTaa                                                                                                                                                 <  2:108345/149‑1 (MQ=255)
      ccGCGAAAGGATGGGTACATCCACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCTTAGGGAgcg                                                                                                                                           >  5:99247/1‑149 (MQ=255)
        gcgAAAGGATGGGTACATCCACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAgcggc                                                                                                                                         >  3:264110/1‑149 (MQ=255)
          gAAAGGATGGGTACATCCACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCTTAGGGAgcggctt                                                                                                                                       <  8:194250/149‑3 (MQ=255)
          gAAAGGATGGGTACATCCACAGAGATCTGGTGCTTGTAACAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGgcgc                                                                                                                                       <  8:47997/149‑1 (MQ=255)
           aaaGGATGGGTACATCCACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCt                                                                                                                                      <  4:264110/149‑1 (MQ=255)
                aTGGGTACATCCACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGcc                                                                                                                                 >  2:323802/1‑149 (MQ=255)
                 tGGGTACATCCACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAg                                                                                                                                             <  4:106190/136‑1 (MQ=255)
                 tGGGTACATCCACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAg                                                                                                                                             >  3:106190/1‑136 (MQ=255)
                    gTACATCCACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGc                                                                                                                            >  3:203951/1‑149 (MQ=255)
                        aTCCACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCGGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTgg                                                                                                                        >  6:265878/1‑149 (MQ=255)
                         tCCACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGc                                                                                                                       >  2:22700/1‑149 (MQ=255)
                          ccACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCt                                                                                                                      <  3:150400/149‑1 (MQ=255)
                          ccACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCt                                                                                                                      >  6:323380/1‑149 (MQ=255)
                               gagaTCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGa                                                                                                                 <  4:228407/149‑1 (MQ=255)
                               gagaTCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGa                                                                                                                 <  4:194630/149‑1 (MQ=255)
                                agaTCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGAt                                                                                                                >  5:129529/1‑149 (MQ=255)
                                     tGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACtt                                                                                                           >  6:298825/1‑149 (MQ=255)
                                      ggTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACtt                                                                                                           <  1:323802/148‑1 (MQ=255)
                                       gTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGc                                                                                                         <  4:244236/149‑1 (MQ=255)
                                       gTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGc                                                                                                         <  4:203951/149‑1 (MQ=255)
                                       gTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGc                                                                                                         >  1:179266/1‑149 (MQ=255)
                                         gCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCTTAGGGAgcggcttt                                                                                                                                      <  2:222619/119‑4 (MQ=255)
                                         gCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCTTAGGGAgcggcttt                                                                                                                                      >  1:222619/1‑116 (MQ=255)
                                          cTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCg                                                                                                       <  1:6181/148‑1 (MQ=255)
                                            tGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGa                                                                                                                                              <  5:164804/108‑1 (MQ=255)
                                            tGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGa                                                                                                                                              >  6:164804/1‑108 (MQ=255)
                                             gTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATg                                                                                                    >  3:88567/1‑148 (MQ=255)
                                              tAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGc                                                                                                   <  5:323380/148‑1 (MQ=255)
                                                  agCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAg                                                                                              >  5:103086/1‑149 (MQ=255)
                                                    cTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCtgtg                                                                                                                   >  1:305861/1‑126 (MQ=255)
                                                       catcaCGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTgg                                                                                          <  6:129529/148‑1 (MQ=255)
                                                            cgatgtcACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGaaaa                                                                                    <  3:118334/143‑1 (MQ=255)
                                                             gggTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGAGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGaaaa                                                                                    >  2:144484/1‑149 (MQ=255)
                                                             gggTCCACCCGCTCCCACAGGTGATTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCGTAAGGAGCGGCGCTTGGCCTGGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGCAAAt                                                                                   >  2:201866/1‑149 (MQ=255)
                                                               gTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAACTGGAGAAAATcc                                                                                 <  2:179266/149‑1 (MQ=255)
                                                                tCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCg                                                                                <  1:22700/149‑1 (MQ=255)
                                                                tCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCg                                                                                >  5:89580/1‑149 (MQ=255)
                                                                   aCCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGaaaa                                                                                    >  4:118334/1‑142 (MQ=255)
                                                                     ccGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGc                                                                               >  7:166391/1‑145 (MQ=255)
                                                                     ccGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGc                                                                               <  8:166391/145‑1 (MQ=255)
                                                                                 gTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGCAGAAAATCCGCCCAAGCAGGTCATcc                                                                >  6:107700/1‑148 (MQ=255)
                                                                                   gCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAg                                                             >  2:256881/1‑149 (MQ=255)
                                                                                    cTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAg                                                             <  6:103086/148‑1 (MQ=255)
                                                                                      tGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCt                                                           >  5:109909/1‑148 (MQ=255)
                                                                                         gTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAgg                                                                      <  6:18365/134‑1 (MQ=255)
                                                                                         gTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAgg                                                                      >  5:18365/1‑134 (MQ=255)
                                                                                                gCCAGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCtg                                                <  5:265878/149‑1 (MQ=255)
                                                                                                      aTTCATCTTCTGCAACGAATTTCAGTCAAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGAtc                                           <  1:144484/149‑1 (MQ=255)
                                                                                                      aTTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGa                                                              <  3:13877/129‑1 (MQ=255)
                                                                                                      aTTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGa                                                              >  4:13877/1‑129 (MQ=255)
                                                                                                         cATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTccc                                       >  6:34670/1‑149 (MQ=255)
                                                                                                          aTCTTCTGCAACGAATTTCAGGCAAAAAATAGCCGCTCCCTAAGGATCGGCGCTTGGAATTGCGTGGCTGTGGATTACTTGCCGATGGAGAAGCTGGAGAAAATCCGCCCATGAAGGTCATCCGAGCTGAATTCTCCTGTGATCTcccc                                      <  4:245937/149‑1 (MQ=255)
                                                                                                            cttctGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTcccc                                      >  3:315607/1‑147 (MQ=255)
                                                                                                               ctGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAgc                                  <  4:88567/148‑1 (MQ=255)
                                                                                                                  cAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAgcgc                                <  6:109909/147‑1 (MQ=255)
                                                                                                                     cGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGAGGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGc                            >  1:80447/1‑149 (MQ=255)
                                                                                                                       aaTTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCtgt                         <  5:298825/149‑1 (MQ=255)
                                                                                                                          ttCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGAtc                                           >  1:23093/1‑128 (MQ=255)
                                                                                                                          ttCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGAtc                                           <  2:23093/128‑1 (MQ=255)
                                                                                                                              ggCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGc                  >  8:60802/1‑149 (MQ=255)
                                                                                                                               gCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCg                                                                                                       >  3:175267/1‑63 (MQ=255)
                                                                                                                               gCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCg                                                                                                       <  4:175267/63‑1 (MQ=255)
                                                                                                                                   aaaaaaaGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGCGTaa              <  6:89580/149‑1 (MQ=255)
                                                                                                                                        gcgCCGCTCCTTAGGGAGCGGCGCTTGGCCTTGCGTGGGTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGCGTAAGTCCt         >  6:208186/3‑149 (MQ=255)
                                                                                                                                          gCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGCGTAAGTCCTCa       <  5:107700/149‑1 (MQ=255)
                                                                                                                                           ccGCTCCTTAGGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGCGTAAGTCCTCAg      >  6:250861/1‑149 (MQ=255)
                                                                                                                                               tCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTg                    >  7:215342/1‑131 (MQ=255)
                                                                                                                                               tCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTg                    <  8:215342/131‑1 (MQ=255)
                                                                                                                                                 ccTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGCGTAAGTCCTCAGCCAg  >  8:111267/1‑147 (MQ=255)
                                                                                                                                                  cTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGCGTAAGTCCTCAGCCAg  >  7:193252/1‑146 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TCGAGCCCGCGAAAGGATGGGTACATCCACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGCGTAAGTCCTCAGCCAG  >  NC_002947/4832‑5122

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