Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,790,103 C→T 6.7% G63D (GGC→GAC)  PP_2442 ← DUF480 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,790,1030CT 6.7% 224.8 / 5.0 75G63D (GGC→GAC) PP_2442DUF480 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (33/37);  new base T (2/3);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

CCCCGCAGGAACATCAGGCCCATCAGCACCAACTGCGCTGGCACCAGCTCCAGTGCCTTGTCTACCCGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGTTGCAGGCCAGGACCAGGGCATTGAGGGTCAGCGGGTAGCTTTCCGGGCTGGTGGCCTGCTTTTCGATCAGCGAGCCGAGAACGCGGATCTCGATGC  >  NC_002947/2789959‑2790250
                                                                                                                                                |                                                                                                                                                   
ccccGCAGGAACATCAGGCCCATCAGCACCAACTGCGCTGGCACCAGCTCCAGTGCCTTGTCTACCCGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGAcc                                                                                                                                                <  2:216675/150‑1 (MQ=255)
  ccGCAGGAACATCAGGCCCATCAGCACCAACTGCGCTGGCACCAGCTCCAGTGCCTTGTCTACCCGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTg                                                                                                                                              >  1:844238/1‑150 (MQ=255)
   cGCAGGAACATCAGGCCCATCAGCACCAACTGCGCTGGCACCAGCTCCAGTGCCTTGTCTACCCGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTgg                                                                                                                                             >  2:1284238/1‑150 (MQ=255)
    gCAGGAACATCAGGCCCATCAGCACCAACTGCGCTGGCACCAGCTCCAGTGCCTTGTCTACCCGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGc                                                                                                                                            >  2:645235/1‑150 (MQ=255)
     cAGGAACATCAGGCCCATCAGCACCAACTGCGCTGGCACCAGCTCCAGTGCCTTGTCTACCCGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGcc                                                                                                                                           <  1:1284238/150‑1 (MQ=255)
       ggAACATCAGGCCCATCAGCACCAACTGCGCTGGCACCAGCTCCAGTGCCTTGTCTACCCGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCtt                                                                                                                                         <  1:736182/150‑1 (MQ=255)
         aaCATCAGGCCCATCAGCACCAACTGCGCTTGCACCAGCTCCAGTGCCTTGTCTACCCGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCGGCCCTTCGAGGGCGCGCGAGGCCTGGCCGACCTGGCCTTgg                                                                                                                                       <  2:559164/150‑1 (MQ=255)
         aaCATCAGGCCCATCAGCACCAACTGCGCTGGTACCAGCTCCAGTGCCTTGTCTACCCGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCAGGTCGGCCAGGCCTTgc                                                                                                                                       <  2:240402/150‑2 (MQ=255)
         aaCATCAGGCCCATCAGCACCAACTGCGCTGGCACCAGCTCCAGTGCCTTGTCTACCCGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTgg                                                                                                                                       <  2:850673/150‑1 (MQ=255)
           cATCAGGCCCATCAGCACCAACTGCGCTGGCACCAGCTCCAGTGCCTTGTCTACCCGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCCCAAGGCCAGGTCGGCCAGGCCTTgcgc                                                                                                                                     <  2:482487/150‑4 (MQ=255)
             tCAGGCCCATCAGCACCAACTGCGCTGGCACCAGCTCCAGTGCCTTGTCTACCCGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGa                                                                                                                                   <  1:251843/150‑1 (MQ=255)
               aGGCCCATCAGCACCAACTGCGCTGGCACCAGCTCCAGTGCCTTGTCTACCCGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAgg                                                                                                                                 <  1:869818/150‑1 (MQ=255)
               aGGCCCATCAGCACCAACTGCGCTGGCACCAGCTCCAGTGCCTTGTCTACCCGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAgg                                                                                                                                 >  2:319673/1‑150 (MQ=255)
                    cATCAGCACCAACTGCGCTGGCACCAGCTCCAGTGCCTTGTCTACCCGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTtcat                                                                                                                            >  2:1127316/1‑150 (MQ=255)
                             cAACTGCGCTGGCACCAGCTCCAGTGCCTTGTCTACCCGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTc                                                                                                                   >  1:1247181/1‑150 (MQ=255)
                                  gcgcTGGCACCAGCTCCAGTGCCTTGTCTACCCGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGc                                                                                                              >  2:795177/1‑150 (MQ=255)
                                  gcgcTGGCACCAGCTCCAGTGCCTTGTCTACCCGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGc                                                                                                              <  2:844238/150‑1 (MQ=255)
                                          aCCAGCTCCAGTGCCTTGTCTACCCGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTtca                                                                                                                             >  1:803564/1‑127 (MQ=255)
                                          aCCAGCTCCAGTGCCTTGTCTACCCGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTtca                                                                                                                             <  2:803564/127‑1 (MQ=255)
                                            cAGCTCCAGTGCCTTGTCTACCCGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCgg                                                                                                               <  2:1408854/139‑1 (MQ=255)
                                            cAGCTCCAGTGCCTTGTCTACCCGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCgg                                                                                                               >  1:1408854/1‑139 (MQ=255)
                                            cAGCTCCAGTGCCTTGTCTACCCGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTg                                                                                                    <  1:273771/150‑1 (MQ=255)
                                            cAGCTCCAGTGCCTTGTCTACCCGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTg                                                                                                    <  1:645235/150‑1 (MQ=255)
                                               cTCCAGTGCCTTGTCTACCCGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGtt                                                                                                 <  1:988676/150‑1 (MQ=255)
                                                 ccAGTGCCTTGTCTACCCGAGGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGTTGc                                                                                               <  1:555484/150‑1 (MQ=255)
                                                  cAGTGCCTTGTCTACCCGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTtca                                                                                                                             >  1:734880/1‑119 (MQ=255)
                                                  cAGTGCCTTGTCTACCCGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTtca                                                                                                                             <  2:734880/119‑1 (MQ=255)
                                                   aGTGCCTTGTCTACCCGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGTTGCAg                                                                                             <  1:1127316/150‑1 (MQ=255)
                                                          tGTCTACCCGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGTTGCAGGCCAGGa                                                                                      <  1:227751/150‑1 (MQ=255)
                                                            tCTACCCGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGTTGCAGGCCAGGAcc                                                                                    >  2:1394392/1‑150 (MQ=255)
                                                              tACCCGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGTTGCAGGCCAGGACCAg                                                                                  >  2:212421/1‑150 (MQ=255)
                                                               aCCCGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGTTGCAGGCCAGGACCAgg                                                                                 <  1:1008864/150‑1 (MQ=255)
                                                                cccGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGcc                                                                                                                                                    >  1:1387294/1‑82 (MQ=255)
                                                                cccGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGcc                                                                                                                                                    <  2:1387294/82‑1 (MQ=255)
                                                                cccGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCAc                                                                                                                         >  2:320341/1‑109 (MQ=255)
                                                                cccGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCAc                                                                                                                         <  1:320341/109‑1 (MQ=255)
                                                                cccGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTg                                                                                                            <  1:1394756/122‑1 (MQ=255)
                                                                cccGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTg                                                                                                            >  2:1394756/1‑122 (MQ=255)
                                                                 ccGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGtcttct                                                                                                     <  1:1275458/128‑1 (MQ=255)
                                                                 ccGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGtcttct                                                                                                     >  2:1275458/1‑128 (MQ=255)
                                                                  cGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCtt                                                                                                                                         <  1:1297440/91‑1 (MQ=255)
                                                                  cGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCtt                                                                                                                                         >  2:1297440/1‑91 (MQ=255)
                                                                        tCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTTAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGTTGCAGGCCAGGACCAGGGCATTGAgg                                                                        <  2:1052411/150‑1 (MQ=255)
                                                                        tCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGTTGCAGGCCAGGACCAGGGCATTGAgg                                                                        <  1:763051/150‑1 (MQ=255)
                                                                          ccAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCAGGTCGGCCAGGCCTTgc                                                                                                                                       <  1:1187893/85‑2 (MQ=37)
                                                                          ccAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCAGGTCGGCCAGGCCTTgc                                                                                                                                       >  2:1187893/1‑84 (MQ=37)
                                                                                   ggCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGTTGCAGGCCAGGACCAGGGCATTGAGGGTCAGCGGGTa                                                             >  1:625735/1‑150 (MQ=255)
                                                                                    gCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGTTGCAGGCCAGGACCAGGGCATTGAGGGTCAGCGGGTAg                                                            >  1:421943/1‑150 (MQ=255)
                                                                                             cccATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGTTGCAGGCCAGGACCAGGGCATTGAGGGTCAGCGGGTAGCTTTCCggg                                                   >  1:961611/1‑150 (MQ=255)
                                                                                              ccATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGTTGCAGGCCAGGACCAGGGCATTGAGGGTCAGCGGGTAGCTTTCCGGGc                                                  >  2:1279447/1‑150 (MQ=255)
                                                                                              ccATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGTTGCAGGCCAGGACCAGGGCATTGAGGGTCAGCGGGTAGCTTTCCGGGc                                                  >  2:304092/1‑150 (MQ=255)
                                                                                                   tGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGTTGCAGGCCAGGACCAGGGCATTGAGGGTCAGCGGGTAGCTTTCCGGGCtggtg                                             <  2:421943/150‑1 (MQ=255)
                                                                                                     cAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGTTGCAGGCCAGGACCAGGGCATTGAGGGTCAGCGGGTAGCTTTCCGGGCTGGTGGc                                           >  1:816161/1‑150 (MQ=255)
                                                                                                      aGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCAc                                                                                                                         >  2:975820/1‑71 (MQ=255)
                                                                                                      aGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCAc                                                                                                                         <  1:975820/71‑1 (MQ=255)
                                                                                                       ggCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGTTGCAGGCCAGGACCAGGGCATTGAGGGTCAGCGGGTAGCTTTCCGGGCTGGTGGCCt                                         <  1:1279447/150‑1 (MQ=255)
                                                                                                       ggCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGTTGCAGGCCAGGACCAGGGCATTGAGGGTCAGCGGGTAGCTTTCCGGGCTGGTGGCCt                                         >  2:69781/1‑150 (MQ=255)
                                                                                                           ggTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGTTGCAGGCCAGGACCAGGGCATTGAGGGTCAGCGGGTAGCTTTc                                                       <  2:4711/132‑1 (MQ=255)
                                                                                                           ggTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGTTGCAGGCCAGGACCAGGGCATTGAGGGTCAGCGGGTAGCTTTc                                                       >  1:4711/1‑132 (MQ=255)
                                                                                                             tCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGTTGCAGGCCAGGACCAGGGCATTGAGGGTCAGCGGGTAGCTTTCCGGGCTGGTGGCCTGCtttt                                   <  2:1427454/150‑1 (MQ=255)
                                                                                                             tCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGTTGCAGGCCAGGACCAGGGCATTGAGGGTCAGCGGGTAGCTTTCCGGGCTGGTGGCCTGCtttt                                   <  2:816161/150‑1 (MQ=255)
                                                                                                                  cccTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGTTGCAGGCCAGGACCAGGGCATTGAGGGTCAGc                                                                  <  2:1283590/114‑1 (MQ=255)
                                                                                                                  cccTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGTTGCAGGCCAGGACCAGGGCATTGAGGGTCAGc                                                                  >  1:1283590/1‑114 (MQ=255)
                                                                                                                   ccTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGTTGCAGGCCAGGACCAGGGCATTGAGGGTCAGCGGGTAGCTTTCCGGGCTGGTGGCCTGCTTTTCGATCa                             <  2:1104713/150‑1 (MQ=255)
                                                                                                                     tGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGTTGCAGGCCAGGACCAGGGCATTGAGGGTCAGCGGGTAGCTTTCCGGGCTGGTGGCCTGCTTTTCGATCAGc                           <  1:1394392/150‑1 (MQ=255)
                                                                                                                         cTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTc                                                                                                                   >  2:930291/1‑58 (MQ=255)
                                                                                                                         cTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTc                                                                                                                   <  1:930291/58‑1 (MQ=255)
                                                                                                                         cTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGTTGCAGGCCAGGACCAGGGCATTGAGGGTCAGCGGGTAGCTTTCCGGGCTGGTGGCCTGCTTTTCGATCAGCGAGc                       >  2:1256638/1‑150 (MQ=255)
                                                                                                                         cTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGTTGCAGGCCAGGACCAGGGCATTGAGGGGCAGCGGGTAGCTTTCCGGGCTGGTGGCCTGCTTTTCGATCAGCGAGc                       >  1:1327195/1‑150 (MQ=255)
                                                                                                                          ttCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGTTGCa                                                                                              <  2:906722/78‑1 (MQ=255)
                                                                                                                          ttCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGTTGCa                                                                                              >  1:906722/1‑78 (MQ=255)
                                                                                                                          ttCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGTTGCAGGTCAGGACCAGGGCATTGAGGGTCAGCGGGTAGCTTTCCGGGCTGGTGGCCTGCTTTTCGATCAGCGAGcc                      >  1:567210/1‑150 (MQ=255)
                                                                                                                              aGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGTTGCAGGCCAGGACCAGGGCATTGAgg                                                                        <  2:567555/96‑1 (MQ=255)
                                                                                                                              aGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGTTGCAGGCCAGGACCAGGGCATTGAgg                                                                        >  1:567555/1‑96 (MQ=255)
                                                                                                                                 tcacccAAGGCCAGGTCGGCCAGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGTTGCAGGCCAGGACCAGGGCATTGAGGGTCAGCGGGTAGCCTCCCGTGCTGGTGGCCTGCTTTTCGATCAGCGAGCCGAGAAcg               >  2:187223/6‑150 (MQ=255)
                                                                                                                                           ccTGGCCGAACTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGTTGCAGGCCAGGACCAGGGCATTGAGGGTCAGCGGGTAGCTTTCCGGGCTGGTGGCCTGCTTTTCGATCAGCGAGCCGAGAACGCGGATCTCGa     <  2:1327195/150‑1 (MQ=255)
                                                                                                                                              ggCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGTTGCAGGCCAGGACCAGGGCATTGAGGGTCAGCGGGTAGCTTTCCGGGCTGGTGGCCTGCTTTTCGATCAGCGAGCCGAGAACGCGGATCTCGATGc  >  1:675141/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CCCCGCAGGAACATCAGGCCCATCAGCACCAACTGCGCTGGCACCAGCTCCAGTGCCTTGTCTACCCGATGTTCCCAGCGGTCGGCACGGCTGCCCATCTGCAGGCGGGTCATGCCCTGACCTTCGAGGGCGCGCAAGGCCTGGCCGACCTGGCCTTGGGTGAGGTTCATCACCGGCTCGCGGCTGGTCTTCTGGTTGCAGGCCAGGACCAGGGCATTGAGGGTCAGCGGGTAGCTTTCCGGGCTGGTGGCCTGCTTTTCGATCAGCGAGCCGAGAACGCGGATCTCGATGC  >  NC_002947/2789959‑2790250

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