Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,064,215 C→T 6.0% M261I (ATG→ATA PP_1841 ← cytochrome c family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,064,2150CT 6.0% 200.8 / 2.4 67M261I (ATG→ATAPP_1841cytochrome c family protein
Reads supporting (aligned to +/- strand):  ref base C (29/34);  new base T (2/2);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

AGCAAGAGCGGCTCACGGCCTGACAGTTTTGTGTTGGCCCGTTCGCGGGTTAACCCGCGCCTACACATTGATGGATGGACACCACCTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCGGCAAACCCCGGCATCTTCGCCTGCGGCCACCGCCGCAGGCTTTGCGGGTCGCGTATCAACTTGCGTAAATAAACCGGCTGAAAATA  >  NC_002947/2064074‑2064358
                                                                                                                                             |                                                                                                                                                
agcaagAGCGGCTCACGGCCTGACAGTTTTGTGTTGGCCCGTTCGCGGGTTAACCCGCGCCTACACATTGATGGATGGACACCACCTGCCGTTGTGGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCTATGTGcc                                                                                                                                           >  5:334327/1‑149 (MQ=255)
agcaagAGCGGCTCACGGCCTGACAGTTTTGTGTTGGCCCGTTCGCGGGTTAACCCGCGCCTACACATTGATGGATGGACACCACCTGCCGTTGTGGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCTATGTGcc                                                                                                                                           >  5:329890/1‑149 (MQ=255)
agcaagAGCGGCTCACGGCCTGACAGTTTTGTGTTGGCCCGTTCGCGGGTTAACCCGCGCCTACACATTGATGGATGGACACCACCTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGcc                                                                                                                                           >  8:190450/1‑149 (MQ=255)
    agagCGGCTCACGGCCTGACAGTTTTGTGTTGGCCCGTTCGCGGGTTAAACCGCGCCTACACATTGATGGATGGACACCACCTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCa                                                                                                                                        >  3:231820/1‑148 (MQ=255)
     gagCGGCTCACGGCCTGACAGTTTTGTGTTGGCCCGTTCGCGGGTTAACCCGCGCCTACACATTGATGGATGGACACCACCTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCaaa                                                                                                                                      <  5:131887/149‑1 (MQ=255)
       gCGGCTCACGGCCTGACAGTTTTGTGTTGGCCCGTTCGCGGGTTAACCCGCGCCTACACATTGATGGATGGACACCACCTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATa                                                                                                                                    <  5:211041/149‑1 (MQ=255)
        cGGCTCACGGCCTGACAGTTTTGTGTTGGCCCGTTCGCGGGTTAACCCGCGCCTACACATTGATGGATGGACACCACCTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGTCCCAAATAg                                                                                                                                   <  5:123372/149‑1 (MQ=255)
         ggCTCACGGCCTGACAGTTTTGTGTTGGCCCGTTCGCGGGTTAACCCGCGCCTACACATTGATGGATGGACACCACCTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAgg                                                                                                                                  <  2:258006/149‑1 (MQ=255)
           cTCACGGCCTGACAGTTTTGTGTTGGCCCGTTCGCGGGTTAACCCGCGCCTACACATTGATGGATGGACACCACCTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGc                                                                                                                                 >  8:43155/1‑148 (MQ=255)
                 gCCTGACAGTTTTGTGTTGGCCCGTTCGCGGGTTAACCCGCGCCTACACATTGATGGATGGACACCACCTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAgc                                                                                                                          <  7:43155/149‑1 (MQ=255)
                    tGACAGTTTTGTGTTGGCCCGTTCGCGGGTTAACCCGCGCCTACACATTGATGGATGGACACCACCTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCa                                                                                                                       >  1:46490/1‑149 (MQ=255)
                       cAGTTTTGTGTTGGCCCGTTCGCGGGTTAACCCGCGCCTACACATTGATGGATGGACACCACCTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATc                                                                                                                     <  3:180644/148‑1 (MQ=255)
                          ttttGTGTTGGCCCGTTCGCGGGTTAACCCGCGCCTACACATTGATGGATGGACACCACCTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAgg                                                                                                                                  <  4:206386/132‑1 (MQ=255)
                          ttttGTGTTGGCCCGTTCGCGGGTTAACCCGCGCCTACACATTGATGGATGGACACCACCTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAgg                                                                                                                                  >  3:206386/1‑132 (MQ=255)
                          ttttGTGTTGGCCCGTTCGCGGGTTAACCCGCGCCTACACATTGATGGATGGACACCACCTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATc                                                                                                                     <  6:122103/145‑1 (MQ=255)
                          ttttGTGTTGGCCCGTTCGCGGGTTAACCCGCGCCTACACATTGATGGATGGACACCACCTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATc                                                                                                                     >  5:122103/1‑145 (MQ=255)
                                gtTGGCCCGTTCGCGGGTTAACCCGCGCCTACACATTGATGGATGGACACCACCTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCt                                                                                                             <  6:77111/147‑1 (MQ=255)
                                gtTGGCCCGTTCGCGGGTTAACCCGCGCCTACACATTGATGGATGGACACCACCTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCt                                                                                                             >  5:77111/1‑147 (MQ=255)
                                 ttGGCCCGTTCGCGGGTTAACCCGCGCCTACACATTGATGGATGGACACCACCTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGtt                                                                                                          <  4:278480/149‑1 (MQ=255)
                                    gCCCGTTCGCGGGTTAACCCGCGCCTACACATTGATGGATGGACACCACCTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCAc                                                                                                       >  3:22052/1‑149 (MQ=255)
                                        gTTCGCGGGTTAACCCGCGCCTACACATTGATGGATGGACACCACCTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCaa                                                                                                  <  3:19141/149‑1 (MQ=255)
                                                 ttAACCCGCGCCTACACATTGATGGATGGACACCACCTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTc                                                                                         <  2:46490/149‑1 (MQ=255)
                                                      ccGCGCCTACACATTGATGGATGGACACCACCTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCGGCa                                                                                     >  2:277662/1‑148 (MQ=255)
                                                      ccGCGCCTACACATTGATGGATGGACACCACCTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCGGCa                                                                                     >  2:127562/1‑148 (MQ=255)
                                                      ccGCGCCTACACATTGATGGATGGACACCACCTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCGGCa                                                                                     >  1:302008/1‑148 (MQ=255)
                                                       cgggcCAACACATTGATGGATGGACACCACCTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCGGCaaa                                                                                   <  4:37837/146‑1 (MQ=255)
                                                       cgcgCCTACACATTGATGGATGGACACCACCTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCGGCaaa                                                                                   <  1:127562/149‑1 (MQ=255)
                                                            cTACACATTGATGGATGGACACCACCTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTc                                                                                         <  2:330810/138‑1 (MQ=255)
                                                            cTACACATTGATGGATGGACACCACCTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTc                                                                                         >  1:330810/1‑138 (MQ=255)
                                                                  aTTGATGGATGGACACCACCTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGtt                                                                                                          <  3:283065/116‑1 (MQ=255)
                                                                  aTTGATGGATGGACACCACCTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGtt                                                                                                          >  4:283065/1‑116 (MQ=255)
                                                                  aTTGATGGATGGACACCACCTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCGGCAAACCCCGGCATCt                                                                        >  6:317385/1‑149 (MQ=255)
                                                                     gatggatgGACACCACCTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCGGCAAACCCCGGCATCTTc                                                                      <  7:190450/148‑1 (MQ=255)
                                                                     gatggaggGACACCACCTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTGGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCGGCAAACCCCGGCATCTTc                                                                      <  8:292181/148‑1 (MQ=255)
                                                                      atggatggACACCACCTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCGGCAAACCCCGGCATCTTCGc                                                                    >  3:127954/1‑149 (MQ=255)
                                                                               caccacCTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCGGCAAACCCCGGCATCTTCGCCTGCGGCCa                                                           <  4:22052/149‑1 (MQ=255)
                                                                               caccacCTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCGGCAAACCCCGGCATCTTCGCCTGCGGCCa                                                           >  6:12195/1‑149 (MQ=255)
                                                                                accaccTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCGGCAAACCCCg                                                                              >  7:11415/1‑129 (MQ=255)
                                                                                accaccTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCGGCAAACCCCg                                                                              <  8:11415/129‑1 (MQ=255)
                                                                                  caccTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCGGCAAACCCCGGCATCTTCGCCTGCGGCCAccg                                                        >  6:11187/1‑149 (MQ=255)
                                                                                    ccTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCGGCAAACCCAGGCATCTTCGCCTGCGGCCAccgcc                                                      >  5:189877/1‑149 (MQ=255)
                                                                                       gCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCGGCAAACCCCGGCATCTTCGCCTGCGGCCAccgccg                                                     <  5:266720/147‑1 (MQ=255)
                                                                                       gCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCGGCAAACCCCGGCATCTTCGCCTGCGGCCAccgccg                                                     >  6:266720/1‑147 (MQ=255)
                                                                                        ccGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCGGCAAACCCCGGCATCTTCGCCTGCGGCCACCGCCGCa                                                   <  4:127954/148‑1 (MQ=255)
                                                                                          gTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCGGCAAACCCCGGCATCTTCGCCTGCGGCGGTGGCCGCAGGc                                                <  5:205303/149‑1 (MQ=255)
                                                                                           ttGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGcc                                                                                                                                >  5:60421/1‑69 (MQ=255)
                                                                                           ttGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGcc                                                                                                                                <  6:60421/69‑1 (MQ=255)
                                                                                                 gAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCGGCAAACCCCGGCATCTTCGCCTGCGGCCACCGCCGCAGGCTTTGCgg                                         >  7:55373/1‑149 (MQ=255)
                                                                                                      ggTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCGGCAAACCCCGGCATCTTCGCCTGCGGCCACCGCCGCAGGCTTTGCGGGTcgc                                    <  6:192492/149‑1 (MQ=255)
                                                                                                        ttCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCGGCAAACCCCGGCATCTTCGCCTGCGGCCACCGCCGCAGGCTTTGCGGGTcgcg                                   >  5:115802/1‑148 (MQ=255)
                                                                                                             ccGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACt                                                                                                      >  2:230806/1‑77 (MQ=255)
                                                                                                             ccGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACt                                                                                                      <  1:230806/77‑1 (MQ=255)
                                                                                                             ccGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCg                                                                                        >  7:171666/1‑90 (MQ=255)
                                                                                                             ccGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCg                                                                                        <  8:171666/90‑1 (MQ=255)
                                                                                                             ccGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCGGCAAACCCCGGCATCTTCGCCTGCGGCCACCGCCGCAGGCTTTGCGGGTCGCGTa                                 >  2:147919/1‑145 (MQ=255)
                                                                                                             ccGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCGGCAAACCCCGGCATCTTCGCCTGCGGCCACCGCCGCAGGCTTTGCGGGTCGCGTa                                 <  1:147919/145‑1 (MQ=255)
                                                                                                             ccGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCGGCAAACCCCGGCATCTTCGCCTGCGGCCACCGCCGCAGGCTTTGCGGGTCGCGTATCaa                             >  5:102670/1‑149 (MQ=255)
                                                                                                               gcgAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCGGCAAACCCCGGCATCTTCGCCTGCGGCCACCGCCGCAGGCTTTGCGGGTCGCGTATCAACt                           >  3:5751/1‑149 (MQ=255)
                                                                                                                   aTAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCGGCAAACCCCGGCATCTTCGCCTGCGGCCACCGCCGCAGGCTTTGCGGGTCGCGTATCAACTTGc                        <  5:12195/148‑1 (MQ=255)
                                                                                                                    tAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCGGCAAACCCCGGCATCTTCGCCTGCGGCCACCGCCGCAGGCTTTGCGGGTCGCGTATCAACTTGCGt                      <  2:302008/149‑1 (MQ=255)
                                                                                                                       gCCCTTAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCGGCAAACCCCGGCATCTTCGCCTGCGGCCACCGCCGCAGGCTTTGCGGGTCGCGTATCAACTTGCGtaaa                   <  1:277662/149‑1 (MQ=255)
                                                                                                                             aaGGCTTGCGCGTTGCTATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCGGCAAACCCCGGCATCTTCGCCTGCGGCCACCGCCGCAGGCTTTGCGGGTCGCGTATCAACTTGCGtaaataaa               <  6:334327/147‑1 (MQ=255)
                                                                                                                             aaGGCTTGCGCGTTGCTATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCGGCAAACCCCGGCATCTTCGCCTGCGGCCACCGCCGCAGGCTTTGCGGGTCGCGTATCAACTTGCGtaaataaa               <  6:329890/147‑1 (MQ=255)
                                                                                                                             aaGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCGGCAAACCCCGGCTTCTTCGCCTGCGGCCACCGCCGCAGGCTTTGCGGGTCGCGTATCAACTTGCGTAAATAAAc              <  8:55373/148‑1 (MQ=255)
                                                                                                                             aaGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCGGCAAACCCCGGCATCTTCGCCTGCGGCCACCGCCGCAGGCTTTGCGGGTCGCGTATCAACTTGCGtaaataaa               <  6:115802/147‑1 (MQ=255)
                                                                                                                                    gcgcGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAGGCACACTTACGGCAAACCCCGCCATATTCGCATGCGGCCACCCCCGCAGGATTTGCGGGTCGGTTATCAACTTGGGTAAATAAACCGGATg        >  1:255428/1‑148 (MQ=255)
                                                                                                                                        gTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCGGCAAACCCCGGCATCTTCGCCTGCGGCCACCGCCGCAGGCTTTGCGGGTCGCGTATCAACTTGCGTAAATAAACCGGCTGAAAATa  <  4:5751/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                
AGCAAGAGCGGCTCACGGCCTGACAGTTTTGTGTTGGCCCGTTCGCGGGTTAACCCGCGCCTACACATTGATGGATGGACACCACCTGCCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCCTCAAGGCTTGCGCGTTGCCATGTGCCCCAAATAGGCCAACAGCGCATCCAGTTCCTGTTCACTCAACACACTTTCGGCAAACCCCGGCATCTTCGCCTGCGGCCACCGCCGCAGGCTTTGCGGGTCGCGTATCAACTTGCGTAAATAAACCGGCTGAAAATA  >  NC_002947/2064074‑2064358

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: