Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,151,972 C→T 6.0% Q231* (CAG→TAG) ‡ edd → phosphogluconate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,151,9720CT 6.0% 250.4 / 4.7 83Q231* (CAG→TAG) ‡eddphosphogluconate dehydratase
Reads supporting (aligned to +/- strand):  ref base C (40/37);  new base T (2/3);  total (43/40)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

TGGTTTCCGGCATCTCCAACAAGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCATGCCGCTGGGC  >  NC_002947/1151823‑1152118
                                                                                                                                                     |                                                                                                                                                  
tGGTTTCCGGCATCTCCAACAAGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAATCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACaccaacc                                                                                                                                                    >  2:629990/1‑150 (MQ=255)
   tttCCGGCATCTCCAACAAGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGt                                                                                                                                                 >  1:748219/1‑150 (MQ=255)
   tttCCGGCATCTCCAACAAGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGt                                                                                                                                                 <  1:37462/150‑1 (MQ=255)
     tCCGGCATCTCCAACAAGCAGAAGGCCGATGTGCGCCAGCGCTATGCTGCAGGCAAGGCCAGCCGCGAAGAGATGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTg                                                                                                                                               <  1:1343647/150‑1 (MQ=255)
       cGGCATCTCCAACAAGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCATCACCAACtggttggt                                                                                                                                             >  2:50283/1‑150 (MQ=255)
         gCATCTCCAACAAGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGa                                                                                                                                           <  2:253862/150‑1 (MQ=255)
         gCATCTCCAACAAGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGa                                                                                                                                           >  1:253862/1‑150 (MQ=255)
           aTCTCCAACAAGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATg                                                                                                                                         <  2:99028/150‑1 (MQ=255)
            tctcCAACAAGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCATCACCAACTGGTTGGTGTTgg                                                                                                                                        <  1:50283/150‑1 (MQ=255)
            tctcCAACAAGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACAGCCAACACCAACCAGTTGGTGATgg                                                                                                                                        >  1:187684/1‑150 (MQ=255)
               ccAACAAGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAg                                                                                                                                     >  2:145580/1‑150 (MQ=255)
                caacaaGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCATCACCAACTGGTTGGTGTTgg                                                                                                                                        <  2:1491530/146‑1 (MQ=255)
                caacaaGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCATCACCAACTGGTTGGTGTTgg                                                                                                                                        >  1:1491530/1‑146 (MQ=255)
                    aaGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGTCGCCAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATg                                                                                                                                <  1:100412/150‑1 (MQ=255)
                    aaGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATg                                                                                                                                <  2:78077/150‑1 (MQ=255)
                      gCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATggg                                                                                                                              <  1:145580/150‑1 (MQ=255)
                         gAAGGCCGCTGTCCGCCAGCGCTATGCTGAAGGCAAGGCCAGGCGCGAAGAGCTGTTGGAGTACGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCt                                                                                                                           >  2:811033/1‑150 (MQ=255)
                          aaGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTg                                                                                                                          <  1:1457744/150‑1 (MQ=255)
                           aGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGc                                                                                                                         >  2:247254/1‑150 (MQ=255)
                           aGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGc                                                                                                                         <  2:1573401/150‑1 (MQ=255)
                                gATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGa                                                                                                                                           >  2:1454235/1‑127 (MQ=255)
                                gATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGa                                                                                                                                           <  1:1454235/127‑1 (MQ=255)
                                   gtgCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCg                                                                                                                 <  2:870014/150‑1 (MQ=255)
                                   gtgCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACAAGTTGGGGGTGGAAGTGATGGGGCTGCACCTGCCg                                                                                                                 >  2:1440426/1‑150 (MQ=255)
                                        ccAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGa                                                                                                                                           <  2:1042413/119‑1 (MQ=255)
                                        ccAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGa                                                                                                                                           >  1:1042413/1‑119 (MQ=255)
                                          aGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATg                                                                                                                                >  2:744171/1‑128 (MQ=255)
                                          aGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATg                                                                                                                                <  1:744171/128‑1 (MQ=255)
                                          aGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCt                                                                                                          >  1:1490764/1‑150 (MQ=255)
                                               tATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCtcgttc                                                                                                     <  1:247254/150‑1 (MQ=255)
                                                 tGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCtcgttcgt                                                                                                   >  2:522419/1‑150 (MQ=255)
                                                    cGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCaa                                                                                                >  1:1629858/1‑150 (MQ=255)
                                                     gAAGGCAAGGCCAGCCGAGAAGGGCGGTTGGAGTGCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAAc                                                                                               <  1:360409/150‑1 (MQ=255)
                                                      aaggcaaggcCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAAcc                                                                                              <  2:193186/150‑1 (MQ=255)
                                                        ggcaaggcCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATggg                                                                                                                              <  2:586962/116‑1 (MQ=255)
                                                        ggcaaggcCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATggg                                                                                                                              >  1:586962/1‑116 (MQ=255)
                                                           aaggcCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCATCACCAACTGGTTGGTGTTggcg                                                                                                                                      <  1:1148276/105‑3 (MQ=255)
                                                           aaggcCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCATCACCAACTGGTTGGTGTTggcg                                                                                                                                      >  2:1148276/1‑103 (MQ=255)
                                                           aaggcCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTac                                                                                         >  1:171071/1‑150 (MQ=255)
                                                               ccAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACAcgc                                                                                     <  2:1490764/150‑1 (MQ=255)
                                                                cAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCtcgttc                                                                                                     >  1:1407096/1‑133 (MQ=255)
                                                                cAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCtcgttc                                                                                                     >  1:1406026/1‑133 (MQ=255)
                                                                cAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCtcgttc                                                                                                     <  2:1407096/133‑1 (MQ=255)
                                                                cAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCtcgttc                                                                                                     <  2:1406026/133‑1 (MQ=255)
                                                                   ccGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAg                                                                                                                                     <  1:57233/98‑1 (MQ=255)
                                                                   ccGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAg                                                                                                                                     >  2:57233/1‑98 (MQ=255)
                                                                    cgcgAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGAt                                                                                                                                          <  1:1237995/92‑1 (MQ=255)
                                                                    cgcgAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGAt                                                                                                                                          >  2:1237995/1‑92 (MQ=255)
                                                                       gAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACAGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCCCCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCga                                                                             <  2:187684/150‑2 (MQ=255)
                                                                           agCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCAGAGCGCCACAGTCGTCAACCCGTACACGCCGCTGCGTGACg                                                                         >  2:543488/1‑150 (MQ=255)
                                                                                  ggAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCAc                                                                  >  2:1150028/1‑150 (MQ=255)
                                                                                      tCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCCGTTGGTAACGGAATTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCAccgcc                                                              <  1:1433332/150‑1 (MQ=255)
                                                                                         gagaTGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATgg                                                                                                                               >  2:689320/1‑82 (MQ=255)
                                                                                         gagaTGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATgg                                                                                                                               <  1:689320/82‑1 (MQ=255)
                                                                                         gagaTGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCAccgccg                                                             >  1:558408/1‑148 (MQ=255)
                                                                                         gagaTGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAg                                                           >  2:1352458/1‑150 (MQ=255)
                                                                                         gagaTGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTAAACGCCGCTTCGTGACGCACTCAccgccg                                                             <  2:558408/148‑1 (MQ=255)
                                                                                          agaTGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGaa                                                                                                                                      <  1:1516255/74‑1 (MQ=255)
                                                                                          agaTGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGaa                                                                                                                                      >  2:1516255/1‑74 (MQ=255)
                                                                                                aTTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTgcgc                                                    <  2:748219/150‑1 (MQ=255)
                                                                                                     tACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGcagcag                                               >  1:1257724/1‑150 (MQ=255)
                                                                                                       ccACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGt                                             >  1:188464/1‑150 (MQ=255)
                                                                                                          caGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGAt                                                                                                                                          >  2:762717/1‑54 (MQ=255)
                                                                                                          caGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGAt                                                                                                                                          <  1:762717/54‑1 (MQ=255)
                                                                                                          caGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGAc                                          >  1:1226239/1‑150 (MQ=255)
                                                                                                          caGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGAc                                          >  2:1029916/1‑150 (MQ=255)
                                                                                                               ccGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGNTGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCa                                     >  1:719142/1‑150 (MQ=255)
                                                                                                                     cctgcacctTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCa                               >  2:484267/1‑150 (MQ=255)
                                                                                                                      ctgcacctTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCaa                              >  2:1180166/1‑150 (MQ=255)
                                                                                                                      ctgcacctTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCaa                              <  1:1273869/150‑1 (MQ=255)
                                                                                                                       tgcacctTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGc                                                                               <  2:931940/100‑1 (MQ=255)
                                                                                                                       tgcacctTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGc                                                                               >  1:931940/1‑100 (MQ=255)
                                                                                                                       tgcacctTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCGCTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAg                             <  2:1629858/150‑1 (MQ=255)
                                                                                                                        gcacctTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAgg                            <  1:1029916/150‑1 (MQ=255)
                                                                                                                          acctTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGcc                          >  2:1348074/1‑150 (MQ=255)
                                                                                                                            ctTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAAcc                                                                                              >  1:1286518/1‑80 (MQ=255)
                                                                                                                            ctTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAAcc                                                                                              <  2:1286518/80‑1 (MQ=255)
                                                                                                                               cTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCaa                              <  2:194073/141‑1 (MQ=255)
                                                                                                                               cTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCaa                              >  1:194073/1‑141 (MQ=255)
                                                                                                                                    gCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCt                >  2:994246/1‑150 (MQ=255)
                                                                                                                                    gCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCt                >  2:503854/1‑150 (MQ=255)
                                                                                                                                    gCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCt                <  1:629990/150‑1 (MQ=255)
                                                                                                                                        ttccAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCAt            >  2:92081/3‑150 (MQ=255)
                                                                                                                                        cGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCACGGCCAGCGGCAGCTTCAt            <  2:1257724/150‑1 (MQ=255)
                                                                                                                                            aaCACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCATGCCg        >  2:340452/1‑150 (MQ=255)
                                                                                                                                              ccccaaccaGTGGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCATGCCGCt      <  1:811033/148‑1 (MQ=255)
                                                                                                                                               accaaccaGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACAcgccgc                                                                                  <  1:1434945/73‑1 (MQ=255)
                                                                                                                                               accaaccaGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACAcgccgc                                                                                  >  2:1434945/1‑73 (MQ=255)
                                                                                                                                               accaaccaGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAg                             <  1:753824/126‑1 (MQ=255)
                                                                                                                                               accaaccaGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAg                             >  2:753824/1‑126 (MQ=255)
                                                                                                                                               accaaccaGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCATGCCGCTg     >  1:684906/1‑150 (MQ=255)
                                                                                                                                                ccaaccgGTTGGTAACTGGAGTGCGGGGTCTGCACCTGCCGGGCGCCTCGTGCGTCAGCCCGTACCCGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCATGCCGCTgg    <  2:719142/150‑1 (MQ=255)
                                                                                                                                                  aacctGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGCCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCTGGGCCAGCGGCAGCTTCATGGCGATGGGc  <  2:405302/150‑1 (MQ=255)
                                                                                                                                                  aaccaGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCATGCCGCTGGGc  >  2:1240574/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TGGTTTCCGGCATCTCCAACAAGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCATGCCGCTGGGC  >  NC_002947/1151823‑1152118

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: