Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 516,285 A→T 5.7% L237Q (CTG→CAG)  PP_0429 ← short‑chain dehydrogenase/reductase family oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947516,2850AT 5.7% 217.1 / 4.2 71L237Q (CTG→CAG) PP_0429short‑chain dehydrogenase/reductase family oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (33/32);  new base T (2/2);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.69e-01

ATAGATCGTCCACAGCACTTGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCCGCGCAATGCGCCGGCGTGCG  >  NC_002947/516136‑516434
                                                                                                                                                     |                                                                                                                                                     
aTAGATCGTCCACAGCACTTGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGgcagca                                                                                                                                                       <  2:387056/150‑1 (MQ=255)
aTAGATCGTCCACAGCACTTGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGgcagca                                                                                                                                                       <  2:348308/150‑1 (MQ=255)
     tCGTCCACAGCACTTGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCACCGCGGCTGCTGCcgc                                                                                                                                                  >  1:626038/1‑150 (MQ=255)
        tCCACAGCACTTGCTCCGTGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGt                                                                                                                                               >  2:285603/1‑150 (MQ=255)
          cacaGCACTTGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCACCGCGGCTGCTGCCGCGCTgg                                                                                                                                             <  2:626038/150‑2 (MQ=255)
                      tCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACgccagcca                                                                                                                                 <  2:1208343/150‑1 (MQ=255)
                        cATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCcagccagc                                                                                                                               <  1:946903/150‑1 (MQ=255)
                        cATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCcagccagc                                                                                                                               >  1:40333/1‑150 (MQ=255)
                        cATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCcagccagc                                                                                                                               <  1:224306/150‑1 (MQ=255)
                        cATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGAGGTGAGAACGCcagccagc                                                                                                                               >  1:1064090/1‑150 (MQ=255)
                         aTGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGaa                                                                                                                                          >  2:1272580/1‑138 (MQ=255)
                         aTGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGaa                                                                                                                                          <  1:1272580/138‑1 (MQ=255)
                         aTGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCg                                                                                                                              <  2:40333/150‑1 (MQ=255)
                              gAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTgc                                                                                                                         <  1:503590/150‑1 (MQ=255)
                                 ggTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCgcc                                                                                                                      >  1:429963/1‑150 (MQ=255)
                                    cGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGAGGTGAGAACGCCAGCCAGCGGTTGCgccggc                                                                                                                   <  2:1064090/150‑1 (MQ=255)
                                       gCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGtt                                                                                                                           <  1:1062619/139‑1 (MQ=255)
                                       gCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGtt                                                                                                                           >  2:1062619/1‑139 (MQ=255)
                                       gCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGccggccgg                                                                                                                <  2:579658/150‑1 (MQ=255)
                                                  tCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGaa                                                                                                                                          >  1:1370362/1‑113 (MQ=255)
                                                  tCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGaa                                                                                                                                          <  2:1370362/113‑1 (MQ=255)
                                                  tCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGt                                                                                                     >  2:1228826/1‑150 (MQ=255)
                                                   cAGCCTTGAGCCGCTCTCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGGCGGGCATGATCAGTg                                                                                                    <  1:49227/150‑1 (MQ=255)
                                                           aGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGt                                                                                                                                               <  2:315095/99‑1 (MQ=255)
                                                           aGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGt                                                                                                                                               >  1:315095/1‑99 (MQ=255)
                                                             ccGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCtt                                                                                          <  1:123879/150‑1 (MQ=255)
                                                              cGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTg                                                                                         >  2:61306/1‑150 (MQ=255)
                                                                cTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTc                                                                                       >  1:1109669/1‑150 (MQ=255)
                                                                 tCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTcc                                                                                      <  1:1159740/150‑1 (MQ=255)
                                                                    cccAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGt                                                                                   <  1:290711/150‑1 (MQ=255)
                                                                     ccAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTg                                                                                  >  1:1397112/1‑150 (MQ=255)
                                                                        gTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCAt                                                                                                            <  2:1181941/121‑1 (MQ=255)
                                                                        gTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCAt                                                                                                            >  1:1181941/1‑121 (MQ=255)
                                                                           gCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCtt                                                                                          >  1:604575/1‑136 (MQ=255)
                                                                           gCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCtt                                                                                          <  2:604575/136‑1 (MQ=255)
                                                                           gCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCa                                                                            <  1:61306/150‑1 (MQ=255)
                                                                            cTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCAc                                                                           <  1:605272/150‑1 (MQ=255)
                                                                              aTTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACgg                                                                         >  1:637334/1‑150 (MQ=255)
                                                                               ttCCATGCCTGAGGGGCAATCGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCGGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGGGCGGTTCTTGTCCAGTGCCCGCACggc                                                                        >  2:158719/1‑149 (MQ=255)
                                                                               ttCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGgt                                                                        >  2:951100/1‑150 (MQ=255)
                                                                                  cATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAAc                                                                                                                                         <  1:1169744/82‑1 (MQ=255)
                                                                                  cATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAAc                                                                                                                                         >  2:1169744/1‑82 (MQ=255)
                                                                                      ccTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAgc                                                                 <  2:637334/150‑1 (MQ=255)
                                                                                      ccTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAgc                                                                 <  2:1342969/150‑1 (MQ=255)
                                                                                       cTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCGAGGGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGAATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTCCAggg                                                                >  2:145367/1‑148 (MQ=255)
                                                                                        tGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAgcgc                                                               >  2:787592/1‑150 (MQ=255)
                                                                                          aGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCa                                                             <  2:429963/150‑1 (MQ=255)
                                                                                                      ggcggTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCGGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAgcgc                                                               <  2:529695/136‑1 (MQ=255)
                                                                                                      ggcggTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCGGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAgcgc                                                               >  1:529695/1‑136 (MQ=255)
                                                                                                            ggATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGc                                           >  2:284994/1‑150 (MQ=255)
                                                                                                                  ccccGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCgg                                                                                               >  1:857316/1‑92 (MQ=255)
                                                                                                                  ccccGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCgg                                                                                               <  2:857316/92‑1 (MQ=255)
                                                                                                                  ccccGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGcgcg                                     >  2:1159099/1‑150 (MQ=255)
                                                                                                                  ccccGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGcgcg                                     >  1:314584/1‑150 (MQ=255)
                                                                                                                    ccGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGc                                   >  2:108156/1‑150 (MQ=255)
                                                                                                                    ccGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGc                                   <  2:1397112/150‑1 (MQ=255)
                                                                                                                      ggCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCtt                                 <  1:1228826/150‑1 (MQ=255)
                                                                                                                          aGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCa                             <  1:105031/150‑1 (MQ=255)
                                                                                                                           gCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCAt                            >  1:487295/1‑150 (MQ=255)
                                                                                                                             ttGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCg                          >  1:229641/1‑150 (MQ=255)
                                                                                                                                cgcACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGc                       >  1:1198703/1‑150 (MQ=255)
                                                                                                                                cgcACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGc                       <  1:108156/150‑1 (MQ=255)
                                                                                                                                 gcACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGcc                      >  1:1139855/1‑150 (MQ=255)
                                                                                                                                   aCCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCTGCGGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCcgc                    <  2:550315/150‑1 (MQ=255)
                                                                                                                                   aCCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCcgc                    >  1:38927/1‑150 (MQ=255)
                                                                                                                                       gccagcGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCCGCGCaa                <  1:145367/150‑1 (MQ=255)
                                                                                                                                       gccagcGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCCGCGCaa                >  2:1039302/1‑150 (MQ=255)
                                                                                                                                        ccagcGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAgcgc                                                               >  2:972142/1‑102 (MQ=255)
                                                                                                                                        ccagcGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAgcgc                                                               <  1:972142/102‑1 (MQ=255)
                                                                                                                                        ccagcGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCCGCGCAAt               <  1:951100/150‑1 (MQ=255)
                                                                                                                                          accgcgGCTGCTGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCtt                                 <  2:216515/128‑1 (MQ=255)
                                                                                                                                          accgcgGCTGCTGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCtt                                 >  1:216515/3‑130 (MQ=255)
                                                                                                                                           gcgcGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCCGCGCAATgcg            >  1:1216808/1‑150 (MQ=255)
                                                                                                                                                   gcagcCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCCGCGCAATGCGCCGGCGTg    <  1:1051504/150‑1 (MQ=255)
                                                                                                                                                    cagcCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCCGCGCAATGCGCCGGCGTGc   <  1:787592/150‑1 (MQ=255)
                                                                                                                                                     agcCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCCGCGCAATGCGCCGGCGTGCg  >  1:396396/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATAGATCGTCCACAGCACTTGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCCGCGCAATGCGCCGGCGTGCG  >  NC_002947/516136‑516434

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: