Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,400,978 C→T G266S (GGT→AGT)  abgT ← p‑aminobenzoyl glutamate:H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,400,9780CT97.9% 154.0 / NA 48G266S (GGT→AGT) abgTp‑aminobenzoyl glutamate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (24/22);  minor base G (0/1);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CCCGCCATCTCTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACCTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTATCAATTACACTGACGTGCATTTGCGGATTGAACGCAGCTGCCGCTTCCGTGCTGATCCCCGACAGCAACA  >  NC_000913/1400740‑1401208
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                      
cccGCCATCTCTTTCAACGGTTCAATCATGAAATGCGGTAAAGCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATGCCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGGTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTaa                                                                                                                                                                                                                            <  2:276335/251‑1 (MQ=255)
 ccGCCATCTCTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTaa                                                                                                                                                                                                                            >  1:288365/1‑250 (MQ=255)
     cATCTCTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGCGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCa                                                                                                                                                                                                                        >  1:306241/1‑250 (MQ=255)
          ctTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAAc                                                                                                                                                                                                                    <  1:344009/249‑1 (MQ=255)
             taatcgaatCAATCATTAAATGCGGTAAATCAGCCTGAAGTCGAGTTGGGCGGGTACCGATGCCATAAGCAAGCGCGACAACAAAGAACACAAGAATGATCAGTGGCACGATACCTTTAATAAGGGGTGATGGCATCACGGCGTGATTAATCGGCTCGCGAAATATCCCGGTTTGCGGGATCACCATCAGCGCAAACGCAGAAATACAAAGAAGCGATACTACACTTGCTATGCGTAAACCAAAACGCTg                                                                                                                                                                                                                >  2:377013/9‑250 (MQ=255)
                   gTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCg                                                                                                                                                                                                         <  1:314009/251‑1 (MQ=255)
                    ttCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCATACGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGAGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCATTCGCAGCACTAAAAAGTAGCGCTACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCgg                                                                                                                                                                                                        >  2:149670/1‑251 (MQ=255)
                     tCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGt                                                                                                                                                                                                       >  1:224890/1‑251 (MQ=255)
                          cATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATg                                                                                                                                                                                                  >  1:401415/1‑251 (MQ=255)
                          cATTAAATGCGGTAAATCCGACTGACGTCGAATTGTGCGCCTAGCGCTGCCCTAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGCTAACTTTAATCAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGATCactg                                                                                                                                                                                                  >  1:150215/1‑248 (MQ=255)
                             tAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCATACGACAGCGACACAACAAAGAAAAAAAGAATGATCAGTGGCACGATAACTTTAATAAAGGGTGATGGCATCACGGGGAGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCCATCGCAGCAATAAAAAGGAGCGAAACGACAATTGCTATGCGTAAACCAAAACGCTAACTTTCGGTCAATGtat                                                                                                                                                                                               >  2:17994/1‑249 (MQ=255)
                                                    gTCGAATTGTGCGGGTAGCGATGCCATAAGCCAGCGAGACAACAAAGAACAAAAGAATGATCAGTGGCACGCTACCTTTAATAAAGGGTGATGGCCTCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGTTCAATGTCTGCAGTTTCTCATCGCTGTTTccc                                                                                                                                                                        >  1:186664/1‑251 (MQ=255)
                                                       gAATTGTGAGGGTAGCGATGCCATAAGACAGCGAGCCCACACAGAAAAAAAGAATGATCAGTGGCACGATAACTTTAATAAAGGGTGATGGCATCACGGTGAGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCACTGTCTGCAGTTTCTCATCGCTGTTTCCCTGc                                                                                                                                                                     >  2:366777/1‑251 (MQ=255)
                                                             ggCGGGTAGCGATGCCAGAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGAGCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTTTTTCCCTGCCATTGa                                                                                                                                                               <  2:395513/250‑1 (MQ=255)
                                                                    aGCGATGCCAGAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATCCATTTAATAAAGGGTGATGGCATCATGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGAGCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCGGACTTTCGGTCAATTTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGGCCTAAc                                                                                                                                                         <  1:125986/250‑1 (MQ=255)
                                                                       gATGCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGAAGCGATACGCCACTTGCTATGCGTAAACCAAAAAGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACATAACCGTg                                                                                                                                                     >  2:376896/1‑251 (MQ=255)
                                                                          gccatccgCCAGCGAGACCACCACGAAAAAACGAATGACGAGTGGCACGATACCTATACTAACGGGCGCTGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAACCGCAGCAATACAAAGAAGCGATACGACACTTGCTATGCGAAAACCAAAACGCTGACTTTCGGACAATGTCTGCAGTTTCTCATCGCTCTTTCCCTGCCATTGACCTATCCGTgtct                                                                                                                                                  >  2:58235/8‑248 (MQ=255)
                                                                              aaaGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGTCATTGACCTAACCGTGGCTCgat                                                                                                                                              <  2:119999/250‑1 (MQ=255)
                                                                               aaGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCAGCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCgat                                                                                                                                              <  1:253538/250‑1 (MQ=255)
                                                                               aaGCCAGCGAGACAACAAAGAAAAAAAGAAGGATCATTGGCACGATACCCTTAATAAAGGGTGATGGCATCGCGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCGCCAGCAGCGCAATCGCATCAATAAAAAGTATCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGTTGTTTCCCTGCCATTGACCTAAGCGTTGCTCgatg                                                                                                                                             <  2:168763/251‑1 (MQ=255)
                                                                                             acaaAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGGCGGTTAt                                                                                                                               <  2:186502/251‑1 (MQ=255)
                                                                                              caaAGAAAACAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTAt                                                                                                                               >  2:285839/1‑250 (MQ=255)
                                                                                                 aGAAAAAAAGAATGAGCAGTGGCACGATTCCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCTCCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACGTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCTTGCCATTGACCTAACCGTGGCGCGATGATTTTGGCGGTTATCatg                                                                                                                           <  2:306241/251‑3 (MQ=255)
                                                                                                    aaaaaaGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGgccc                                                                                                                        >  2:429529/1‑250 (MQ=255)
                                                                                                       aaaGAATGATCAGTGGCACGCTACCTTTCATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATCTCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTCAACCAAAACGCTGACTTTCGGTCAATGTCTGAAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGTTTAACAGgccgcc                                                                                                                      >  2:268179/1‑250 (MQ=255)
                                                                                                         aGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGCGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCcaac                                                                                                                   >  1:284991/1‑251 (MQ=255)
                                                                                                                  aGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAg                                                                                                          <  1:358376/251‑1 (MQ=255)
                                                                                                                   gTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGt                                                                                                         <  1:262055/251‑1 (MQ=255)
                                                                                                                        aCGATACCTTTAATAAAGGTTGATGGCATCACGGTGTGATTAAGCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACGCTTGCTATGCGTAAACCAAAGCGCTGGCTTGCGGTCGATGTCTGCAGTTTCTCATCGCTGGTTGCCGGCCTTGGACTTAACCTTGGGTCGATGATTTTGTTGTTTATCAGGCTGCCGACAATCGTCATTACGAc                                                                                                    <  2:174526/251‑1 (MQ=255)
                                                                                                                          gATACCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACt                                                                                                   >  2:84797/1‑250 (MQ=255)
                                                                                                                          gATACCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACt                                                                                                   >  1:111966/1‑250 (MQ=255)
                                                                                                                                  taataaAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCt                                                                                           <  2:158070/250‑1 (MQ=255)
                                                                                                                                  taataaAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCt                                                                                           <  1:323677/250‑1 (MQ=255)
                                                                                                                                      aaaGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGcc                                                                                       >  1:90044/1‑250 (MQ=255)
                                                                                                                                      aaaGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACCTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGcc                                                                                       >  2:241922/1‑250 (MQ=255)
                                                                                                                                                ggCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGAGCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCa                                                                            <  1:149670/251‑1 (MQ=255)
                                                                                                                                                             gATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTATCAATTaca                                                               <  1:285839/251‑1 (MQ=255)
                                                                                                                                                                                           gCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTATCAATTACACTGACGTGCATTTGCGGATTGAACGCAGCt                                 <  1:252002/251‑1 (MQ=255)
                                                                                                                                                                                              ggATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGAAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTATCAATTACACTGACGTGCATTTGCGGATTGAACGCAGCTGcc                              >  1:176462/1‑251 (MQ=255)
                                                                                                                                                                                                 tCACCATCAGCGCAATCGCAGCAATAAAACGTAGCGCTACGACACTTGCTATTCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGACCTGGCCATAACATACCCGTTATCAATTACACTGACGTGCATTTGCGGATTGAAAGCAGCTGCCGcc                           >  2:22530/1‑250 (MQ=255)
                                                                                                                                                                                                                 cGCAACAATAAAAAGTAGCGATACGCCATTAGCTATGCGTAATCCAAAACGCTGACTTTCGGTCAAGGTCTACAGTTTCTCATCGCTGTTTCCCTGCCATTGCCCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTATCAATTACACTGACGTGCATTTGCGGATTGAACGCAGCTGCCGCTTCCGTGCTGATCCCCg           <  2:224890/251‑1 (MQ=255)
                                                                                                                                                                                                                   cagcaATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTGTCTCAGCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTGAGCAGGCCGCGAACAATCGTCAGTACTACTACGGAGCTGGCCATAAAATACCAGTTATCAATTACACTGACGTGCATGTGCGGATTGAACGCAGCTGCCGCTTCCGTGCTGATCCCCGAc         <  1:423167/251‑1 (MQ=255)
                                                                                                                                                                                                                        aTAAAAAGTAGCGATACGACAATTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGATGTTTCCATGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTATCAATTACACTGACGTGCATTTGCGGATTGAACGCAGCTGCCGATTCCGTGCTGATCCACGACAGCa     >  2:55093/1‑250 (MQ=255)
                                                                                                                                                                                                                         tAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAACACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTATCAATTACACTGACGTGCATTCGCGGATTGAACGCAGCGGCCGCTTCCGTGATGATCCCCGACAGCAac   >  2:90624/1‑251 (MQ=255)
                                                                                                                                                                                                                         aaaaaatgaGCGATACGACTCTTGCTATGCGTAAACCAAAACGCTGACGTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTGATCTATTACACGGACGTGCATTTGCGGATTGAACGCAGCTGCCGCTTCCGTGCTGAGCCCCGTCAGCAac   <  2:401415/243‑1 (MQ=255)
                                                                                                                                                                                                                          aaaaaGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTATCAATTACACTGACGTGCATTTGCGGATTGAACGCAGCTGCCGCTTCCGTGCTGATCCCCGACAGCAaca  <  1:42349/251‑1 (MQ=255)
                                                                                                                                                                                                                                         gACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTATCAATTACACTGACGTGCATTTGCGGATTGAACGCAGCTGCCGCTTcc                        >  1:56748/1‑214 (MQ=255)
                                                                                                                                                                                                                                         gACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTATCAATTACACTGACGTGCATTTGCGGATTGAACGCAGCTGCCGCTTcc                        <  2:56748/214‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                      
CCCGCCATCTCTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACCTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTATCAATTACACTGACGTGCATTTGCGGATTGAACGCAGCTGCCGCTTCCGTGCTGATCCCCGACAGCAACA  >  NC_000913/1400740‑1401208

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑