Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,089,301 A→T 8.5% K53* (AAG→TAG) ‡ fumC‑I → class 2 fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,089,3010AT 8.5% 212.0 / 8.2 72K53* (AAG→TAG) ‡fumC‑Iclass 2 fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (34/31);  new base T (2/4);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

ACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGGCCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTATCGTCAAGGCTGTCGAGCAATTGCTGGCCGAGGACTTCATCCAGCACTTCCCGGTGGATGTGTTCCAGACAGGCTCTGGCACCAGT  >  NC_002947/1089157‑1089450
                                                                                                                                                |                                                                                                                                                     
aCCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGGCCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTCCACTGTTGTTGGCCAAGGcc                                                                                                                                                  <  2:1128858/150‑1 (MQ=255)
 cccGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGGCCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCTTGGCCg                                                                                                                                                 >  2:211808/1‑150 (MQ=255)
  ccGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGGCCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGc                                                                                                                                                >  2:1387551/1‑150 (MQ=255)
            cGTGACAGCATGGGTGAGCTGCAGGTGCCGGCCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCggccaag                                                                                                                                      <  1:934905/150‑1 (MQ=255)
            cGTGACAGCATGGGTGAGCTGCAGGTGCCGGCCCAGGCCCTCTACGGCGCCCACACCCAGCGCGCGGTCGACATCTTCTCTTTCCGTCGCCTACTTCGAGTTGCGGACGTTGGTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCggccaag                                                                                                                                      <  1:26918/150‑1 (MQ=255)
             gTGACAGCATGGGTGAGCTGCAGGTGCCGGCCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCTTGGCCGCAGCggccttgg                                                                                                                                     <  1:211808/150‑5 (MQ=255)
             gTGAAAGCATGGGTGAGCTGCAGGTCCCGGCCCAGGCCCTCTACGGCGCCCACACCCCGCGCGCGGTCTACAACTTACCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCACTCGTGCTTTGTTGCTGGGCTTGGCCTCAGCggccctgg                                                                                                                                     <  2:520569/150‑5 (MQ=255)
              tGACAGCATGGGTGAGCTGCAGGTGCCGGCCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCTTGGCCGCAGCggccttggc                                                                                                                                    <  1:557069/150‑6 (MQ=255)
                aCAGCATGGGTGAGCTGCAGGTGCCGGCCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCTTGGCCGCAGCggccttggcca                                                                                                                                  <  2:588066/150‑8 (MQ=255)
                     aTGGGTGAGCTGCAGGTGCCGGCCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTc                                                                                                                             <  1:289075/150‑1 (MQ=255)
                         gTGAGCTGCAGGTGCCGGCCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAAc                                                                                                                         >  1:165956/1‑150 (MQ=255)
                          tGAGCTGCAGGTGCCGGCCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACt                                                                                                                        >  1:146630/1‑150 (MQ=255)
                           gAGCTGCAGGTGCCGGCCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTg                                                                                                                       <  1:1008498/150‑1 (MQ=255)
                              cTGCAGGTGCCGGCCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAg                                                                                                                    >  1:1068035/1‑150 (MQ=255)
                                 cAGGTGCCGGCCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCaa                                                                                                                 <  1:613722/150‑1 (MQ=255)
                                 cAGGTGCCGGCCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCaa                                                                                                                 >  1:25179/1‑150 (MQ=255)
                                       ccGGCCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCATCGGTCAGCGTATGCCGGCCCAGTTCATTCGTTCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTcg                                                                                                           >  2:125886/1‑149 (MQ=255)
                                         ggCCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCgg                                                                                                                                           <  2:84046/116‑1 (MQ=255)
                                         ggCCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCgg                                                                                                                                           >  1:84046/1‑116 (MQ=255)
                                           cccAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCg                                                                                                                                            <  2:507009/113‑1 (MQ=255)
                                           cccAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCg                                                                                                                                            >  1:507009/1‑113 (MQ=255)
                                            ccAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCTTGGCCGCAGCggccttggccag                                                                                                                                 >  1:371738/1‑122 (MQ=255)
                                               ggCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCa                                                                                                   <  2:146630/150‑1 (MQ=255)
                                                 cccTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAgg                                                                                                 >  1:1286745/1‑150 (MQ=255)
                                                  ccTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAggg                                                                                                >  2:1392416/1‑150 (MQ=255)
                                                   ctctACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGTCGAAGGCCTCTGCGGCCAAGCCCAACATCGAACTGGAGCAACCGTCTGCCTGGCAGGGc                                                                                               <  2:90567/150‑1 (MQ=255)
                                                      tACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAg                                                                                            <  2:96847/150‑1 (MQ=255)
                                                      tACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAg                                                                                            <  2:957709/150‑1 (MQ=255)
                                                          gcgcCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTa                                                                                        >  1:390161/1‑150 (MQ=255)
                                                            gcCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTATc                                                                                      <  2:1188232/150‑1 (MQ=255)
                                                             cccAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCggc                                                                                                                                          <  1:831717/97‑1 (MQ=255)
                                                             cccAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCggc                                                                                                                                          >  2:831717/1‑97 (MQ=255)
                                                              ccAGACCGAGCGCGCGGTCGACAACTTCCCGAGCAGCGGCCAGCGTATGCCGGCCCAGTTCATGCGTGCCGTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTATCGt                                                                                    >  2:447050/1‑150 (MQ=255)
                                                              ccAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTATCGt                                                                                    >  2:1101239/1‑150 (MQ=255)
                                                                  aCCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTATCGTCAAg                                                                                <  2:25179/150‑1 (MQ=255)
                                                                   cccAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTATCGTCAAgg                                                                               <  1:1392416/150‑1 (MQ=255)
                                                                             ggTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGt                                                                                                             <  1:1184136/110‑1 (MQ=255)
                                                                             ggTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGt                                                                                                             >  2:1184136/1‑110 (MQ=255)
                                                                             ggTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTATCGTCAAGGCTgagatcgg                                                                     >  1:228241/1‑143 (MQ=255)
                                                                             ggTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCCGCCGGGCAGGGCCAGGCTATCGTCAAGGCTGTCGAgga                                                                     >  1:435143/1‑148 (MQ=255)
                                                                                   cAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTATCGTCAAGGCTGTCGAGCAATTGCt                                                               >  1:1318543/1‑150 (MQ=255)
                                                                                    aaCTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTATCGTCAAGGCTGTCGAGCAATTGCTg                                                              >  2:909099/1‑150 (MQ=255)
                                                                                     aCTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTATCGTCAAGGCTGTCGAGCAATTGCTgg                                                             >  2:689576/1‑150 (MQ=255)
                                                                                       ttCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTATCGTCAAGGCTGTCGAGCAATTGCTGGcc                                                           >  2:478388/1‑150 (MQ=255)
                                                                                          ccGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGc                                                                                               <  1:99272/111‑1 (MQ=255)
                                                                                          ccGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGc                                                                                               >  2:99272/1‑111 (MQ=255)
                                                                                            gATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTATCGTCAAGGCTGTCGAGCAATTGCTGGCCGAGGa                                                      <  1:1387551/150‑1 (MQ=255)
                                                                                                aGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTATCGTCAAGGCTGTCGAGCAATTGCTGGCCGAGGACTTc                                                  <  2:1318543/150‑1 (MQ=255)
                                                                                                  cGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTATCGTCAAGGCTGTCGAGCAATTGCTGGCCGAGGACTTCAt                                                >  2:376135/1‑150 (MQ=255)
                                                                                                   ggCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTATCGTCAAGGCTGTCGAGCAATTGCTGGCCGAGGACTTCATc                                               <  2:1038755/150‑1 (MQ=255)
                                                                                                   ggCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTATCGTCAAGGCTGTCGAGCAATTGCTGGCCGAGGACTTCATc                                               >  1:1038755/1‑150 (MQ=255)
                                                                                                     ccAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTATCGTCAAGGCTGTCGAGCAATTGCTGGCCGAGGACTTCATCCa                                             >  1:1349317/1‑150 (MQ=255)
                                                                                                           tATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTATCGTCAAGGCTGTCGAGCAATTGCTGGCCGAGGACTTCATCCAGCACtt                                       <  1:447050/150‑1 (MQ=255)
                                                                                                                  tccccgctCATTCGTGGCTTCTCGCTGCCCGAGGCCGCTGCGGCCAAGGCGAACGTCGAACTGGAGCAACTGTCTGCCGAGCAGGGCCAGGCTATCGGCAAGGCTGTCGAGCAATTGCGGGCCGAGGACTTCATCCAGCACTGCCCGGTg                                <  2:435143/143‑1 (MQ=255)
                                                                                                                  gCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTATCGTCAAGGCTGTCGAGCAATTGCTGGCCGAGGACTTCATCCAGCACTTCCCGGTg                                <  2:1286745/150‑1 (MQ=255)
                                                                                                                       gTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTATCGTCAAGGCTGTCGAGCAATTGCTGGCCGAGGa                                                      <  2:21615/123‑1 (MQ=255)
                                                                                                                       gTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTATCGTCAAGGCTGTCGAGCAATTGCTGGCCGAGGa                                                      >  1:21615/1‑123 (MQ=255)
                                                                                                                        ttCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTATCGTCAAGGCTGTCGAGCAATTGc                                                                >  1:1052925/1‑112 (MQ=255)
                                                                                                                        ttCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTATCGTCAAGGCTGTCGAGCAATTGc                                                                <  2:1052925/112‑1 (MQ=255)
                                                                                                                          cATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTATCGTCAAGGCTGTCGAGCAATTGCTGGCCGAGGACTTCATCCAGCACTTCCCGGTGGATGTGtt                        <  1:867170/150‑1 (MQ=255)
                                                                                                                                    ttgttgCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTATCGTCAAGGCTGTCGAGCAATTGCTGGCCGAGGACTTCATCCAGCACTTCCCGGTGGATGTGTTCCAGCCAGGc              >  2:658942/1‑150 (MQ=255)
                                                                                                                                       ttgCTGTCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGCAACAACTGTCTGCCGTTCTGTGCCAGGCTATCGTCAAGGCTGTCGAGCAATTGCTGGCCGAGGACTTCATCCAGCACTTCCCGGTGGATGTGTTCCAGACAGGctct           >  2:409932/1‑150 (MQ=255)
                                                                                                                                         gCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTATCGTCAAGGCTGTCGAGCAATTGCTGGCCGAGGACTTCATCCAGCACTTCCCGGTGGATGTGTTCCAGACAGGCTCTgg         <  1:909099/150‑1 (MQ=255)
                                                                                                                                          cTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTATCGTCAAgg                                                                               <  2:230200/79‑1 (MQ=255)
                                                                                                                                          cTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTATCGTCAAgg                                                                               >  1:230200/1‑79 (MQ=255)
                                                                                                                                            ggCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTATCGTCAAGGCTGTCGAGCAATTGCTGGCCGAGGACTTCATCCAGCACTTCCCGGTGGATGTGTTCCAGACAGGCTCTGGCAc      <  1:478388/150‑1 (MQ=255)
                                                                                                                                            ggCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTATCGTCAAGGCTGTCGAGCAATTGCTGGCCGAGGACTTCATCCAGCACTTCCCGGTGGATGTGTTCCAGACAGGCTCTGGCAc      >  1:832316/1‑150 (MQ=255)
                                                                                                                                               cAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTATc                                                                                      <  2:689413/67‑1 (MQ=255)
                                                                                                                                               cAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTATc                                                                                      >  1:689413/1‑67 (MQ=255)
                                                                                                                                               cAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTATCGTCAAGGCTGTCGAGCAATTGCTGGCCGAGGACTTCATCCAGCACTTCCCGGTGGATGTGTTCCAGACAGGCTCTGGCAc      <  1:645574/147‑1 (MQ=255)
                                                                                                                                               cAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTATCGTCAAGGCTGTCGAGCAATTGCTGGCCGAGGACTTCATCCAGCACTTCCCGGTGGATGTGTTCCAGACAGGCTCTGGCAc      >  2:645574/1‑147 (MQ=255)
                                                                                                                                               cAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTATCGTCAAGGCTGTCGAGCAATTGCTGGCCGAGGACTTCATCCAGCACTTCCCGGTGGATGTGTTCCAGACAGGCTCTGGCACCAg   >  2:1114660/1‑150 (MQ=255)
                                                                                                                                                aaGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTATCGTCAAGGCTGTCGAGCAATTGCTGGCCGAGGACTTCATCCAGCACTTCCCGGGGGATGTGGTCGATACAGGCTCTGGcctccgt  >  2:299938/1‑144 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
ACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGGCCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTATCGTCAAGGCTGTCGAGCAATTGCTGGCCGAGGACTTCATCCAGCACTTCCCGGTGGATGTGTTCCAGACAGGCTCTGGCACCAGT  >  NC_002947/1089157‑1089450

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