Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,539,991 A→G 6.3% L275P (CTG→CCG)  cysA ← sulfate/thiosulfate transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,539,9910AG 6.3% 197.1 / 2.9 65L275P (CTG→CCG) cysAsulfate/thiosulfate transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (31/29);  new base G (2/2);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

TCACGACCGTCAGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCAACATGGAACTGCCCGCCGCGAATGGTTCCCTGCAGGCGGTTCACTTCGCCCATAAATTCGAGCACAAAACGGGTCGCCGGTTCGCG  >  NC_000913/2539845‑2540139
                                                                                                                                                  |                                                                                                                                                    
tCACGACCGTCAGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATAGGGTGTAGTGACCTTACGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAg                                                                                                                                                     >  7:195980/1‑148 (MQ=255)
    gACCGTCAGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGCGTGTAGCGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAaga                                                                                                                                                <  5:123846/149‑1 (MQ=255)
      ccGGCAGCGGTTCGTTGAACCACCCCCACGGCTGCACCACCAAATGGGGGTAGTGACCTTTCGGGCTGGCTTCCAGGACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAgaga                                                                                                                                              <  2:274366/149‑1 (MQ=255)
           aGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCATGAAGAGATCCAc                                                                                                                                         <  8:395735/149‑1 (MQ=255)
           aGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCAc                                                                                                                                         <  6:186352/149‑1 (MQ=255)
           aGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGCAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCAc                                                                                                                                         <  1:28666/149‑1 (MQ=255)
            gCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCGGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCAcc                                                                                                                                        >  5:43085/1‑149 (MQ=255)
                  cGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGccc                                                                                                                                   <  3:336056/148‑1 (MQ=255)
                   gTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGCCCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCt                                                                                                                                 >  8:394010/1‑149 (MQ=255)
                           accCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCa                                                                                                                         <  2:21235/149‑1 (MQ=255)
                             cccAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGCCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAgg                                                                                                                       >  5:278611/1‑149 (MQ=255)
                              ccAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATc                                                                                                                                            <  8:168119/127‑1 (MQ=255)
                              ccAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATc                                                                                                                                            >  7:168119/1‑127 (MQ=255)
                              ccAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGgt                                                                                                                      <  4:119141/149‑1 (MQ=255)
                              ccAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGgt                                                                                                                      <  4:164530/149‑1 (MQ=255)
                                   ggCTGCACCACTATTTGGGTTTTGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTa                                                                                                                 <  6:278611/149‑1 (MQ=255)
                                     cTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCGGGCTTCCAGTACCGGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGACGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGc                                                                                                               <  2:206311/149‑1 (MQ=255)
                                         accacTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCGGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAg                                                                                                           >  1:300634/1‑149 (MQ=255)
                                           cacTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCGGGAAgaga                                                                                                                                              >  3:17367/1‑112 (MQ=255)
                                           cacTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCGGGAAgaga                                                                                                                                              <  4:17367/112‑1 (MQ=255)
                                           cacTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCg                                                                                                         >  2:269218/1‑149 (MQ=255)
                                           cacTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAAACGAGGCTGTTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCg                                                                                                         >  6:108093/1‑149 (MQ=255)
                                                cttGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCa                                                                                                     <  8:316630/147‑1 (MQ=255)
                                                      tgtAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATgc                                                                                              <  3:186014/149‑1 (MQ=255)
                                                         aGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATcc                                                                                                                                           >  2:41825/1‑101 (MQ=255)
                                                         aGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATcc                                                                                                                                           <  1:41825/101‑1 (MQ=255)
                                                          gTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTg                                                                                                                                >  8:200195/1‑111 (MQ=255)
                                                          gTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTg                                                                                                                                <  7:200195/111‑1 (MQ=255)
                                                            gACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGcc                                                                                         >  7:41438/1‑148 (MQ=255)
                                                             aCCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCaa                                                                                       >  4:26827/1‑149 (MQ=255)
                                                                attCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGAACCACCGGCCCCTGATACGCAGTTGTGTAGCCCAGCGGCCAGCGATGCGCGCCAAc                                                                                      >  5:341413/2‑147 (MQ=255)
                                                                  tCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCAACATg                                                                                   <  4:87553/148‑1 (MQ=255)
                                                                    gggCTGGCTTCCAGTACCTGTACCGGCATCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCAACATGGaa                                                                                <  8:107306/149‑1 (MQ=255)
                                                                    gggCTGGCTTCCAGTACCTGTACCGGCAGCGGCGACTCGAGGCTGGGACGGCGGCTGATATCCACTTCCCAACGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCAACATGGaa                                                                                <  6:70964/149‑1 (MQ=255)
                                                                     ggctggctTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCAACATGGAAc                                                                               >  1:24707/1‑149 (MQ=255)
                                                                      gctggctTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCAACATGGAACt                                                                              >  8:62458/1‑149 (MQ=255)
                                                                          gctTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCAACATGGAACTGc                                                                            <  7:394010/147‑1 (MQ=255)
                                                                           ctTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCAACATGGAACTGCccg                                                                         >  4:153831/1‑149 (MQ=255)
                                                                             tCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCAACATGGAACTGCccgcc                                                                       >  3:177045/1‑149 (MQ=255)
                                                                               cAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGTAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCAACATGGAACTGCCcgccgc                                                                     >  1:197643/1‑149 (MQ=255)
                                                                                aGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCAACATGGAACTGCCcgccgc                                                                     >  8:200238/1‑148 (MQ=255)
                                                                                 gTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCa                                                                                                            >  5:259797/1‑108 (MQ=255)
                                                                                 gTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCa                                                                                                            <  6:259797/108‑1 (MQ=255)
                                                                                 gTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCAACATGGAACTGCCCGCcgcg                                                                    >  3:57633/1‑148 (MQ=255)
                                                                                    ccTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGa                                                                                                                               <  4:343451/86‑1 (MQ=255)
                                                                                    ccTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGa                                                                                                                               >  3:343451/1‑86 (MQ=255)
                                                                                        tACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCAACATGGAACTGCCCGCCGCGAATGGTTc                                                            <  4:351686/149‑1 (MQ=255)
                                                                                          ccGGCAGCGGCGTATCGAGGCTGGTACGGCGGGTGATATCGACTTCCCAAGGGCGCTGGAAGAGATCCACCGGCCACTGATACGCAGGTGTTTAGCCCAGCGGCCAGCGATGCGCGCCAACATGGAACTGCCCGCCGCGAATGGTTcc                                                           <  2:261690/148‑1 (MQ=255)
                                                                                           cGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCTCAGGAAGAGATACACCGGCCCCTGATACGCAGGGGTGTAGCCCAGCGGCCAGCGATGCGCGCCAACATGGAACTGCCCGCCGCGAATGGTTCCCt                                                         <  3:130118/149‑1 (MQ=255)
                                                                                            ggCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCAACATGGAACTGCCCGCCGCGAATGGTTCCCTg                                                        >  1:285418/1‑149 (MQ=255)
                                                                                                      aaTCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCCGCGGCCAGCGATGCGCGCCAACATGGAACTGCCCGCCGCGAATGGTTACCTGAcggcgg                                                 >  4:222819/1‑146 (MQ=255)
                                                                                                      aaTCGAGGCTGGTACGGCGGCTGATATCCACTTCCAACGGGAGCAGGAAGAGAGCCACCGGCCCCTGAGACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCAACATGGAACTGCCCGCCGCGAATGGTTCCCTGCAGGCgg                                                 <  3:222819/146‑1 (MQ=255)
                                                                                                        tCGAGGCTGGTACGGCGGCTGATATCCACTTCCCTAGGGCTCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCAACATGGAACTGCCCGCCGCGAATGGTTCCCTGCAGGCGGTTCa                                             <  8:41438/148‑1 (MQ=255)
                                                                                                        tCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGAGCGCCAACATGGAACTGCCCGCCGCGAATGGTTCCCTGCAGGCGGTTCAc                                            >  8:206145/1‑149 (MQ=255)
                                                                                                          gAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCGGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCAACATGGAACTGCCCGCCGCGAATGGTTCCCTGCAGGCGGTTCACt                                           <  2:300634/148‑1 (MQ=255)
                                                                                                                         gcTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCAACATGGAACTGCCCGCCGCGAATGGTTCCCTGCAGGCGGTTCACTTCGCCCATAAATTCGa                           <  7:62458/149‑1 (MQ=255)
                                                                                                                         gcTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCAACATGGAACTGCCCGCCGCGAATGGTTCCCTGCAGGCGGTTCACTTCGCCCATAAATTCGa                           <  3:26827/149‑1 (MQ=255)
                                                                                                                                 ccACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCAACATGGAACTGCCCGCCGCGAATGGTTCCCTGCAGGCGGTTCACTTCGCCCATAAATTCGAGCACaaaa                   >  7:247777/1‑149 (MQ=255)
                                                                                                                                  cACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCAACATGGAACTGCCCGCCGCGAATGGTTCCCTGCAGGCGGTTCACTTCGCCCATAAATTCGAGCACAAAAc                  >  7:320479/1‑149 (MQ=255)
                                                                                                                                    cTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCAACATGGAACTGCCCGCCGCGAATGGTTCCCTGCAGGCGGTTCACTTCGCCCATAAATTCGAGCACAAAACgg                <  8:320479/149‑1 (MQ=255)
                                                                                                                                        ccAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCAACATGGAACTGCCCGCCGCGAATGGTTCCCTGCAGGCGGTTCACTTCGCCCATAAATTCGAGCACAAAACGGGTCg            >  4:288425/1‑149 (MQ=255)
                                                                                                                                         cAAGGGCGCAGGAAGAGCTCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCAACATGGCACTGCCCGCCGCGAATGGTTCCCTGCAGGCGGTTCACGTCGCCCATAAATTCGCGCACAAAACGGGTCGc           >  1:100732/1‑149 (MQ=255)
                                                                                                                                          aaGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCAACATGGCACTGCCCGCCGCGAATGGTTCCCTGCAGGCGGTTCACTTCGCCCATAAATTCGAGCACAAAACGGGTCGcc          >  5:307182/1‑149 (MQ=255)
                                                                                                                                          aaGGGCGCAGGAAGAGATCCACCGGCCCCTGAAACGCAGGTGTGTAGCACAGCGGCCAGCGATGCGCGCCAACATGGAACTGCCCGCCGCGAATGGTTCCCTGCAGGCCGTTCCCTGCGCGCATAAATTCGAGCACAAACCGGCACGcc          >  2:353218/1‑149 (MQ=255)
                                                                                                                                                gcAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCAACATGGAACTGCCCGCCGCGAATGGTTCCCTGCAGGCGGTTCACTTCGCCCATAAATTCGAGCACAAAACGGGTCGCCGGTTcg    >  5:24823/1‑149 (MQ=255)
                                                                                                                                                  aGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCAACATGGAACTGCCCGCCGCGAATGGTTCCCTGCAGGCGGTTCACTTCGCCCATATATTCGAGCACAACACGGGTCGCCGGTTcgcg  >  6:99025/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TCACGACCGTCAGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCAACATGGAACTGCCCGCCGCGAATGGTTCCCTGCAGGCGGTTCACTTCGCCCATAAATTCGAGCACAAAACGGGTCGCCGGTTCGCG  >  NC_000913/2539845‑2540139

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