Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,522,285 C→G 5.3% F30L (TTC→TTG) ‡ idh → isocitrate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,522,2850CG 5.3% 239.8 / 3.0 75F30L (TTC→TTG) ‡idhisocitrate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (32/39);  new base G (2/2);  total (34/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

TGCGCCTACCCTGCGCACCACGAGACTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTGGCTACCACCCCTGAA  >  NC_002947/4522136‑4522429
                                                                                                                                                     |                                                                                                                                                
tGCGCCTACCCTGCGCACCACGAGACTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGACGCGTTc                                                                                                                                                  >  2:624321/1‑150 (MQ=255)
  cgcCTACCCTGCGCACCACGAGACTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCAc                                                                                                                                                >  1:1382834/1‑150 (MQ=255)
  cgcCTACCCTGCGCACCACGAGACTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCAc                                                                                                                                                <  2:1803825/150‑1 (MQ=255)
   gcCTACCCTGCGCACCACGAGACTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCAcc                                                                                                                                               >  2:620788/1‑150 (MQ=255)
   gcCTACCCTGCGCACCACGAGACTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTAGCTGCTGCCGATCATCGAAGCGTTCAcc                                                                                                                                               >  1:1268669/1‑150 (MQ=255)
          cTGCGCACCACGAGACTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGGTGAACGCTTCga                                                                                                                                        <  2:967559/150‑2 (MQ=255)
                accacGAGACTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACa                                                                                                                                  >  2:1551900/1‑150 (MQ=255)
                  cacGAGACTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCACCGATCATCTATACCTTCACCCGCTAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATc                                                                                                                                >  1:1205120/1‑150 (MQ=255)
                   acGAGACTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCg                                                                                                                               <  1:729779/150‑1 (MQ=255)
                    cGAGACTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGc                                                                                                                              >  1:1158844/1‑150 (MQ=255)
                     gagaCTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCt                                                                                                                             <  1:1598651/150‑1 (MQ=255)
                             aGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAAc                                                                                                                     >  1:1086339/1‑150 (MQ=255)
                                aCGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCGTCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGGTGAACGCTTCgat                                                                                                                                       <  1:749818/129‑3 (MQ=255)
                                aCGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACTAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTcg                                                                                                                  >  1:276275/1‑150 (MQ=255)
                                aCGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGGTGAACGCTTCgat                                                                                                                                       >  2:749818/1‑127 (MQ=255)
                                  gCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTcgcg                                                                                                                >  1:1518468/1‑150 (MQ=255)
                                       gCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACAtt                                                                                                           >  2:1824836/1‑150 (MQ=255)
                                         acacaGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTc                                                                                                         >  2:1258228/1‑150 (MQ=255)
                                           acaGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTccc                                                                                                       >  2:915512/1‑150 (MQ=255)
                                              gagagTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCctgg                                                                                                    <  2:1086339/150‑1 (MQ=255)
                                                gagTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGACACTCGCGACATTTCCctggct                                                                                                  <  1:537054/150‑1 (MQ=255)
                                                gagTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCctggct                                                                                                  <  2:1132329/150‑1 (MQ=255)
                                                  gTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCctggctgg                                                                                                <  1:1824836/150‑1 (MQ=255)
                                                   ttAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTccc                                                                                                       <  2:721144/142‑1 (MQ=255)
                                                   ttAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTccc                                                                                                       >  1:721144/1‑142 (MQ=255)
                                                   ttAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCtggctggc                                                                                               >  2:1702749/1‑150 (MQ=255)
                                                   ttAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCtggctggc                                                                                               <  1:1772849/150‑1 (MQ=255)
                                                    tAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGcc                                                                                              <  1:1702749/150‑1 (MQ=255)
                                                    tAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACCCGCTCCATATCCCTggccggcc                                                                                              <  2:567320/150‑1 (MQ=255)
                                                         cGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATc                                                                                         <  1:684794/150‑1 (MQ=255)
                                                            aTGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCTCTTCGGCTGACATCGCTGTCGAAACTAGCGACATTTCCCTGGCTGGCCGTATCCtt                                                                                      <  1:1551900/150‑1 (MQ=255)
                                                              gCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTgc                                                                                    <  2:1056681/150‑1 (MQ=255)
                                                                    ccGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCtt                                                                              <  1:1434043/150‑1 (MQ=255)
                                                                      gTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTcc                                                                            <  1:1258228/150‑1 (MQ=255)
                                                                       ttCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTccc                                                                           <  2:1814711/150‑1 (MQ=255)
                                                                        tCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCg                                                                          <  1:1722695/150‑1 (MQ=255)
                                                                          cAAGATCATCTATACCTTCACCGACTAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGa                                                                        <  2:276275/150‑1 (MQ=255)
                                                                          cAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGa                                                                        <  1:240970/150‑1 (MQ=255)
                                                                          cAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGa                                                                        <  1:517168/150‑1 (MQ=255)
                                                                            aGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGGCATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGc                                                                      <  2:1107174/150‑1 (MQ=255)
                                                                               tcatcTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAAc                                                                   >  2:1540102/1‑150 (MQ=255)
                                                                                   cTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTggg                                                               <  1:624321/150‑1 (MQ=255)
                                                                                       aCCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCa                                                           <  1:620788/150‑1 (MQ=255)
                                                                                        ccTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCag                                                          >  2:1831973/1‑150 (MQ=255)
                                                                                         cTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCct                                                                                                      <  2:372802/105‑1 (MQ=255)
                                                                                         cTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCct                                                                                                      >  1:372802/1‑105 (MQ=255)
                                                                                         cTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCaga                                                         >  1:1704127/1‑150 (MQ=255)
                                                                                          ttCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCagag                                                        >  2:1756395/1‑150 (MQ=255)
                                                                                            cACCGACGAAGCCCCCGCCCTCGCCACCTCCTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTNGCTGGCCGTNTCCTTGCCNCCTTCCCGGAGCAACTGGGCGCAgagan                                                      <  1:743874/150‑2 (MQ=255)
                                                                                               cgacgaAGCCCCCGCCCTCGCCACTTACTCGCTGCTTCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCa                                                   <  1:186023/150‑1 (MQ=255)
                                                                                                aacgcggCCCGCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTTACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCaa                                                  >  2:965147/7‑150 (MQ=255)
                                                                                                    aaGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAg                                              <  1:915512/150‑1 (MQ=255)
                                                                                                     aGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAgg                                             <  2:322031/150‑1 (MQ=255)
                                                                                                        ccccGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCaga                                                         >  2:8731/1‑135 (MQ=255)
                                                                                                        ccccGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCaga                                                         <  1:8731/135‑1 (MQ=255)
                                                                                                                tCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCAc                                      <  2:943129/146‑1 (MQ=255)
                                                                                                                tCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCAc                                      >  1:943129/1‑146 (MQ=255)
                                                                                                                tCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTgg                                  <  2:987992/150‑1 (MQ=255)
                                                                                                                  gCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCg                                >  2:90086/1‑150 (MQ=255)
                                                                                                                   ccACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCaa                                                                    <  1:971739/113‑1 (MQ=255)
                                                                                                                   ccACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCaa                                                                    >  2:971739/1‑113 (MQ=255)
                                                                                                                         aCTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTgg                         <  2:1518468/150‑1 (MQ=255)
                                                                                                                          cTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTggg                        >  1:490035/1‑150 (MQ=255)
                                                                                                                           tCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGc                       >  1:472065/1‑150 (MQ=255)
                                                                                                                             gctgctGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTgg                                                                <  2:1936956/107‑1 (MQ=255)
                                                                                                                             gcggctgcCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTgg                                                                >  1:1936956/4‑107 (MQ=255)
                                                                                                                                  tgcCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTgg                <  2:1917027/150‑1 (MQ=255)
                                                                                                                                      gATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTGGCTac            >  2:31468/1‑150 (MQ=255)
                                                                                                                                      gATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTGGCTac            >  2:1579300/1‑150 (MQ=255)
                                                                                                                                        cagccGAAGCGGTGAACGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTGGCTacca          >  1:1282368/5‑150 (MQ=255)
                                                                                                                                           tcGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTccc                                                                           >  2:744433/1‑82 (MQ=255)
                                                                                                                                           tcGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTccc                                                                           <  1:744433/82‑1 (MQ=255)
                                                                                                                                             gAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTGGCTACCACCCCt     <  1:31468/150‑1 (MQ=255)
                                                                                                                                             gAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTGGCTACCACCCCt     >  2:978782/1‑150 (MQ=255)
                                                                                                                                               aGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTGGCTACCACCCCTGa   >  2:1569814/1‑150 (MQ=255)
                                                                                                                                                gCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTGGCTACCACCCCTGaa  <  2:1704127/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
TGCGCCTACCCTGCGCACCACGAGACTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTGGCTACCACCCCTGAA  >  NC_002947/4522136‑4522429

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: ‑