Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,128,148 A→T 6.1% Q158L (CAG→CTG)  htrG → signal transduction protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,128,1480AT 6.1% 241.1 / 4.3 78Q158L (CAG→CTG) htrGsignal transduction protein
Reads supporting (aligned to +/- strand):  ref base A (40/33);  new base T (3/2);  total (43/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 = 9.98e-01

GCTCGATGCCCAGGTCACCGAACTGACCGGGCAGTTGAAGACCATCGACGACAGCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCG  >  NC_002947/4127999‑4128295
                                                                                                                                                     |                                                                                                                                                   
gCTCGATGCCCAGGTCACCGAACTGACCGGGCAGTTGAAGACCATCGACGACAGCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGagca                                                                                                                                                     <  2:518569/150‑1 (MQ=255)
         ccAGGTCACCGAACTGACCGGGCAGTTGAAGACCATGGACGACAGCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCa                                                                                                                                            >  1:1935014/1‑150 (MQ=255)
         ccAGGTCACCGAACTGACCGGGCAGTTGAAGACCATCGACGACAGCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCa                                                                                                                                            >  2:1632132/1‑150 (MQ=255)
         ccAGGTCACCGAACTGACCGGGCAGTTGAAGACCATCGACGACAGCTGGAAGAACCGCGTACAGGGCAAGCAGGAAACGCTGGATTCGCGGAAGCAGCTGATTGATGAGTTGGAGGCACGCGACAAGGCTCTGATGGAGCTGCTTGGTCa                                                                                                                                            >  2:963180/1‑150 (MQ=255)
         ccAGGTCACAGAACTGACCGGGCAGTTGAAGACCATCGACGACAGCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCa                                                                                                                                            >  2:31881/1‑150 (MQ=255)
           aGGTCACCGAACTGACCGGGCAGTTGAAGACCATCGACGACAGCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCaga                                                                                                                                          <  1:1603277/150‑1 (MQ=255)
             gTCACCGAACTGACCGGGCAGTTGAAGACCATCGACGACAGCGGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGAAGCAAAGAGAGTATTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGc                                                                                                                                        >  2:1817109/1‑150 (MQ=255)
                aCCGAACTGACCGGGCAGTTGAAGACCATCGACGACAGCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCa                                                                                                                                      >  2:1140820/1‑149 (MQ=255)
                aCCGAACTGACCGGGCAGTTGAAGACCATCGACGACAGCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCa                                                                                                                                      <  1:1140820/149‑1 (MQ=255)
                 ccGAACTGACCGGGCAGTTGAAGACCATCGACGACAGCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCTCTGATCGAGCTGCTCGTTCAGTGccttg                                                                                                                                    <  1:963180/150‑4 (MQ=255)
                    gaCTGACCGGGCAGTTGAAGACCATCGACGACAGCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCa                                                                                                                                            >  2:2181958/2‑139 (MQ=255)
                    gaCTGACCGGGCAGTTGAAGACCATCGACGACAGCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCa                                                                                                                                            <  1:2181958/138‑1 (MQ=255)
                        gACCGGGCAGTTGAAGACCATCGACGACAGCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCt                                                                                                                             >  2:1102528/1‑150 (MQ=255)
                           cGGGCAGTTGAAGACCATCGACGACAGCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTgcg                                                                                                                          >  1:2059372/1‑150 (MQ=255)
                            gggCAGTTGAAGACCATCGACGACAGCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACACAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTgcgc                                                                                                                         >  1:739847/1‑150 (MQ=255)
                             ggCAGTTGAAGACCATCGACGACAGCTGGAAGAAACTCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGcgcg                                                                                                                        >  1:2300357/1‑150 (MQ=255)
                                aGTTGAAGACCATCGACGACAGCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCaga                                                                                                                                          <  2:1321456/129‑1 (MQ=255)
                                aGTTGAAGACCATCGACGACAGCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCAATCAGAGCCAGTCGGACCTGCGCGATa                                                                                                                     <  1:1817109/150‑1 (MQ=255)
                                aGTTGAAGACCATCGACGACAGCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGATCATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCaga                                                                                                                                          >  1:1321456/1‑129 (MQ=255)
                                    gAAGACCATCGACGACAGCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCg                                                                                                                                  <  2:827795/133‑1 (MQ=255)
                                    gAAGACCATCGACGACAGCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCg                                                                                                                                  >  1:827795/1‑133 (MQ=255)
                                        aCCATCGACGACAGCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAg                                                                                                                                     <  1:2002302/126‑1 (MQ=255)
                                        aCCATCGACGACAGCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAg                                                                                                                                     >  2:2002302/1‑126 (MQ=255)
                                         ccATCGACGACAGCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGccc                                                                                                            <  2:1519168/150‑1 (MQ=255)
                                            tCGACGACAGCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCTGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCGCGcc                                                                                                         >  1:2004724/1‑150 (MQ=255)
                                                cgacAGCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTggg                                                                                                     <  1:31881/150‑1 (MQ=255)
                                                 gacAGCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACACAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTgggg                                                                                                    <  2:739847/150‑1 (MQ=255)
                                                     gCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGAGGATg                                                                                                <  2:276512/150‑1 (MQ=255)
                                                      cTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATga                                                                                               <  1:469779/150‑1 (MQ=255)
                                                        ggAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATgaga                                                                                             >  1:261696/1‑150 (MQ=255)
                                                         gaagaaCCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCTCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGaa                                                                                            <  2:2300357/150‑1 (MQ=255)
                                                              aCCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGc                                                                                       <  2:261696/150‑1 (MQ=255)
                                                               ccGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAg                                                                                                                >  2:1938493/1‑124 (MQ=255)
                                                               ccGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAg                                                                                                                <  1:1938493/124‑1 (MQ=255)
                                                               ccGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGAt                                                                                                 <  2:901766/139‑1 (MQ=255)
                                                               ccGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGAt                                                                                                 >  1:901766/1‑139 (MQ=255)
                                                                        gggCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATga                                                                                               <  2:1531560/132‑1 (MQ=255)
                                                                        gggCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATga                                                                                               >  1:1531560/1‑132 (MQ=255)
                                                                            aTGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATgcgc                                                                         >  2:1000984/1‑150 (MQ=255)
                                                                            aTGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATgcgc                                                                         >  1:298550/1‑150 (MQ=255)
                                                                            aTGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATgcgc                                                                         >  1:365507/1‑150 (MQ=255)
                                                                             tGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGtga                                                                                 <  1:360326/141‑1 (MQ=255)
                                                                             tGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGtga                                                                                 >  2:360326/1‑141 (MQ=255)
                                                                                aGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGt                                                                                                                                    >  2:2218072/1‑87 (MQ=255)
                                                                                aGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGt                                                                                                                                    <  1:2218072/87‑1 (MQ=255)
                                                                                        cTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTAt                                                             >  1:1824863/1‑150 (MQ=255)
                                                                                            aTTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATggcg                                                         >  1:475923/1‑150 (MQ=255)
                                                                                               cgcgCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGgca                                                      <  1:1632132/150‑1 (MQ=255)
                                                                                                gcgcAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCAATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGgcag                                                     >  2:664050/1‑150 (MQ=255)
                                                                                                 cgcAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGgcagc                                                    >  2:1995370/1‑150 (MQ=255)
                                                                                                      gcagcTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGc                                               >  2:1946491/1‑150 (MQ=255)
                                                                                                          cTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGt                                           <  1:1816265/150‑1 (MQ=255)
                                                                                                          cTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGt                                           >  1:1384017/1‑150 (MQ=255)
                                                                                                            gattgatGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGc                                         >  2:591753/1‑150 (MQ=255)
                                                                                                             attgatGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCa                                        >  2:1582885/1‑150 (MQ=255)
                                                                                                              ttgatGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAg                                       <  2:1824863/150‑1 (MQ=255)
                                                                                                                       tGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTgg                              <  2:365507/150‑1 (MQ=255)
                                                                                                                              aCGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGc                                     >  2:326447/1‑136 (MQ=255)
                                                                                                                              aCGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGc                                     <  1:326447/136‑1 (MQ=255)
                                                                                                                                gCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGa                     >  2:820297/1‑150 (MQ=255)
                                                                                                                                gCAACAAGGCACTGAACGAGCAGCTCGATCAAAGCCAGTCTGATCTGCGCGATACCCAGGCCCGCCTGGGGGAGGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGTCGGTGCAGGCTTGCTGGCGGGGCTGa                     <  1:1102528/150‑1 (MQ=255)
                                                                                                                                 caacaaGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGtga                                                                                 <  2:176288/89‑1 (MQ=255)
                                                                                                                                 caacaaGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGtga                                                                                 >  1:176288/1‑89 (MQ=255)
                                                                                                                                  aacaaGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGc                                                                                       <  2:1372668/82‑1 (MQ=255)
                                                                                                                                  aacaaGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGc                                                                                       >  1:1372668/1‑82 (MQ=255)
                                                                                                                                    caaGGCACTGAACGAGCAGCTCGATCCGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCt                 <  1:591753/150‑1 (MQ=255)
                                                                                                                                      aGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGGATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGc               >  1:1077685/1‑150 (MQ=255)
                                                                                                                                          aCTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGgtg                                                                >  1:1671598/1‑97 (MQ=255)
                                                                                                                                          aCTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGgtg                                                                <  2:1671598/97‑1 (MQ=255)
                                                                                                                                            tGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATggcggc                                                       >  1:1752913/1‑104 (MQ=255)
                                                                                                                                            tGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATggcggc                                                       <  2:1752913/104‑1 (MQ=255)
                                                                                                                                            tGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCgg            >  2:308686/1‑147 (MQ=255)
                                                                                                                                            tGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCgg            <  1:308686/147‑1 (MQ=255)
                                                                                                                                             gatcGAGCTGCTCGTTCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCgg                           <  2:1289521/128‑1 (MQ=255)
                                                                                                                                             gatcGAGCTGCTCGTTCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCgg                           >  1:1289521/4‑131 (MQ=255)
                                                                                                                                             gAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACt        <  1:1582885/150‑1 (MQ=255)
                                                                                                                                              aaCGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCAGTTGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTg       >  1:1919978/1‑150 (MQ=255)
                                                                                                                                                 gAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGc                                               >  2:1343083/1‑107 (MQ=255)
                                                                                                                                                 gAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGc                                               <  1:1343083/107‑1 (MQ=255)
                                                                                                                                                   gcagcTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCg  >  1:2311501/1‑150 (MQ=255)
                                                                                                                                                   gcagcTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCg  <  1:1946491/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCTCGATGCCCAGGTCACCGAACTGACCGGGCAGTTGAAGACCATCGACGACAGCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCG  >  NC_002947/4127999‑4128295

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