Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451 1,126,647 +C intergenic (‑101/+48) tdcG‑II ← / ← gcvP‑I L‑serine dehydratase/glycine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154511,126,6451.C97.0% 632.5 / 14.3 166intergenic (‑99/+50)tdcG‑II/gcvP‑IL‑serine dehydratase/glycine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (1/4);  new base C (87/73);  total (88/78)
Fisher's exact test for biased strand distribution p-value = 1.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 176 total aligned reads displayed.

AGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTAC  >  AE015451/1126495‑1126788
                                                                                                                                                        |                                                                                                                                               
aGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTaaa                                                                                                                                                  <  1:958929/151‑1 (MQ=255)
 gCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAc                                                                                                                                                <  1:958930/151‑1 (MQ=255)
   tGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAccc                                                                                                                                              <  1:958936/151‑2 (MQ=255)
   tGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAc                                                                                                                                                <  1:958934/149‑1 (MQ=255)
       ggggCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAcccgctc                                                                                                                                          <  1:958939/151‑6 (MQ=255)
         ggCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCt                                                                                                                                         <  1:958942/151‑1 (MQ=255)
               tgccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTAcaca                                                                                                                                    <  2:958944/150‑1 (MQ=255)
                 ccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACagag                                                                                                                                 <  1:958945/151‑1 (MQ=255)
                        tGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCa                                                                                                                          <  2:958947/151‑1 (MQ=255)
                          aaCAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCAcc                                                                                                                        >  2:959212/1‑151 (MQ=255)
                            cAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAcccc                                                                                                                                              <  2:958949/126‑2 (MQ=255)
                             aGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAACCCCGCTCCTACACAGAGCCCGCCACCCaa                                                                                                                     <  1:958950/151‑1 (MQ=255)
                              ggTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCaa                                                                                                                     <  2:958954/150‑1 (MQ=255)
                              ggTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCCCCTACACAGAGCCCGCCACCCaa                                                                                                                     >  1:959074/1‑150 (MQ=255)
                                     gACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGt                                                                                                              <  2:958957/150‑1 (MQ=255)
                                     gACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTg                                                                                                             <  1:958958/151‑1 (MQ=255)
                                      aCGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTgg                                                                                                            <  1:958959/151‑1 (MQ=255)
                                       cGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAgg                                                                                                               >  2:959238/1‑147 (MQ=255)
                                        gCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGGATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCg                                                                                                          >  2:959066/1‑151 (MQ=255)
                                         cTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCg                                                                                                          <  2:958961/150‑1 (MQ=255)
                                          tCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCgg                                                                                                         <  2:958962/150‑1 (MQ=255)
                                             gTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAg                                                                                                     <  1:958964/151‑1 (MQ=255)
                                              tGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAgg                                                                                                    >  2:959152/1‑151 (MQ=255)
                                              tGACATGTGGACCTCCTGCCAGGGTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGCTATAGCCACCATCGCAGTCGCGGGTAAACCCGCTCCTACGCAGAGCCCGCCACCCAAGCCAGGTGGCg                                                                                                          >  2:959137/1‑145 (MQ=38)
                                               gACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAggg                                                                                                   <  2:958966/151‑1 (MQ=255)
                                                     tgGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGCGGGCTATc                                                                                              >  1:959147/1‑150 (MQ=255)
                                                     tgGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCCAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCa                                                                                             <  1:958969/151‑1 (MQ=255)
                                                         cctcctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCg                                                                                         <  2:958971/151‑1 (MQ=255)
                                                          ctcctGCCGGGTTGTTGGGTGGACGTTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGGGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACaga                                                                                                                                  <  2:958973/109‑1 (MQ=25)
                                                          ctcctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGGTTACCCGCGAATGCGGAGGGATAGCCACCATCGCATTCGCGGGGAAACCCGCTCCTac                                                                                                                                       >  2:959099/1‑104 (MQ=255)
                                                           tcctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCATACACAGAGCCCGCCACCCAAGCCAGGTGGCgggg                                                                                                       >  2:959289/1‑135 (MQ=255)
                                                            cctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCCGGTGGCGGGGAGGGCTATCAGGCGTcc                                                                                      >  2:959183/1‑151 (MQ=255)
                                                             ctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCt                                                                                     >  1:959149/1‑151 (MQ=255)
                                                             ctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCt                                                                                     >  1:959112/1‑151 (MQ=255)
                                                               gCCAGGTTGTTGGGGGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGa                                                                                   <  1:958975/151‑1 (MQ=255)
                                                                   ggTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGc                                                                               <  1:958976/151‑1 (MQ=255)
                                                                      tgttgGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGctc                                                                             >  1:959339/1‑150 (MQ=255)
                                                                       gttgGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCtctc                                                                           <  1:958978/151‑1 (MQ=255)
                                                                       gttgGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCtctc                                                                           <  1:958980/151‑1 (MQ=255)
                                                                        ttgGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGAGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAACCCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCtcg                                                                          <  2:958982/151‑1 (MQ=255)
                                                                           ggTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCtcgatc                                                                       <  1:958983/151‑1 (MQ=255)
                                                                                   ttGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACaga                                                                                                                                  <  2:958986/84‑1 (MQ=255)
                                                                                   ttGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCa                                                               <  1:958984/151‑1 (MQ=255)
                                                                                    tGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCtcgatcg                                                                      >  1:959203/1‑143 (MQ=255)
                                                                                         ccTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGca                                                         <  2:958989/151‑1 (MQ=255)
                                                                                             tAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTaca                                                                                                                                      >  2:959269/1‑70 (MQ=255)
                                                                                             tAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACAcc                                                     >  1:959096/1‑151 (MQ=255)
                                                                                              aGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCa                                                    >  1:959237/1‑151 (MQ=255)
                                                                                              aGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGGCCTGATAGCTCTCGATCGACGGGCAGGCGCACAcc                                                     >  1:959285/1‑150 (MQ=255)
                                                                                               ggCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAg                                                   <  2:958992/151‑1 (MQ=255)
                                                                                               ggCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAg                                                   <  2:958991/151‑1 (MQ=255)
                                                                                                 cgcgGGTTTACCCGTGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCTGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGTGCAGGCGCACACCAGGt                                                 <  1:958994/151‑1 (MQ=255)
                                                                                                  gcgGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGtt                                                >  2:959199/1‑151 (MQ=255)
                                                                                                  gcgGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGtt                                                >  1:959182/1‑151 (MQ=255)
                                                                                                    gggTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGc                                              >  1:959254/1‑151 (MQ=255)
                                                                                                       tttACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCgg                                            >  1:959211/1‑150 (MQ=255)
                                                                                                        ttACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTc                                          >  2:959282/1‑151 (MQ=255)
                                                                                                        ttACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTc                                          <  2:958996/151‑1 (MQ=255)
                                                                                                         tACCCGCGAATGCGGAGGTATAGCCACCATCGCATGCGCGGGTAAGCCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCtcgatc                                                                       <  1:958998/121‑1 (MQ=255)
                                                                                                            ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACgg                                                                  >  2:959302/1‑123 (MQ=255)
                                                                                                            ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCg                                      >  2:959185/1‑151 (MQ=255)
                                                                                                            ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCg                                      >  2:959215/1‑151 (MQ=255)
                                                                                                             cgcgAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGgcg                                                           <  2:959000/129‑1 (MQ=255)
                                                                                                              gcgAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGctc                                                                             <  2:959002/110‑1 (MQ=255)
                                                                                                              gcgAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGctc                                                                             >  1:959002/1‑110 (MQ=255)
                                                                                                              gcgAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGaa                                    >  1:959222/1‑151 (MQ=255)
                                                                                                              gcgAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCGCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGa                                     >  2:959301/1‑150 (MQ=255)
                                                                                                               cgAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGt                                           <  1:959003/143‑1 (MQ=255)
                                                                                                               cgAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAac                                   >  1:959398/1‑151 (MQ=255)
                                                                                                                gAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAaca                                  >  1:959414/1‑151 (MQ=255)
                                                                                                                gAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAaca                                  <  2:959006/151‑1 (MQ=255)
                                                                                                                gAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAaca                                  <  1:959005/151‑1 (MQ=255)
                                                                                                                gAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAaca                                  >  1:959125/1‑151 (MQ=255)
                                                                                                                 aaTGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAacac                                 <  1:959008/151‑1 (MQ=255)
                                                                                                                  aTGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACg                                >  2:959249/1‑151 (MQ=255)
                                                                                                                   tGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACg                                >  2:959477/1‑150 (MQ=255)
                                                                                                                    gCGGAGGTATAGCCACCATCGCATTCGCGGGTAGACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGt                               <  2:959011/150‑1 (MQ=255)
                                                                                                                     cGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGtt                              >  2:959162/1‑150 (MQ=255)
                                                                                                                      ggAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACg                                                                   <  2:959014/112‑1 (MQ=255)
                                                                                                                      ggAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGca                                                         <  1:959013/122‑1 (MQ=255)
                                                                                                                      ggAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGt                            >  1:959242/1‑151 (MQ=255)
                                                                                                                      ggAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGt                            >  2:959397/1‑151 (MQ=255)
                                                                                                                         ggTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAg                                                                                >  1:959423/1‑96 (MQ=255)
                                                                                                                         ggTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCg                          >  2:959273/1‑150 (MQ=255)
                                                                                                                            ataGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAccc                      >  2:959352/1‑151 (MQ=255)
                                                                                                                            ataGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAccc                      <  1:959017/151‑1 (MQ=255)
                                                                                                                            ataGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAccc                      >  2:959119/1‑151 (MQ=255)
                                                                                                                             taGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGcc                                                                                                                               <  2:959018/45‑1 (MQ=14)
                                                                                                                              aGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCgg                    >  1:959382/1‑151 (MQ=255)
                                                                                                                               gCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGc                   >  1:959214/1‑151 (MQ=255)
                                                                                                                                ccaccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGcc                  >  1:959228/1‑151 (MQ=255)
                                                                                                                                 caccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCg                 >  1:959262/1‑151 (MQ=255)
                                                                                                                                   ccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGAc               <  1:959019/151‑1 (MQ=255)
                                                                                                                                         cATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggcgg         >  1:959265/1‑151 (MQ=255)
                                                                                                                                           ttCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCggcggc        <  1:959021/150‑1 (MQ=255)
                                                                                                                                           ttCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGcc       >  1:959239/1‑151 (MQ=255)
                                                                                                                                               cgGGTACACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTg                             <  2:959023/125‑1 (MQ=255)
                                                                                                                                                gggTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTAc  >  1:959106/1‑151 (MQ=255)
                                                                                                                                                gggTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTAc  >  2:959196/1‑151 (MQ=255)
                                                                                                                                                gggTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTAc  >  2:959139/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                               
AGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTAC  >  AE015451/1126495‑1126788

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 43 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: