Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 5,822,725 T→G 5.2% I85S (ATC→AGC) ‡ rdgB → dITP/XTP pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154515,822,7250TG 5.2% 268.1 / 3.7 72I85S (ATC→AGC) ‡rdgBdITP/XTP pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base T (23/45);  new base G (2/2);  total (25/47)
Fisher's exact test for biased strand distribution p-value = 6.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

GTATTGGCCAGCCACAATGCCGGCAAACTGAAGGAACTCCAGGCCATGCTCGGCCATTCCGTGCAGCTGCACTCGATTGGTGAGTTCAGCCAGGTGGAGCCGGAAGAAACCGGCCTGTCGTTTGTCGAGAACGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGCTTTGGCTACGATCCGCTGT  >  AE015451/5822493‑5822952
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                   
gTATTGGCCAGCCACAATGCCGGCAAACTGAAGGAACTCCAGGCCATGCTCGGCCATTCCGTGCAGCTGCACTCGATTGGTGAGTTCAGCCAGGTGGAGCCGGAAGAAACCGGCCTGTCGTTTGTCGAGAACGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTa                                                                                                                                                                                                                                  >  1:1263284/1‑236 (MQ=255)
 tATTGGCCAGCCACAATGCCGGCAAACTGAAGGAACTCCAGGCCATGCTCGGCCATTCCGTGCAGCTGCACTCGATTGGTGAGTTCAGCCAGGTGGAGCCGGAAGAAACCGGCCTGTCGTTTGTCGAGAACGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTAc                                                                                                                                                                                                                                 <  1:784228/236‑1 (MQ=255)
 tATTGGCCAGCCACAATGCCGGCAAACTGAAGGAACTCCAGGCCATGCTCGGCCATTCCGTGCAGCTGCACTCGATTGGTGAGTTCAGCCAGGTGGAGCCGGAAGAAACCGGCCTGTCGTTTGTCGAGAACGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTAc                                                                                                                                                                                                                                 <  1:447605/236‑1 (MQ=255)
      gccagccaCAATGCCGGCAAACTGAAGGAACTCCAGGCCATGCTCGGCCATTCCGTGCAGCTGCACTCGATTGGTGAGTTCAGCCAGGTGGAGCCGGAAGAAACCGGCCTGTCGTTTGTCGAGAACGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCg                                                                                                                                                                                                                              <  1:9509/234‑1 (MQ=255)
                        aaaCTGAAGGAACTCCAGGCCATGCTCGGCCATTCCGTGCAGCTGCACTCGATTGGTGAGTTCAGCCAGGTGGAGCCGGAAGAAACCGGCCTGTCGTTTGTCGAGAACGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCaa                                                                                                                                                                                                          >  1:763484/1‑236 (MQ=255)
                                   aCTCCAGGCCATGCTCGGCCATTCCGTGCAGCTGCACTCGATTGGTGAGTTCAGCCAGGTGGAGCCGGAAGAAACCGGCCTGTCGTTTGTCGAGAACGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTAc                                                                                                                                                                                                                                 <  2:739560/202‑1 (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                                                                                                                                                                                                   <  1:708378/170‑1 (MQ=255)
                                                                                                 aGCCGGAAGAAACCGGCCTGTCGTTTGTCGAGAACGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGAc                                                                                                                                                                                                   >  2:708375/1‑170 (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                                                                                                                                              >  1:1154871/1‑190 (MQ=255)
                                                                                                                                  aCGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGa                                                                                                                                              <  2:1154868/190‑1 (MQ=255)
                                                                                                                                    gCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGt                                                                                                                            >  1:782860/1‑206 (MQ=255)
                                                                                                                                    gCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGt                                                                                                                            <  2:782857/206‑1 (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                                                                                                                                                                                                >  1:1106842/1‑65 (MQ=255)
                                                                                                                                                                                                             tGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGAcgcg                                                                                                                                                                                                <  2:1106839/65‑1 (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              >  2:681341/1‑225 (MQ=255)
                                                                                                                                                                                                                               cgcCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGCTTTGGCt              <  1:681344/225‑1 (MQ=255)
                                                                                                                                                                                                                                 ccGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGCTTTGGCTACGATCCGCTGt  <  1:775630/235‑1 (MQ=255)
                                                                                                                                                                                                                                  cGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTg                                                                                                                                               >  2:1226812/1‑93 (MQ=255)
                                                                                                                                                                                                                                       aTCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGac                                                                     >  2:244352/1‑161 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                   
GTATTGGCCAGCCACAATGCCGGCAAACTGAAGGAACTCCAGGCCATGCTCGGCCATTCCGTGCAGCTGCACTCGATTGGTGAGTTCAGCCAGGTGGAGCCGGAAGAAACCGGCCTGTCGTTTGTCGAGAACGCCATCCTCAAGGCGCGCAACGCCGCCCGCATTTCCGGCCTGCCGGCCTTGGCCGATGACTCTGGACTGGCGGTGGACTTCCTCGGTGGCGCGCCGGGCATCTACTCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGCTTTGGCTACGATCCGCTGT  >  AE015451/5822493‑5822952

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