Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,252,936 G→A 5.8% R144H (CGC→CAC)  ubiC → chorismate‑‑pyruvate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,252,9360GA 5.8% 217.2 / 2.5 69R144H (CGC→CAC) ubiCchorismate‑‑pyruvate lyase
Reads supporting (aligned to +/- strand):  ref base G (34/31);  new base A (2/2);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

TTCCTGTGTCAACGTTAAGCGGGCCGGAGCTGGCGTTACAAAAATTGGGTAAAACGCCGTTAGGACGCTATCTGTTCACATCATCGACATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGAAAAAAATATGGAGTGGAGTCTGACGCAGAATAAGCTGCTGGCGTTTCATCGCTTAATGCGTACGGATAAGC  >  NC_000913/4252789‑4253079
                                                                                                                                                   |                                                                                                                                               
ttCCTGTGTCAACGTTAAGCGGGCCGGAGCTGGCGTTACAAAAATTGGGTAAAACGCCGTTAGGACGCTATCTGTTCACATCATCGACATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGc                                                                                                                                                >  6:335756/1‑149 (MQ=255)
        tCAACGTTAAGCGGGCCGGAGCTGGCGTTACAAAAATTGGGTAAAACGCCGTTAGGACGCTATCTGTTCACATCATCGACATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGAt                                                                                                                                         >  2:62031/1‑148 (MQ=255)
        tCAACGTTAAGCGGGCCGGAGCTGGCGTTACAAAAATTGGGTAAAACGCCGTTAGGACGCTATCTGTTCACATCATCGACATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCACCTGCGAt                                                                                                                                         >  1:436840/1‑148 (MQ=255)
          aaCGTTAAGCGGGCCTGAGCTGGCGTTACAAAAATTGGGTAAAACGCCGTTAGGACGCTATCTGTTCACATCATCGACATTAACCCGGGACGTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTaa                                                                                                                                      >  4:18151/1‑149 (MQ=255)
           aCGTTAAGCGGGCCGGAGCTGGCGTTACAAAAATTGGGTAAAACGCCGTTAGGACGCTATCTGTTCACATCATCGACATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAg                                                                                                                                     <  8:334572/149‑1 (MQ=255)
             gTTAAGCGGGCCGGAGCTGGCGTTACAAAAATTGGGTAAAACGCCGTTAGGACGCTATCTGTTCACATCATCGACATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCg                                                                                                                                   >  4:410123/1‑149 (MQ=255)
                             cTGGCGTTACAAAAATTGGGTAAAACGCCGTTAGGACGCTATCTGTTCACATCATCGACATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGc                                                                                                                    >  1:53229/1‑148 (MQ=255)
                             cTGGCGTTACAAAAATTGGGTAAAACGCCGTTAGGACGCTATCTGTTCACATCATCGACATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCACCTGCGATTAAGCGGTAAACCGCTGTTGc                                                                                                                    <  2:436840/148‑1 (MQ=255)
                                gCGTTACAAAAATTGGGTAAAACGCCGTTAGGACGCTATCTGTTCACATCATCGACATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTaaa                                                                                                                              >  8:203123/1‑135 (MQ=255)
                                gCGTTACAAAAATTGGGTAAAACGCCGTTAGGACGCTATCTGTTCACATCATCGACATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTaaa                                                                                                                              <  7:203123/135‑1 (MQ=255)
                                gCGTTACAAAAATTGGGTAAAACGCCGTTAGGACGCTATCTGTTCACATCATCGACATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTTTACCGCttaatcgc                                                                                                                  <  3:410123/147‑8 (MQ=255)
                                gCGTTACAAAAATTGGGTAAAACGCCGTTAGGACGCTATCTGTTCACATCATCGACATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTaa                                                                                                                 >  4:344397/1‑148 (MQ=255)
                                 cGTTACAAAAATTGGGTAAAACGCCGTTAGGACGCTATCTGTTCACATCATCGACATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACa                                                                                                               <  2:53229/149‑1 (MQ=255)
                                   ttACAAAAATTGGGTAAAACGCCGTTAGGACGCTATCTGTTCACATCATCGACATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAg                                                                                                              >  7:345698/1‑148 (MQ=255)
                                     aCAAAAATTGGGTAAAACGCCGTTAGGACGCTATCTGTTCACATCATCGACATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGaa                                                                                                            <  7:88146/148‑1 (MQ=255)
                                         aaaTTGGGTAAAACGCCGTTAGGACGCTATCTGTTCACATCATCGACATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTACGCGGTAAACCGCTGTTGCTAACAGAACTGtt                                                                                                       <  1:441840/149‑1 (MQ=255)
                                                   aaaCGCCGTTAGGACGCTATCTGTTCACATCATCGACATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGc                                                                                              <  5:335756/148‑1 (MQ=255)
                                                    aaCGCCGTTAGGACGCTATCTGTTCACATCATCGACATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCg                                                                                             >  1:139323/1‑148 (MQ=255)
                                                    aaCGCCGTTAGGACGCTATCTGTTCACATCATCGACATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGt                                                                                            >  6:278976/1‑149 (MQ=255)
                                                     aCGCCGTTAGGACGCTATCTGTTCACATCATCGACATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGt                                                                                            >  1:278788/1‑148 (MQ=255)
                                                       gcagTTAGGACGCTATCTGTTCACATCATCGACATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGTCGTCAc                                                                                         >  6:291122/4‑149 (MQ=255)
                                                           ttAGGACGCTATCTGTTCACATCATCGACATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGtt                                                                                     <  8:345698/149‑1 (MQ=255)
                                                             aGGACGCTATCTGTTCACATCATCGACATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTg                                                                                    <  3:344397/148‑1 (MQ=255)
                                                                 cGCTATCTGTTCACATCATCGACATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACt                                                                                >  7:249792/1‑148 (MQ=255)
                                                                    tATCTGTTCACATCATCGACATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTaa                                                                              >  1:279525/1‑147 (MQ=255)
                                                                        tGTTCACATCATCGACATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGa                                                                         >  5:55735/1‑148 (MQ=255)
                                                                          ttCACATCATCGACATTAACCCGGGACTTTATTGAGATAGGCCGGGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAgg                                                                          <  5:64081/145‑1 (MQ=255)
                                                                              catcatCGACATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGaaaaaaa                                                                   >  5:46316/1‑148 (MQ=255)
                                                                               atcatcGACATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGaaaaaaa                                                                   >  4:97340/1‑147 (MQ=255)
                                                                               atcatcGACATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGaaaaaaa                                                                   <  3:97340/147‑1 (MQ=255)
                                                                               atcatcGACATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGAAAAAAata                                                                 <  1:62031/149‑1 (MQ=255)
                                                                                tcatcGACATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGAAAAAAatat                                                                <  6:46316/149‑1 (MQ=255)
                                                                                  atcGACATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGTCGTCACCGTTGTACTAAGAGGAAAAAAATAtg                                                               <  2:279525/148‑1 (MQ=255)
                                                                                  atcGACATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGAAAAAAATAtg                                                               >  1:19641/1‑148 (MQ=255)
                                                                                   tcGACATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGAAAAAAATAtg                                                               >  3:289098/1‑147 (MQ=255)
                                                                                     gACATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCACCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGAAAAAAATAtggagt                                                           >  1:365490/1‑149 (MQ=255)
                                                                                      aCATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGAAAAAAATAtggagtg                                                          >  8:71444/1‑149 (MQ=255)
                                                                                      aCATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGAAAAAAATAtggagtg                                                          <  2:98804/149‑1 (MQ=255)
                                                                                        aTTAACCCGGGACTTTATTTAGATAGGCCGTTATGCAGTGCTGTGGGGGCGACGTTCCCGCATGCTATTAAGCGTTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCTTCACCTTTGTACTAAGAGGAAAAAAATAtggagtgga                                                        <  5:291122/149‑1 (MQ=255)
                                                                                        aTTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGAAAAAAATAtggagtgga                                                        <  3:18151/149‑1 (MQ=255)
                                                                                             cccGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACt                                                                                >  1:376168/1‑120 (MQ=255)
                                                                                             cccGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACt                                                                                <  2:376168/120‑1 (MQ=255)
                                                                                              ccGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGAAAAAAATAtggagtggag                                                       <  1:59087/144‑1 (MQ=255)
                                                                                              ccGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGAAAAAAATAtggagtggag                                                       >  2:59087/1‑144 (MQ=255)
                                                                                              ccGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGAAAAAAATATGGAGTGGAGTCTGa                                                  <  7:71444/149‑1 (MQ=255)
                                                                                                  gACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGAAAAAAATATGGAGTGGAGTCTGACg                                                <  8:249792/147‑1 (MQ=255)
                                                                                                     tttATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGAAAAAAATATGGAGTGGAGTCTGACGCa                                              <  7:211368/146‑1 (MQ=255)
                                                                                                     tttATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGAAAAAAATATGGAGTGGAGTCTGACGCa                                              >  8:211368/1‑146 (MQ=255)
                                                                                                            agaTAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGAAAAAAATATGGAGTGGAGTCTGACGCAGAATAAgctg                                    >  8:83400/1‑149 (MQ=255)
                                                                                                             gaTAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGAAAAAAATATGGAGTGGAGTCTGACGCAGAATAAgctgc                                   >  3:175638/1‑149 (MQ=255)
                                                                                                               tAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGtt                                                                                     <  6:197672/97‑1 (MQ=255)
                                                                                                               tAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGtt                                                                                     >  5:197672/1‑97 (MQ=255)
                                                                                                                    cGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGAAAAAAATATGGAGTGGAGTCTGACGCAGAATAAGCTGCTGGCGtt                            >  6:226671/1‑149 (MQ=255)
                                                                                                                    cGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGAAAAAAATATGGAGTGGAGTCTGACGCAGAATAAGCTGCTGGCGtt                            >  3:128381/1‑149 (MQ=255)
                                                                                                                    cGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGGCACCGTTGTACTAAGAGGAAAAAAATATGGAGTGGAGTCTGACGCAGAATAAGCTGCTGGCGtt                            <  2:139323/149‑1 (MQ=255)
                                                                                                                      tGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGAAAAAAATATGGAGTGGAGTCTGACGCa                                              <  2:85650/129‑1 (MQ=255)
                                                                                                                      tGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGAAAAAAATATGGAGTGGAGTCTGACGCa                                              >  1:85650/1‑129 (MQ=255)
                                                                                                                         tGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGAAAAAAATATGGAGTGGAGTCTGACGCAGAATAAGCTGCTGGCGTTTCATc                       >  1:265280/1‑149 (MQ=255)
                                                                                                                             gggCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGAAAAAAATATGGAGTGGAGTCTGACGCAGAATAAGCTGCTGGCGTTTCATCGCtt                   <  5:226671/149‑1 (MQ=255)
                                                                                                                                cTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCTGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGAAAAAAATATGGAGGGGAGTCTGACGCAGAATAAGCTGCTGGCGTTTCATCGCTTAAt                <  6:55735/149‑1 (MQ=255)
                                                                                                                                cTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGAAAAAAATATGGAGTGGAGTCTGACGCAGAATAAGCTGCTGGCGTTTCATCGCTTAAt                <  2:19641/149‑1 (MQ=255)
                                                                                                                                   tgGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGAAAAAAATATGGAGTGGAGTCTGACGCAGAATAAGCTGCTGGCGTTTCATCGCTTAATGCg             >  5:255710/1‑149 (MQ=255)
                                                                                                                                     ggggCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGAAAAAAATATGGAGTGGAGTCTGACGCAGAATAAGCTGCTGGCGTTTCATCGCTTAATGCGt            <  2:278788/148‑1 (MQ=255)
                                                                                                                                     ggggCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGATGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGAAAAAAATATGGAGTGGAGTCTGACGCAGAATAAGCTGCTGGCGTTTCATCGCTTAATGCGTa           >  2:109014/1‑149 (MQ=255)
                                                                                                                                       ggCGACGTTCCCACCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGAAAAAAATATGGAGTGGAGTCTGACGCAGAATAAGCTGCTGGCGTTTCATCGCTTAATGCGTACg         <  2:365490/149‑1 (MQ=255)
                                                                                                                                            cGTTCCCGCCTGCGATTAAGCTGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGAAAAAAATATGGAGTGGAGTCTGACGCAGAATAAGCTGCTGGCGTTTCATCGCTTAATGCGTACGGATaa    <  4:128381/149‑1 (MQ=255)
                                                                                                                                              ttCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGAAAAAAATATGGAGTGGAGTCTGACGCAGAATAAGCTGCTGGCGTTTCATCGCTTAATGCGTACGGATAAGc  >  8:382078/1‑149 (MQ=255)
                                                                                                                                               tCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGAAAAAAATAtggagtg                                                          >  1:241999/1‑92 (MQ=255)
                                                                                                                                               tCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGAAAAAAATAtggagtg                                                          <  2:241999/92‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
TTCCTGTGTCAACGTTAAGCGGGCCGGAGCTGGCGTTACAAAAATTGGGTAAAACGCCGTTAGGACGCTATCTGTTCACATCATCGACATTAACCCGGGACTTTATTGAGATAGGCCGTGATGCCGGGCTGTGGGGGCGACGTTCCCGCCTGCGATTAAGCGGTAAACCGCTGTTGCTAACAGAACTGTTTTTACCGGCGTCACCGTTGTACTAAGAGGAAAAAAATATGGAGTGGAGTCTGACGCAGAATAAGCTGCTGGCGTTTCATCGCTTAATGCGTACGGATAAGC  >  NC_000913/4252789‑4253079

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