Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,039,262 G→C 9.3% R45G (CGC→GGC)  PP_0897 ← class 1 fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,039,2620GC 9.3% 278.3 / 14.9 94R45G (CGC→GGC) PP_0897class 1 fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base G (41/44);  new base C (6/3);  total (47/47)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.81e-01

TGTCGTCGACGCTCATGGTGGCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCG  >  NC_002947/1039115‑1039402
                                                                                                                                                   |                                                                                                                                            
tGTCGTCGACGCTCATGGTGGCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGt                                                                                                                                             >  1:301141/1‑149 (MQ=255)
  tcgtcgACGCTCATGGTGGCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTg                                                                                                                                            >  2:157830/1‑148 (MQ=255)
    gtcgACGCTCATGGTGGCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTgcag                                                                                                                                         >  7:101324/1‑149 (MQ=255)
        aCGCTCATGGTGGCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCGGCAgg                                                                                                                                     >  4:97586/1‑149 (MQ=255)
         cGCTCATGGTGGCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGcagcag                                                                                                                                      <  1:157830/147‑1 (MQ=255)
           cTCATGGTGGCGCCATCCGAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCTGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGa                                                                                                                                  <  6:242586/149‑1 (MQ=255)
              aTGGTGGCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTc                                                                                                                               >  6:251095/1‑149 (MQ=255)
               tggtggCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTcc                                                                                                                              <  7:210399/149‑1 (MQ=255)
                 gtggCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCCTGCTGCACGCGATTccat                                                                                                                            <  7:124550/149‑3 (MQ=255)
                  tggCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCCTGCTGCACGCGATTccatc                                                                                                                           <  6:97508/149‑4 (MQ=255)
                    gcgcCATCCCAGCGCACGTCCACGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACg                                                                                                                         >  8:176582/1‑149 (MQ=255)
                     cgcCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGc                                                                                                                        <  2:339832/149‑1 (MQ=255)
                          tCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGa                                                                                                                     >  2:144404/1‑147 (MQ=255)
                          tCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATa                                                                                                                   <  8:148610/149‑1 (MQ=255)
                            ccAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAgg                                                                                                                 >  2:309861/1‑149 (MQ=255)
                            ccAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAg                                                                                                                  <  7:203827/148‑1 (MQ=255)
                            ccAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAg                                                                                                                  <  2:255989/148‑1 (MQ=255)
                            ccAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAg                                                                                                                  <  5:213787/148‑1 (MQ=255)
                                cgcACGTCCATGCCCCCGCGGATGAACACGGTGACGATACCGGTCTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTc                                                                                                             >  6:324190/1‑149 (MQ=255)
                                 gcACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCg                                                                                                            >  1:340479/1‑149 (MQ=255)
                                  cACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGt                                                                                                           <  5:150455/149‑1 (MQ=255)
                                  cACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGt                                                                                                           <  4:47723/149‑1 (MQ=255)
                                    cGTCCACGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGc                                                                                                         <  7:176582/149‑1 (MQ=255)
                                     gTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCa                                                                                                        <  1:221894/149‑1 (MQ=255)
                                      tCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGAAt                                                                                                       >  8:104497/1‑149 (MQ=255)
                                       ccATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATg                                                                                                      >  7:1608/1‑149 (MQ=255)
                                       ccATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATg                                                                                                      >  7:376547/1‑149 (MQ=255)
                                       ccATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCAGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATg                                                                                                      >  1:341596/1‑149 (MQ=255)
                                            cccACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGTCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCt                                                                                                  <  5:183530/148‑1 (MQ=255)
                                              cACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGa                                                                                               >  3:257924/1‑149 (MQ=255)
                                              cACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGa                                                                                               <  2:340479/149‑1 (MQ=255)
                                              cACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGa                                                                                               <  2:249465/149‑1 (MQ=255)
                                               aCGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGAt                                                                                              <  5:336003/149‑1 (MQ=255)
                                                  cgGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGaa                                                                                           <  1:318506/149‑1 (MQ=255)
                                                   ggATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGaaa                                                                                          <  7:15377/149‑1 (MQ=255)
                                                     aTGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATc                                                                                        <  1:76711/149‑1 (MQ=255)
                                                              gTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGggtg                                                                               >  1:344802/1‑149 (MQ=255)
                                                                 aCGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGgta                                                                            <  8:1608/149‑1 (MQ=255)
                                                                   gATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTg                                                                                                 <  2:266179/126‑1 (MQ=255)
                                                                   gATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTg                                                                                                 >  1:266179/1‑126 (MQ=255)
                                                                   gATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGgtag                                                                           >  1:299613/1‑148 (MQ=255)
                                                                   gATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGgtag                                                                           >  5:84891/1‑148 (MQ=255)
                                                                   gATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGgtag                                                                           >  1:301160/1‑148 (MQ=255)
                                                                   gATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGgtag                                                                           >  7:148507/1‑148 (MQ=255)
                                                                     tACAGGTGTCCTGGCAGATTGGACGAGGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTAc                                                                        <  8:256842/149‑1 (MQ=255)
                                                                       ccGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACga                                                                      <  8:101324/149‑1 (MQ=255)
                                                                       ccGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACga                                                                      <  8:160271/149‑1 (MQ=255)
                                                                       ccGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACga                                                                      <  8:376547/149‑1 (MQ=255)
                                                                        cGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTAc                                                                        >  1:386649/1‑146 (MQ=255)
                                                                        cGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTAc                                                                        >  3:63338/1‑146 (MQ=255)
                                                                         ggTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACgaga                                                                    <  2:301141/149‑1 (MQ=255)
                                                                          gtgtCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGAt                                                                   >  3:100748/1‑149 (MQ=255)
                                                                                tgtcAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCCTGCAGCAGGCGATTCCTCACGCAGATAGGCCTAGTGCATGGCCTGCATGAAATCGAAGGGGTGGTAGTACGAGATGAAtt                                                              >  1:227064/4‑148 (MQ=255)
                                                                                tGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTg                                                             >  1:227060/1‑149 (MQ=255)
                                                                                tGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTg                                                             >  3:200005/1‑149 (MQ=255)
                                                                                tGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTg                                                             >  7:410191/1‑149 (MQ=255)
                                                                                  gCAGATTGGACGATTGCCGGTGGCGCACATGCGCGAGTTGATCTGGTTCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCa                                                           <  4:257924/149‑1 (MQ=255)
                                                                                   cAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGGTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAg                                                          <  3:240893/149‑1 (MQ=255)
                                                                                    aGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAga                                                         <  1:144404/149‑2 (MQ=255)
                                                                                               tGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGAc                                              >  8:372778/1‑149 (MQ=255)
                                                                                               tGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGAc                                              >  1:275844/1‑149 (MQ=255)
                                                                                                 gCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGc                                            >  1:266618/1‑149 (MQ=255)
                                                                                                 gCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATAGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGc                                            >  3:340583/1‑149 (MQ=255)
                                                                                                   cGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGCGAGGAATTGCAGGGCGTCGGCGACGctc                                          >  8:72625/1‑149 (MQ=255)
                                                                                                   cGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGctc                                          >  8:78459/1‑149 (MQ=255)
                                                                                                    ggTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGctct                                         <  2:386649/149‑1 (MQ=255)
                                                                                                      actcgcgCATGTGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTg                                        >  4:118794/4‑148 (MQ=255)
                                                                                                       ggCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGaa                                      <  2:344802/149‑1 (MQ=255)
                                                                                                           cacaTGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGTGGTGGTAGTTGGATATTAATTGCAGGGCGTCGGCGACGCTCTGAATCAg                                  <  2:243280/149‑1 (MQ=255)
                                                                                                            acaTGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAgg                                 >  2:268178/1‑149 (MQ=255)
                                                                                                             caTGCGCGAGTTGATCAGGATCTGGGCGATGGACTCGCGGGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGAGGAAATCGACGGGGGGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGgt                                <  7:72625/149‑1 (MQ=255)
                                                                                                             caTGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGgt                                <  7:78459/149‑1 (MQ=255)
                                                                                                                 cgcgAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTTCGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGgtcgtc                            <  1:309861/149‑1 (MQ=255)
                                                                                                                 cgcgAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGgtcgtc                            <  4:100748/149‑1 (MQ=255)
                                                                                                                     aGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTcttgc                        <  1:268178/149‑1 (MQ=255)
                                                                                                                       ttGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTcttgctt                      <  5:251095/149‑1 (MQ=255)
                                                                                                                        tGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTcttgctt                      >  7:133946/1‑148 (MQ=255)
                                                                                                                         gATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTcttgcttg                     <  8:148507/148‑1 (MQ=255)
                                                                                                                            cAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCa                 >  2:240403/1‑149 (MQ=255)
                                                                                                                              ggATCTGGTCGATTGAATCGCGTTCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACg               <  6:84891/149‑1 (MQ=255)
                                                                                                                              ggATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCATATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGGAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACg               <  4:63338/149‑1 (MQ=255)
                                                                                                                               gATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACgg              >  4:279781/1‑149 (MQ=255)
                                                                                                                                     ggCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATg         <  7:372778/148‑1 (MQ=255)
                                                                                                                                     ggCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGAAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTATTACGAGAAGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATg         <  5:324190/148‑1 (MQ=255)
                                                                                                                                      gCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATAGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCa       <  4:340583/149‑1 (MQ=255)
                                                                                                                                        tgaggAATCGCCTGCTGCACGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAgc     >  7:339598/4‑149 (MQ=255)
                                                                                                                                         gaggAATCGCCTGCTGCACGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAgcg    >  8:11874/3‑149 (MQ=255)
                                                                                                                                         gaggAATCGCCTGCTGCACGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAgcg    >  7:218373/3‑149 (MQ=255)
                                                                                                                                         gaggAATCGCCTGCTGCACGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAgcg    >  5:69579/3‑149 (MQ=255)
                                                                                                                                          tGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAgcgc   >  3:332690/1‑149 (MQ=255)
                                                                                                                                           ggAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAgcgcg  >  3:248105/1‑149 (MQ=255)
                                                                                                                                           ggAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGCCGTCGGCGACGATCTGAATCAGGTCTACTTGCTTGAACACGGACATGCAgcgcg  >  2:246382/1‑149 (MQ=255)
                                                                                                                                            gAATCGCCTGCTGCACGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCg                                                   >  2:88226/1‑99 (MQ=255)
                                                                                                                                            gAATCGCCTGCTGCACGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCg                                                   <  1:88226/99‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
TGTCGTCGACGCTCATGGTGGCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCG  >  NC_002947/1039115‑1039402

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