Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,895,825 G→A 6.2% V341I (GTC→ATC)  pabB → aminodeoxychorismate synthase, subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,895,8250GA 6.2% 201.5 / 2.5 64V341I (GTC→ATC) pabBaminodeoxychorismate synthase, subunit I
Reads supporting (aligned to +/- strand):  ref base G (31/29);  new base A (2/2);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.35e-01

GAACTCAGCGAAAGATCGTGCCGAAAATCTGATGATTGTCGATTTAATGCGTAATGATATCGGTCGTGTTGCCGTAGCAGGTTCGGTAAAAGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCTGGTGGCTCAATAACCGGGGCTCCGAAAGTACGGGCTATGGAAATTATCGACGAACTGGAACCGCAGCGACGCAATGCCT  >  NC_000913/1895686‑1895972
                                                                                                                                           |                                                                                                                                                   
gAACTCAGCGAAAGATCGTGCCGAAAATCTGATGATTGTCGATTTAATGCGTAATGATATCGGTCGTGTTGCCGTAGCAGGTTCGGTAAAAGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCAcc                                                                                                                                             >  2:285387/1‑148 (MQ=255)
   cTCAGCGAAAGATCGTGCCGAAAATCTGATGATTGTCGATTTAATGCGTAATGATATCGGTCGTGTTGCCGTAGCAGGTTCGGTAAAAGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATaa                                                                                                                                         >  3:265495/1‑149 (MQ=255)
           aaGATCGTGCCGAAAATCTGATGATTGTCGATTTAATGCGTAATGATATCGGTCGTGTTGCCGTAGCAGGTTCGGTAAAAGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCaa                                                                                                                                 >  3:339197/1‑149 (MQ=255)
                  tGCCGAAAATCTGATGATTGTCGATTTAATGCGTAATGATATCGGTCGTGTTGCCGTAGCAGGTTCGGTAAAAGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGt                                                                                                                                                    <  6:173254/123‑1 (MQ=255)
                  tGCCGAAAATCTGATGATTGTCGATTTAATGCGTAATGATATCGGTCGTGTTGCCGTAGCAGGTTCGGTAAAAGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGt                                                                                                                                                    >  5:173254/1‑123 (MQ=255)
                    ccGAAAATCTGATGATTGTCGATTTAATGCGTAATGATATCGGTCGTGTTGCCGTAGCAGGTTCGGTAAAAGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGaa                                                                                                                        <  7:221019/149‑1 (MQ=255)
                    ccGAAAATCTGATGATTGTCGATTTAATGCGTAATGATATCGGTCGTGTTGCCGTAGCAGGTTCGGTAAAAGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGaa                                                                                                                        <  7:227241/149‑1 (MQ=255)
                      gAAAATCTGATGATTGTCGATTTAATGCGTAATGATATCGGTCGTGTTGCCGTAGCAGGTTCGGTAAAAGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGaa                                                                                                                        >  3:312922/1‑147 (MQ=255)
                                 gatTGTCGATTTAATGCGTAATGATATCGGTCGTGTTGCCGTAGCAGGTTCGGTAAAAGTAAAAGAGCTGTTCGTGGTGGAACACTTCCCTGCCGTGCATCATCATGTCAGAACCATAAAGGCGCAACTACCAGAACAGTTAAACGc                                                                                                             >  1:16893/1‑147 (MQ=255)
                                   ttGTCGATTTAATGCGTAATGATATCGGTCGTGTTGCCGTAGCAGGTTCGGTAAAAGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGc                                                                                                         >  5:180744/1‑149 (MQ=255)
                                             aaTGCGTAATGATATCGGTCGTGTTGCCGTAGCAGGTTCGGTAAAAGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCa                                                                                                                           >  8:225152/1‑121 (MQ=255)
                                             aaTGCGTAATGATATCGGTCGTGTTGCCGTAGCAGGTTCGGTAAAAGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCa                                                                                                                           <  7:225152/121‑1 (MQ=255)
                                             aaTGCGTAATGATATCGGTCGTGTTGCCGTAGCAGGTTCGGTAAAAGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCtgctgc                                                                                               >  8:211633/1‑149 (MQ=255)
                                              aTGCGTAATGATATCGGTCGTGTTGCCGTAGCAGGTTCGGTAAAAGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCtgctgc                                                                                               <  5:9200/148‑1 (MQ=255)
                                               tGCGTAATGATATCGGTCGGGTTGCCGTAGCAGGTTCGGTAAAAGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGGCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTgcgc                                                                                             <  5:201975/149‑1 (MQ=255)
                                                  gTAATGATATCGGTCGTGTTGCCGTAGCAGGTTCGGTAAAAGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCa                                                                                          <  7:50382/149‑1 (MQ=255)
                                                    aaTGATATCGGTCGTGTTGCCGTAGCAGGTTCGGTAAAAGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCa                                                                                          >  2:324808/1‑147 (MQ=255)
                                                    aaTGATATCGGTCGTGTTGCCGTAGCAGGTTCGGTAAAAGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCa                                                                                          <  1:324808/147‑1 (MQ=255)
                                                      tGATATCGGTCCTGATGCCGTAGCATGTTCGGAAAAAGTACCAGAGCTGTTCGTGGTGGAACACTTCCCTGCCGTGCATAATCTGGTCAGCCCCATAACGGCGCAACTACCAGAACAGTTACCCGCCAGCGATCTGCTGCGCGCCGCtt                                                                                      >  2:194708/1‑149 (MQ=255)
                                                           tCGGTCGTGTTGCCGTAGCAGGTTCGGTAAAAGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTcc                                                                                  <  7:208768/148‑1 (MQ=255)
                                                             ggTCGTGTTGTCGTAGCAGGTTCGGTAAAAGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCtgg                                                                               <  5:266702/149‑1 (MQ=255)
                                                             ggTCGTGTTGCCGTAGCAGGTTCGGTAAAAGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCtgg                                                                               <  3:279865/149‑1 (MQ=255)
                                                              gTCGTGTTGCCGTAGCAGGTTCGGTAAAAGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGt                                                                                                                                                    <  6:334742/79‑1 (MQ=255)
                                                              gTCGTGTTGCCGTAGCAGGTTCGGTAAAAGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGt                                                                                                                                                    >  5:334742/1‑79 (MQ=255)
                                                                   gtTGCCGTAGCAGGTTCGGTAAAAGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCTGGTGGCt                                                                          >  1:58842/1‑148 (MQ=255)
                                                                       ccGTAGCAGGTTCGGTAAAAGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCTGGTGGCTCAAt                                                                      <  8:393681/148‑1 (MQ=255)
                                                                         gTAGCAGGTTCGGTAAAAGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCTGGTGGCTCAATAAc                                                                   <  1:285387/149‑1 (MQ=255)
                                                                          tAGCAGGTTCGGTAAAAGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCTGGTGGCTCAATAAcc                                                                  >  6:164701/1‑149 (MQ=255)
                                                                            gCAGGTTCGGTAAAAGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCTGGTGGCTCAATAACCgg                                                                <  1:384848/149‑1 (MQ=255)
                                                                              aGGTTCGGTAAAAGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCTGGTGGCTCAATAACCgggg                                                              <  4:339197/149‑1 (MQ=255)
                                                                              aGGTTCGGTAAAAGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCTGGTGGCTCAATAACCgggg                                                              <  1:311469/149‑1 (MQ=255)
                                                                                   cGGTAAAAGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCTGGTGGCTCAATAACCGGGGCTcc                                                          >  7:100033/1‑148 (MQ=255)
                                                                                     gTAAAAGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCTGGTGGCTCAATAACCGGGGCTCCg                                                         >  5:323241/1‑147 (MQ=255)
                                                                                     gTAAAAGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCTGGTGGCTCAATAACCGGGGCTCCg                                                         <  6:323241/147‑1 (MQ=255)
                                                                                     gTAAAAGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCTGGTGGCTCAATAACCGGGGCTCCGa                                                        >  7:254626/1‑148 (MQ=255)
                                                                                     gTAAAAGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCTGGTGGCTCAATAACCGGGGCTCCGa                                                        >  7:407651/1‑148 (MQ=255)
                                                                                      tAAAAGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCCCCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCGTTTCCTGGTGGCTCAATAACCGGGGCTCCGa                                                        >  8:180633/1‑147 (MQ=255)
                                                                                        aaaGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCTGGTGGCTCAATAACCgggg                                                              <  5:11413/139‑1 (MQ=255)
                                                                                        aaaGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCTGGTGGCTCAATAACCgggg                                                              >  6:11413/1‑139 (MQ=255)
                                                                                        aaaGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCTGGTGGCTCAATAACCGGGGCTCCGAAAg                                                     >  6:344693/1‑148 (MQ=255)
                                                                                        aaaGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCAGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGGTTTTCCTGGTGGCTCAATAACCGGGGCTCCGAAAg                                                     >  8:97661/1‑148 (MQ=255)
                                                                                         aaGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCTGGTGGCTCAATAACCGGGGCTCCGAAAGTa                                                   >  8:24432/1‑149 (MQ=255)
                                                                                           gTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTAGAAGAACAGTTACACGCCTGCGATCTTCTGCGCGCAGCTTTTCCTGGTGGCTCAATAACCGGGGCTCCGAAAGTACg                                                 <  2:16893/149‑1 (MQ=255)
                                                                                           gTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCTGGTGGCTCAATAACCGGGGCTCCGAAAGTACg                                                 <  1:137069/149‑1 (MQ=255)
                                                                                            tACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCTGGTGGCTCAATAACCGGGGCTCCGAAAGTACgg                                                <  4:312922/149‑1 (MQ=255)
                                                                                                agagCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGATCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTc                                                                                   >  4:165830/1‑110 (MQ=255)
                                                                                                agagCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGATCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTc                                                                                   >  4:160241/1‑110 (MQ=255)
                                                                                                agagCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGATCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTc                                                                                   <  3:165830/110‑1 (MQ=255)
                                                                                                agagCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGATCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTc                                                                                   <  3:160241/110‑1 (MQ=255)
                                                                                                       ttCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCTGGTGGCTCAATAACCGGGGCTCCGAAAGTACGGGCTATGGaa                                       <  7:211633/147‑1 (MQ=255)
                                                                                                           tggtggCACCCTTCCCTGCCGTGCATCATCTGGTCAGCCCCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCTGGTGGCTCAATAACCGGGGCTCCGAATGTACGGGCTATGGAAATTATc                                 >  1:326517/1‑149 (MQ=255)
                                                                                                           tggtggAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCTGGTGGCTCAATAACCGGGGCTCCGAAAGTACGGGCTATGGAAATTATc                                 <  7:76706/149‑1 (MQ=255)
                                                                                                               ggAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAg                                                                                                                     >  8:41287/1‑61 (MQ=255)
                                                                                                                gAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCTGGTGGCTCAATAACCGGGGCTCCGAAAGTACGGGCTATGGAAATTATcgacg                             >  8:347853/1‑148 (MQ=255)
                                                                                                                  aCCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCTGGTGGCTCAATAACCGGGGCTCCGAAAGTACGGGCTATGGAAATTATCGACGAAc                          >  8:29110/1‑149 (MQ=255)
                                                                                                                   cccTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCTGGTGGCTCAATAACCGGGGCTCCGAAAGTACGGGCTATGGAAATTATCGACGAACt                         >  7:193572/1‑149 (MQ=255)
                                                                                                                        cccTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCTGGTGGCTCAATACCCGGGGCTCCGAAAGTACGGGCTATGGAAATTATCGACGAACTGGaa                     <  4:265495/148‑1 (MQ=255)
                                                                                                                        cccTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCTGGTGGCTCAATAACCGGGGCTCCGAAAGTACGGGCTATGGAAATTATCGACGAACTGGAAc                    >  6:330545/1‑149 (MQ=255)
                                                                                                                         ccTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCTGGTGGCTCAATAACCGGGGCTCCGAAAGTACGGGCTATGGAAATTATCGACGAACTGGAAcc                   >  1:153661/1‑149 (MQ=255)
                                                                                                                          cTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCTGGTGGCTCAATAACCGGGGCTCCGAAAGTACGGGCTATGGAAATTATCGACGAACTGGAACCg                  <  5:363526/149‑1 (MQ=255)
                                                                                                                              cGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCTGGTGGCTCAATAACCGTGGCTCCGAAAGTACGGGCTATGGAAATTATCGACGAACTGGAACCGCAGc              <  5:164701/149‑1 (MQ=255)
                                                                                                                                      atcTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCTGGTGGCTCAATAACCGGGGCTCCGAAAGTACGGGCTATGGAAATTATCGACGAACTGGAACCGCAGCGACGCaa       <  8:100033/148‑1 (MQ=255)
                                                                                                                                      atcTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCTGGTGGCTCAATAACCGGGGCTCCGAAAGTACGGGCTATGGAAATTATCGACGAACTGGAACCGCAGCGACGCaa       <  7:24432/148‑1 (MQ=255)
                                                                                                                                         tGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCTGGTGGCTCAATTACCGGGGCTCCGAATGTACGGGCTATGGAAATTATCGACGAACTGGAACCGCAGCGACGCAATGcg   >  8:327495/1‑148 (MQ=255)
                                                                                                                                           gTCAGCACCATATCGGCGCAACTTCGAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCTGGTGGCTCAAAAACCGGGGCTCCGAAAGTACGGGCTATGGAAATTATCGACGAACTGGAACCGCAGCGACGCAATGCCt  <  1:194708/148‑1 (MQ=255)
                                                                                                                                           gTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCTGGTGGCTCAATAACCGGGGCTCCGAAAGTACGGGCTATGGAAATTATCGACGAACTGGAACCGCAGCGACGCAATGc    <  8:254626/146‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                   
GAACTCAGCGAAAGATCGTGCCGAAAATCTGATGATTGTCGATTTAATGCGTAATGATATCGGTCGTGTTGCCGTAGCAGGTTCGGTAAAAGTACCAGAGCTGTTCGTGGTGGAACCCTTCCCTGCCGTGCATCATCTGGTCAGCACCATAACGGCGCAACTACCAGAACAGTTACACGCCAGCGATCTGCTGCGCGCAGCTTTTCCTGGTGGCTCAATAACCGGGGCTCCGAAAGTACGGGCTATGGAAATTATCGACGAACTGGAACCGCAGCGACGCAATGCCT  >  NC_000913/1895686‑1895972

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