Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 5,822,730 T→C 5.7% S87P (TCG→CCG)  rdgB → dITP/XTP pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154515,822,7300TC 5.7% 238.1 / 3.4 70S87P (TCG→CCG) rdgBdITP/XTP pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base T (22/44);  new base C (2/2);  total (24/46)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01

GCCAGCCACAATGCCGGCAAACTGAAGGAACTCCAGGCCATGCTCGGCCATTCCGTGCAGCTGCACTCGATTGGTGAGTTCAGCCAGGTGGAGCCGGAAGAAACCGGCCTGTCGTTTGTCGAGAACGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGCTTTGGCTACGATCCGCTGT  >  AE015451/5822499‑5822952
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                              
gccagccaCAATGCCGGCAAACTGAAGGAACTCCAGGCCATGCTCGGCCATTCCGTGCAGCTGCACTCGATTGGTGAGTTCAGCCAGGTGGAGCCGGAAGAAACCGGCCTGTCGTTTGTCGAGAACGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCg                                                                                                                                                                                                                              <  1:9509/234‑1 (MQ=255)
                  aaaCTGAAGGAACTCCAGGCCATGCTCGGCCATTCCGTGCAGCTGCACTCGATTGGTGAGTTCAGCCAGGTGGAGCCGGAAGAAACCGGCCTGTCGTTTGTCGAGAACGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCaa                                                                                                                                                                                                          >  1:763484/1‑236 (MQ=255)
                                         gCTCGGCCATTCCGTGCAGCTGCAATCGATTGGTGAGTTCAGCCAGGTGGAGCCGGAAGAAACCGGCCTGTCGTTTGTCGAGAACGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGcc                                                                                                                                                                                    <  1:75476/235‑1 (MQ=255)
                                            cGGCCATTCCGTGCAGCTGCACTCGATTGGTGAGTTCAGCCAGGTGGAGCCGGAAGAAACCGGCCTGTCGTTTGTCGAGAACGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAg                                                                                                                                                                                 >  1:422901/1‑235 (MQ=255)
                                                             tGCACTCGATTGGTGAGTTCAGCCAGGTGGAGCCGGAAGAAACCGGCCTGTCGTTTGTCGAGAACGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGaaa                                                                                                                                                               <  1:681914/236‑1 (MQ=255)
                                                             tGCACTCGATTGGTGAGTTCAGCCAGGTGGAGCCGGAAGAAACCGGCCTGTCGTTTGTCGAGAACGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGaaa                                                                                                                                                               <  1:454311/236‑1 (MQ=255)
                                                                aCTCGATTGGTGAGTTCAGCCAGGTGGAGCCGGAAGAAACCGGCCTGTCGTTTGTCGAGAACGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTACTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCACAGGTGACGCTGCAAACAACGCCAAGCTGCTCGAAGCCTTGaaa                                                                                                                                                               <  2:763481/233‑1 (MQ=255)
                                                                          tGAGTTCAGCCAGGTGGAGCCGGAAGAAACCGGCCTGTCGTTTGTCGAGAACGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTAt                                                                                                                                                                                                                     <  1:401251/169‑1 (MQ=255)
                                                                             gTTCAGCCAGGTGGAGCCGGAAGAAACCGGCCTGTCGTTTGTCGAGAACGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTg                                                                                                                                               <  1:524323/236‑1 (MQ=255)
                                                                                cAGCCAGGTGGAGCCGGAAGAAACCGGCCTGTCGTTTGTCGAGAACGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTTGAAGCCTTGAAAGACGTGCCTGAAGCTGAgc                                                                                                                                            <  1:410665/236‑1 (MQ=255)
                                                                                    cAGGTGGAGCCGGAAGAAACCGGCCTGTCGTTTGTCGAGAACGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGa                                                                                                                                              <  1:1145276/230‑1 (MQ=255)
                                                                                          gAGCCGGAAGAAACCGGCCTGTCGTTTGTCGAGAACGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGAcgcg                                                                                                                                                                                                >  1:815740/1‑174 (MQ=255)
                                                                                          gAGCCGGAAGAAACCGGCCTGTCGTTTGTCGAGAACGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACGCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGa                                                                                                                                              <  1:41345/224‑1 (MQ=255)
                                                                                           aGCCGGAAGAAACCGGCCTGTCGTTTGTCGAGAACGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGAc                                                                                                                                                                                                   >  2:708375/1‑170 (MQ=255)
                                                                                           aGCCGGAAGAAACCGGCCTGTCGTTTGTCGAGAACGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGAc                                                                                                                                                                                                   <  1:708378/170‑1 (MQ=255)
                                                                                              cGGAAGAAACCGGCCTGTCGTTTGTCGAGAACGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGa                                                                                                                                              <  1:832306/220‑1 (MQ=255)
                                                                                                    aaaCCGGCCTGTCGTTTGTCGAGAACGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGaaa                                                                                                                                                               >  1:1111040/1‑197 (MQ=255)
                                                                                                     aaCCGGCCTGTCGTTTGTCGAGAACGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGt                                                                                                                            <  1:868807/231‑1 (MQ=255)
                                                                                                       ccGGCCTGTCGTTTGTCGAGAACGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGGCTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGaaa                                                                                                                                                               <  1:277650/194‑1 (MQ=255)
                                                                                                                 gTTTGTCGAGAACGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGAGGCAAACAACGCCAAGCTGCTCGa                                                                                                                                                                         <  1:461744/174‑1 (MQ=255)
                                                                                                                  tttGTCGAGAACGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTc                                                                                                                                                                           <  1:838223/171‑1 (MQ=255)
                                                                                                                  tttGTCGAGAACGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTAGGCCGATGGCAAAGGTGACGCGGCAAACAATGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTgg                                                                                                           >  1:676003/1‑235 (MQ=255)
                                                                                                                     gTCGAGAACGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGaaa                                                                                                                                                               <  2:572838/180‑1 (MQ=255)
                                                                                                                      tCGAGAACGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTgcg                                                                                                       <  1:1161572/235‑1 (MQ=255)
                                                                                                                            aCGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGa                                                                                                                                              <  2:1154868/190‑1 (MQ=255)
                                                                                                                            aCGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGa                                                                                                                                              >  1:1154871/1‑190 (MQ=255)
                                                                                                                              gCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGt                                                                                                                            <  2:782857/206‑1 (MQ=255)
                                                                                                                              gCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGt                                                                                                                            >  1:782860/1‑206 (MQ=255)
                                                                                                                                cATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCcgac                                                                                             <  1:1242716/235‑1 (MQ=255)
                                                                                                                                   ccTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCa                                                                                                                                                                                   <  1:506405/146‑1 (MQ=255)
                                                                                                                                   ccTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCa                                                                                                                                                                                   >  2:506402/1‑146 (MQ=255)
                                                                                                                                          gcgcgcAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCcgacga                                                                                           <  1:1281298/227‑1 (MQ=255)
                                                                                                                                           cgcgcAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTc                                                                                                                                                                                                                               <  1:1183821/94‑1 (MQ=255)
                                                                                                                                           cgcgcAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTc                                                                                                                                                                                                                               >  2:1183818/1‑94 (MQ=255)
                                                                                                                                               cAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGc                                                                                                                        >  1:42405/1‑193 (MQ=255)
                                                                                                                                                aaCGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATcc                                                                             >  1:654869/1‑235 (MQ=255)
                                                                                                                                                aaCGCCGCCCGCATTGCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTgc                                                                                                        <  1:75921/208‑1 (MQ=255)
                                                                                                                                                 aCGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATcc                                                                             >  2:630664/1‑234 (MQ=255)
                                                                                                                                                 aCGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGACTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCtg                                                                           >  1:676390/1‑236 (MQ=255)
                                                                                                                                                 aCGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGAGTAGATGCCCGGCGCGCTATGCCGATgg                                                                                                                                                                                                             <  1:663252/106‑1 (MQ=255)
                                                                                                                                                 aCGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGAGTAGATGCCCGGCGCGCTAAGCCGATgg                                                                                                                                                                                                             >  2:663249/1‑106 (MQ=255)
                                                                                                                                                            aTTTCCGGCCTGCCGGCCGTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCa                                                                                                                                         <  2:302321/163‑1 (MQ=255)
                                                                                                                                                                ccGGCCTGCCGGCCTTGGCCGATGACTCTGGACGGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGa                                                                                                                                              <  2:264523/154‑1 (MQ=255)
                                                                                                                                                                     cTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGc                                                                                                                                                                                                                   <  2:165133/80‑1 (MQ=255)
                                                                                                                                                                          ggCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCtgtg                                                             <  1:116517/225‑1 (MQ=255)
                                                                                                                                                                           gCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTgg                                                                                                           >  1:837433/1‑178 (MQ=255)
                                                                                                                                                                              ttGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGCTCTGGCGCTGGTGCGCCACGCcgacga                                                                                           <  1:1206074/191‑1 (MQ=255)
                                                                                                                                                                                 gCCGATGACTCTGGACTGGCGGGGGACTGCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGc                                                                                                                                                                                              <  2:16449/89‑1 (MQ=255)
                                                                                                                                                                                 gCCGACGACTCTGGACTGGCGGTGGACTGCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGaaa                                                                                                                                                               <  2:163280/120‑1 (MQ=255)
                                                                                                                                                                                        aCTCGGGACTGGCGGTGGACTTCCTCGGGGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCgg                                                                                                                                                                                               <  2:19876/81‑1 (MQ=255)
                                                                                                                                                                                              gACTGGCGGTGGCCTTCCTCGGTGGCGCGCCGGGCATCTACTCGGGGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAg                                                                    <  2:1220709/198‑1 (MQ=255)
                                                                                                                                                                                              gACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAgcggc                              >  1:1013995/1‑236 (MQ=255)
                                                                                                                                                                                                 tGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAg                                                                    <  2:545988/195‑1 (MQ=255)
                                                                                                                                                                                                       tGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGAcgcg                                                                                                                                                                                                <  2:1106839/65‑1 (MQ=255)
                                                                                                                                                                                                       tGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGAcgcg                                                                                                                                                                                                >  1:1106842/1‑65 (MQ=255)
                                                                                                                                                                                                       tGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAgg                       >  2:775627/1‑234 (MQ=255)
                                                                                                                                                                                                           cTTCCTCGGTGGCGCGCCGAGTAGATGCCCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTgc                                                                                                        >  2:494341/1‑149 (MQ=255)
                                                                                                                                                                                                           cTTCCTCGGTGGCGCGCCGAGTAGATGCCCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTgc                                                                                                        <  1:494344/149‑1 (MQ=255)
                                                                                                                                                                                                              ccTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGCTTTGGCt              <  1:392842/236‑1 (MQ=255)
                                                                                                                                                                                                                 cGGTGGCGCGCCGGGCATCTCCTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCtgtg                                                             <  1:356458/186‑1 (MQ=255)
                                                                                                                                                                                                                  ggTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGCTTTGGCTa             >  2:412337/1‑233 (MQ=255)
                                                                                                                                                                                                                   gTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGCTTTGGCCACGAt         <  1:858436/236‑1 (MQ=255)
                                                                                                                                                                                                                      gcgcgcCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGCTTTGGCTACGAt         <  2:643654/233‑1 (MQ=255)
                                                                                                                                                                                                                         cgcCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGCTTTGGCt              <  1:681344/225‑1 (MQ=255)
                                                                                                                                                                                                                         cgcCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGCTTTGGCt              >  2:681341/1‑225 (MQ=255)
                                                                                                                                                                                                                           ccGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGCTTTGGCTACGATCCGCTGt  <  1:775630/235‑1 (MQ=255)
                                                                                                                                                                                                                            cGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTg                                                                                                                                               >  2:1226812/1‑93 (MQ=255)
                                                                                                                                                                                                                                 aTCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGac                                                                     >  2:244352/1‑161 (MQ=255)
                                                                                                                                                                                                                                    tACTCGGCGGGCGATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCCGACCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCtgtg                                                             <  1:1170944/167‑1 (MQ=255)
                                                                                                                                                                                                                                    tACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCCGCAAACAACGCCACGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCtgtg                                                             <  1:248195/167‑1 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                              
GCCAGCCACAATGCCGGCAAACTGAAGGAACTCCAGGCCATGCTCGGCCATTCCGTGCAGCTGCACTCGATTGGTGAGTTCAGCCAGGTGGAGCCGGAAGAAACCGGCCTGTCGTTTGTCGAGAACGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGCTTTGGCTACGATCCGCTGT  >  AE015451/5822499‑5822952

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