Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,980 A→G 8.4% 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 8.4% 202.8 / 6.9 71intergenic (‑599/+32)trmF/trmEtRNA uridine 5‑carboxymethylaminomethyl modification protein/GTPase
Reads supporting (aligned to +/- strand):  ref base A (34/31);  new base G (4/2);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

TCGAGCCCGCGAAAGGATGGGTACATCCACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGCGTAAGTCCTCAGCCAGAAG  >  NC_002947/4832‑5125
                                                                                                                                                     |                                                                                                                                                  
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                                                                                                                                                >  3:106190/1‑136 (MQ=255)
                 tGGGTACATCCACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAg                                                                                                                                                <  4:106190/136‑1 (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:203951/149‑1 (MQ=255)
                                       gTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGc                                                                                                            <  4:244236/149‑1 (MQ=255)
                                       gTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGc                                                                                                            >  1:179266/1‑149 (MQ=255)
                                         gCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCTTAGGGAgcggcttt                                                                                                                                         >  1:222619/1‑116 (MQ=255)
                                         gCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCTTAGGGAgcggcttt                                                                                                                                         <  2:222619/119‑4 (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                                                                                   >  5:89580/1‑149 (MQ=255)
                                                                tCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCg                                                                                   <  1:22700/149‑1 (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                                                                         >  5:18365/1‑134 (MQ=255)
                                                                                         gTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAgg                                                                         <  6:18365/134‑1 (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                                                                                                          <  4:175267/63‑1 (MQ=255)
                                                                                                                               gCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCg                                                                                                          >  3:175267/1‑63 (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)
                                                                                                                                                   tAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGCGTAAGTCCTCAGCCAGAAg  >  7:195363/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TCGAGCCCGCGAAAGGATGGGTACATCCACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGCGTAAGTCCTCAGCCAGAAG  >  NC_002947/4832‑5125

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: