Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,686,436 C→T 7.7% L282L (CTG→TTG)  PP_3251 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,686,4360CT 7.7% 227.6 / 6.4 78L282L (CTG→TTG) PP_3251hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (39/33);  new base T (3/3);  total (42/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 = 6.58e-01

CCACGCCCGCCTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCGCATTACAGT  >  NC_002947/3686297‑3686583
                                                                                                                                             |                                                                                                                                                   
ccACGCCCGCCTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGc                                                                                                                                             >  5:369913/1‑148 (MQ=255)
  aCGCCCGCCTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGaa                                                                                                                                          <  4:305121/149‑1 (MQ=255)
  aCGCCCGCCTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGaa                                                                                                                                          >  3:255009/1‑149 (MQ=255)
  aCGCCCGCCTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGaa                                                                                                                                          >  7:286755/1‑149 (MQ=255)
     cccgccTGGCAATTCAGGTGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGCCTTGCAGGTCGAACTg                                                                                                                                       <  8:286517/149‑1 (MQ=255)
      ccgccTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGt                                                                                                                                      >  1:216328/1‑149 (MQ=255)
       cgccTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGCCTTGCAGGTCGAACTgg                                                                                                                                      >  4:406548/1‑147 (MQ=255)
       cgccTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGCCTTGCAGGTCGAACTgg                                                                                                                                      <  3:406548/148‑2 (MQ=255)
          cTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCaa                                                                                                                                                >  8:284905/1‑135 (MQ=255)
          cTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTg                                                                                                                                   >  3:244780/1‑148 (MQ=255)
               aaTTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGCCTTGCAGGTCGAACTgg                                                                                                                                      <  8:333440/140‑2 (MQ=255)
               aaTTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGCCTTGCAGGTCGAACTgg                                                                                                                                      >  7:333440/1‑139 (MQ=255)
               aaTTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGccc                                                                                                                             >  1:100903/1‑149 (MQ=255)
               aaTTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGACAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGccc                                                                                                                             >  7:113179/1‑149 (MQ=255)
                   cAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTg                                                                                                                          >  4:276207/1‑148 (MQ=255)
                    aGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGaa                                                                                                                                          >  7:171932/1‑131 (MQ=255)
                    aGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGaa                                                                                                                                          <  8:171932/131‑1 (MQ=255)
                        gcTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGtt                                                                                                                                     >  6:81843/1‑132 (MQ=255)
                        gcTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGtt                                                                                                                                     <  5:81843/132‑1 (MQ=255)
                         cTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCaa                                                                                                                   <  8:286755/149‑1 (MQ=255)
                                 gCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTc                                                                                                            >  1:390469/1‑148 (MQ=255)
                                       ccccGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGcc                                                                                                     <  1:271045/149‑1 (MQ=255)
                                       ccccGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGc                                                                                                      >  7:52120/1‑148 (MQ=255)
                                       ccccGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTCTTTGCAGCCCTTGAACGCAAGCGCGTCTATGGc                                                                                                      >  5:389699/1‑148 (MQ=255)
                                          cgcgCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAAt                                                                                                  <  3:276207/149‑1 (MQ=255)
                                              cATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCTTCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTg                                                                                                <  1:328592/149‑1 (MQ=255)
                                              cATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGat                                                                                              <  6:369913/149‑1 (MQ=255)
                                               aTGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGatg                                                                                             <  8:52120/149‑1 (MQ=255)
                                                        aGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGAc                                                                                    >  5:228470/1‑149 (MQ=255)
                                                            gTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCAtt                                                                                <  6:278556/149‑1 (MQ=255)
                                                             ttCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTg                                                                               <  5:365432/149‑1 (MQ=255)
                                                             ttCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑AGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTg                                                                               >  2:186178/1‑149 (MQ=255)
                                                                tACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCAtt                                                                                <  1:389072/145‑1 (MQ=255)
                                                                tACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCAtt                                                                                >  2:389072/1‑145 (MQ=255)
                                                                 aCGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTgccgc                                                                           >  4:277732/1‑149 (MQ=255)
                                                                   gCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCt                                                                         >  4:391765/1‑149 (MQ=255)
                                                                      aGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTgca                                                                      <  4:255009/149‑1 (MQ=255)
                                                                         gcgcgcTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTgcagca                                                                   <  6:228470/149‑1 (MQ=255)
                                                                           gcgcTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCagcagc                                                                 <  4:37475/149‑1 (MQ=255)
                                                                              cTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCa                                                              <  4:341517/149‑1 (MQ=255)
                                                                                  ggCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGa                                                           <  4:244780/148‑1 (MQ=255)
                                                                                  ggCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGa                                                           <  3:391765/148‑1 (MQ=255)
                                                                                     tACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGatcaa                                                       >  5:191859/1‑149 (MQ=255)
                                                                                       ccGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGatcaatc                                                     <  8:113179/149‑1 (MQ=255)
                                                                                         gCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGAtcaatcaa                                                   <  2:390469/149‑1 (MQ=255)
                                                                                         gCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTCTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGAtcaatcaa                                                   >  8:140338/1‑149 (MQ=255)
                                                                                          ccTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATcaatcaat                                                  <  6:191859/149‑1 (MQ=255)
                                                                                              aaTGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAAtgg                                                <  7:232436/147‑1 (MQ=255)
                                                                                              aaTGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAAtgg                                                <  6:6788/147‑1 (MQ=255)
                                                                                                tGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAAtggctgg                                            <  2:100903/149‑1 (MQ=255)
                                                                                                  aCCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATcaatcaat                                                  <  5:106828/141‑1 (MQ=255)
                                                                                                  aCCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATcaatcaat                                                  >  6:106828/1‑141 (MQ=255)
                                                                                                     ggggTATCGCGACGCCAGAGCTCATC‑‑CGGCAGTTCGACCTGCTATGCGAACTGTTTGCAGCCCTTGAACGCATGCGGGTCTATGTCCAATTGATGTATGCCGACCATTGCCGCCTgctt                                                                     <  6:386589/119‑3 (MQ=255)
                                                                                                     ggggTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGTATGCCGACCATTGCCGCCTgcag                                                                     >  5:386589/1‑119 (MQ=255)
                                                                                                     ggggTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGaa                                       <  2:216328/149‑1 (MQ=255)
                                                                                                             gccacgccaGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGgcgc                               >  2:379397/1‑149 (MQ=255)
                                                                                                                acgccaGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCAt                            >  2:131301/1‑149 (MQ=255)
                                                                                                                acgccaGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCAt                            >  7:257789/1‑149 (MQ=255)
                                                                                                                  gccaGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAAtggc                                               <  3:172707/128‑1 (MQ=255)
                                                                                                                  gccaGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGc                          <  2:377880/149‑1 (MQ=255)
                                                                                                                   ccaGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGcc                         >  3:77419/1‑149 (MQ=255)
                                                                                                                    caGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGccc                        <  6:389699/149‑1 (MQ=255)
                                                                                                                       aCCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCa                     >  4:390379/1‑149 (MQ=255)
                                                                                                                           cATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCg                 >  2:345176/1‑149 (MQ=255)
                                                                                                                             tC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCACCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGtt               >  5:368739/1‑149 (MQ=255)
                                                                                                                                caaacAGTTCGCCTCGCAGGTCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTg              >  4:140568/5‑149 (MQ=255)
                                                                                                                                  gCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCATTTGATGCATGCCGCCCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCTTTCGTTGcc            >  6:115247/1‑149 (MQ=255)
                                                                                                                                  gCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCagcagc                                                                 >  6:374059/1‑96 (MQ=255)
                                                                                                                                  gCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCagcagc                                                                 <  5:374059/96‑1 (MQ=255)
                                                                                                                                  gCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGcc            >  3:401756/1‑149 (MQ=255)
                                                                                                                                   ccAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCg           >  8:348239/1‑149 (MQ=255)
                                                                                                                                      gTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCGCAt        >  1:30771/1‑149 (MQ=255)
                                                                                                                                       ttCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCGCAtt       <  3:277732/149‑1 (MQ=255)
                                                                                                                                       ttCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCGCAt        >  3:80643/1‑148 (MQ=255)
                                                                                                                                       ttCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCGCAt        >  3:86203/1‑148 (MQ=255)
                                                                                                                                          gACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGc             <  7:45356/140‑1 (MQ=255)
                                                                                                                                          gACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGc             >  8:45356/1‑140 (MQ=255)
                                                                                                                                           aCCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCGCATTACAg   >  2:187764/1‑149 (MQ=255)
                                                                                                                                            ccTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCACCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCGCATTACAGt  <  6:368739/149‑1 (MQ=255)
                                                                                                                                             cTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCGCATTACAGt  >  5:408977/1‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
CCACGCCCGCCTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCGCATTACAGT  >  NC_002947/3686297‑3686583

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