Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,039,262 G→C 5.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 5.3% 308.7 / 6.3 95R45G (CGC→GGC) PP_0897class 1 fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base G (43/47);  new base C (2/3);  total (45/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.49e-01

TGTCGTCGACGCTCATGGTGGCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGC  >  NC_002947/1039115‑1039403
                                                                                                                                                   |                                                                                                                                              
tGTCGTCGACGCTCATGGTGGCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGt                                                                                                                                               <  3:205484/149‑1 (MQ=255)
 gtcgtcGACGCTCATGGTGGCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTg                                                                                                                                              >  5:56288/1‑149 (MQ=255)
 gtcgtcGACGCTCATGGTGGCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTg                                                                                                                                              >  7:121007/1‑149 (MQ=255)
  tcgtcgACGCTCATGGTGGCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTg                                                                                                                                              <  3:116270/148‑1 (MQ=255)
   cgtcgACGCTCATGGTGGCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTg                                                                                                                                              >  1:52232/1‑147 (MQ=255)
    gtcgACGCTCATGGTGGCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTgca                                                                                                                                            <  3:9519/148‑1 (MQ=255)
    gtcgACGCTCATGGTGCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGGgcagc                                                                                                                                          >  4:224572/1‑149 (MQ=255)
        aCGCTCATGGTGGCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAgg                                                                                                                                       <  8:59667/149‑1 (MQ=255)
          gCTCATGGTGGCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGc                                                                                                                                      >  1:334774/1‑148 (MQ=255)
           cTCATGGTGGCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGa                                                                                                                                    <  4:53854/149‑1 (MQ=255)
            tCATGGTGGCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGAt                                                                                                                                   >  7:156466/1‑149 (MQ=255)
               tggtggCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCCTGCTGCACGCGATTcc                                                                                                                                <  7:327881/149‑1 (MQ=255)
                ggtggCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGGACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCt                                                                                                                              >  6:14851/1‑149 (MQ=255)
                ggtggCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCt                                                                                                                              >  1:320116/1‑149 (MQ=255)
                 gtggCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTc                                                                                                                             <  3:311473/149‑1 (MQ=255)
                  tggCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCa                                                                                                                            >  2:353861/1‑149 (MQ=255)
                  tggCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCa                                                                                                                            <  5:364689/149‑1 (MQ=255)
                  tggCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCa                                                                                                                            <  6:22425/149‑1 (MQ=255)
                     cgcCATCCCAGCGCACGTCCATGCCGATGGGGATGAACTCGGGGTGGTTACGGGTGTGCTGGCAGATGGGACGATTGCAGGTGGCGCAGATGCTCGAGTTGATCTGGATCTGGGCGTTGGATTCGCGTGCAGCAGGCGATTCCTCAGGc                                                                                                                         <  4:288407/149‑1 (MQ=255)
                      gcCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGt                                                                                                                                               >  7:123742/1‑127 (MQ=255)
                      gcCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGt                                                                                                                                               <  8:123742/127‑1 (MQ=255)
                        cATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACg                                                                                                                          >  3:321512/1‑145 (MQ=255)
                         aTCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGAt                                                                                                                     <  2:320116/149‑1 (MQ=255)
                          tCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATa                                                                                                                    <  3:136909/149‑1 (MQ=255)
                                     gTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCa                                                                                                         <  4:68655/149‑1 (MQ=255)
                                     gTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCa                                                                                                         >  4:359319/1‑149 (MQ=255)
                                        cATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATgg                                                                                                      >  6:285617/1‑149 (MQ=255)
                                          tGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGcc                                                                                                    <  1:353861/149‑1 (MQ=255)
                                             ccACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTg                                                                                                  >  2:158174/1‑148 (MQ=255)
                                              cACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGa                                                                                                >  6:378751/1‑149 (MQ=255)
                                              cACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGa                                                                                                >  8:170848/1‑149 (MQ=255)
                                                   ggATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGaaa                                                                                           <  4:321512/149‑1 (MQ=255)
                                                   ggATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGaaa                                                                                           <  2:52232/149‑1 (MQ=255)
                                                         acacGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGa                                                                                                >  4:84439/1‑138 (MQ=255)
                                                         acacGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGa                                                                                                <  3:84439/138‑1 (MQ=255)
                                                         acacGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGAt                                                                                               <  3:324279/139‑1 (MQ=255)
                                                         acacGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGAt                                                                                               >  4:324279/1‑139 (MQ=255)
                                                         acacGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACg                                                                                     <  2:334774/149‑1 (MQ=255)
                                                         acaaGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACg                                                                                     >  4:8310/1‑149 (MQ=255)
                                                            cGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACgggg                                                                                  <  4:69924/149‑1 (MQ=255)
                                                             ggTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGggt                                                                                 >  8:310034/1‑149 (MQ=255)
                                                              gTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCa                                                                                                                        <  2:92999/109‑1 (MQ=255)
                                                              gTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCa                                                                                                                        >  1:92999/1‑109 (MQ=255)
                                                                 aCGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGgta                                                                             <  5:14851/149‑1 (MQ=255)
                                                                   gATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGgtag                                                                            >  6:1030/1‑148 (MQ=255)
                                                                       ccGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACga                                                                       <  8:156466/149‑1 (MQ=255)
                                                                       ccGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACga                                                                       <  7:310034/149‑1 (MQ=255)
                                                                        cGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTAc                                                                         >  6:259743/1‑146 (MQ=255)
                                                                         ggTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACgaga                                                                     <  3:8310/149‑1 (MQ=255)
                                                                          gtgtCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGAt                                                                    <  5:378751/149‑1 (MQ=255)
                                                                           tgtCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATg                                                                   <  6:56288/149‑1 (MQ=255)
                                                                             tCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGaa                                                                 >  7:60582/1‑149 (MQ=255)
                                                                               cTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAAtt                                                               >  5:173317/1‑149 (MQ=255)
                                                                                  gCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTg                                                              >  6:167260/1‑147 (MQ=255)
                                                                                    aGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAgg                                                          <  1:158174/149‑1 (MQ=255)
                                                                                      aTTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACg                                                                        <  6:194438/133‑1 (MQ=255)
                                                                                      aTTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACg                                                                        >  5:194438/1‑133 (MQ=255)
                                                                                      aTTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGc                                                        <  3:359319/149‑1 (MQ=255)
                                                                                         ggACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTc                                                     >  8:210320/1‑149 (MQ=255)
                                                                                           aCGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCgg                                                   >  6:110028/1‑149 (MQ=255)
                                                                                               tGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGAc                                               >  7:194716/1‑149 (MQ=255)
                                                                                                   cGGTGGCGCACATGCGCGAGTTGCTCAGGATCGGGGCGATGGACTCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGAGCAAATCGAAGGGGTGGTAGTCCGAGATGAATTGCAGGGCGTCTGCGACGctc                                           >  6:365244/1‑149 (MQ=255)
                                                                                                   cGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGctc                                           <  5:110028/149‑1 (MQ=255)
                                                                                                          gcACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCATGGCGTCGGCGACGCTCTGCAt                                      >  4:34939/1‑147 (MQ=255)
                                                                                                          gcACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGctct                                          >  7:34958/1‑143 (MQ=255)
                                                                                                          gcACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGctct                                          <  8:34958/143‑1 (MQ=255)
                                                                                                          gcACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAAt                                      >  2:141633/1‑147 (MQ=255)
                                                                                                           cacaTGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAg                                   <  7:210320/149‑1 (MQ=255)
                                                                                                              aTGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGTGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGgtc                                <  5:259743/149‑1 (MQ=255)
                                                                                                              aTGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGgtag                                                                            >  7:281507/1‑105 (MQ=255)
                                                                                                              aTGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGgtag                                                                            <  8:281507/105‑1 (MQ=255)
                                                                                                               tGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGgtcg                               >  1:154434/1‑149 (MQ=255)
                                                                                                                 cgcgAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGgtcgtc                             <  2:154434/149‑1 (MQ=255)
                                                                                                                  gcgAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACgaga                                                                     >  7:105255/1‑108 (MQ=255)
                                                                                                                  gcgAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACgaga                                                                     <  8:105255/108‑1 (MQ=255)
                                                                                                                     aGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTcttgc                         <  5:167260/149‑1 (MQ=255)
                                                                                                                      gTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTcttgct                        <  6:368623/149‑1 (MQ=255)
                                                                                                                        tGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTcttgctt                       <  5:285617/148‑1 (MQ=255)
                                                                                                                             aGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAg                                   >  7:117799/1‑131 (MQ=255)
                                                                                                                             aGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAg                                   <  8:117799/131‑1 (MQ=255)
                                                                                                                             aGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCa                  <  7:170848/148‑1 (MQ=255)
                                                                                                                             aGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCAc                 >  3:84373/1‑149 (MQ=255)
                                                                                                                               gATCTGGGCGATGGAATCGCGTGCAGCTGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACgg               >  6:253471/1‑149 (MQ=255)
                                                                                                                                 tCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACgg               >  1:113924/1‑147 (MQ=255)
                                                                                                                                 tCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACgg               <  2:113924/147‑1 (MQ=255)
                                                                                                                                     tgctttGGAATCGCGTGCTGCAGTCGATTACTCACTCAGATAGTCCTAGTGCATGTCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGCATTGCCGGGCGTCGGCGACGCTCTGAATCTGGGCGTCTTGCTTGATCACGGt              <  5:365244/139‑1 (MQ=255)
                                                                                                                                     ggCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATg          <  8:194716/148‑1 (MQ=255)
                                                                                                                                     gcgtgaggAATCGCCTGCTGCACGCGATTCCATCGCCCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTcttgctt                       >  7:52724/7‑136 (MQ=37)
                                                                                                                                     gcgtgaggAATCGCCTGCTGCACGCGATTCCATCGCCCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTcttgctt                       <  8:52724/130‑1 (MQ=37)
                                                                                                                                      cgtgaggAATCGCCTGCTGCACGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTg                                         <  4:24897/111‑1 (MQ=39)
                                                                                                                                        gATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAgc      <  1:141633/149‑1 (MQ=255)
                                                                                                                                         gaggAATCGCCTGCTGCACGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAgcg     >  3:154475/3‑149 (MQ=255)
                                                                                                                                           ggAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAgcgcg   >  7:66205/1‑149 (MQ=255)
                                                                                                                                            gAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAgcgcgc  <  4:84373/149‑1 (MQ=255)
                                                                                                                                               tCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGgtcgtc                             >  7:149273/1‑119 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
TGTCGTCGACGCTCATGGTGGCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGC  >  NC_002947/1039115‑1039403

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: