Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 324,133 C→T 5.5% Q126* (CAG→TAG)  ykgG → LutC family protein; putative electron transport chain YkgEFG component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913324,1330CT 5.5% 238.9 / 2.8 73Q126* (CAG→TAG) ykgGLutC family protein; putative electron transport chain YkgEFG component
Reads supporting (aligned to +/- strand):  ref base C (34/35);  new base T (2/2);  total (36/37)
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.00e-01

AAGAGCTGGGAGATCAGTCGGTCGTGATTAGCGGTGACACGAGGCTGGAGGAATTGGGGATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTGGTTCTTTTTTCCGCCGCCGAGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTAAAAGCACTATCCTG  >  NC_000913/323993‑324276
                                                                                                                                            |                                                                                                                                               
aaGAGCTGGGAGATCAGTCGGTCGTGATTAGCGGTGACACGAGGCTGGAGGAATTGGGGATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTaaa                                                                                                                                         >  1:201173/1‑149 (MQ=255)
aaGAGCTGGGAGATCAGTCGGTCGTGATTAGCGGTGACACGAGGCTGGAGGAATTGGGGATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTaaa                                                                                                                                         >  3:64608/1‑149 (MQ=255)
 agagCTGGGAGATCAGTCGGTCGTGATTAGCGGTGACACGAGGCTGGAGGAATTGGGGATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAg                                                                                                                                        >  4:433179/1‑149 (MQ=255)
  gagCTGGGAGATCAGTCGGTCGTGATTAGCGGTGACACGAGGCTGGAGGAATTGGGGATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGt                                                                                                                                       <  1:218740/149‑1 (MQ=255)
  gagCTGGGAGATCAGTCGGTCGTGATTAGCGGTGACACGAGGCTGGAGGAATTGGGGATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGt                                                                                                                                       >  2:440081/1‑149 (MQ=255)
    gCTGGGAGATCAGTCGGTCGTGATTAGCGGTGACACGAGGCTGGAGGAATTGGGGATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTgg                                                                                                                                     <  3:90651/149‑1 (MQ=255)
        ggAGATCAGTCGGTCGTGATTAGCGGTGACACGAGGCTGGAGGAATTGGGGATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGgtgt                                                                                                                                 >  3:16925/1‑149 (MQ=255)
         gagaTCAGTCGGTCGTGATTAGCGGTGACACGAGGCTGGAGGAATTGGGGATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGtgtt                                                                                                                                <  4:107888/149‑1 (MQ=255)
           gaTCAGTCGGTCGTGATTAGCGGTGACACGAGGCTGGAGGAATTGGGGATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGtgttgt                                                                                                                              <  2:54262/149‑1 (MQ=255)
           gaTCAGTCGGTCGTGATTAGCGGTGACACGAGGCTGGAGGAATTGGGGATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGtgttgt                                                                                                                              <  3:220936/149‑1 (MQ=255)
           gaTCAGTCGGTCGTGATTAGCGGTGACACGAGGCTGGAGGAATTGGGGATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGtgttgt                                                                                                                              <  4:229514/149‑1 (MQ=255)
             tCAGTCGGTCGTGATTAGCGGTGACACGAGGCTGGAGGAATTGGGGATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTgtgt                                                                                                                            <  2:210431/149‑1 (MQ=255)
             tCAGTCGGTCGTGATTAGCGGTGACACGAGGCTGGAGGAATTGGGGATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTgtgt                                                                                                                            <  1:299397/149‑1 (MQ=255)
             tCAGTCGGTCGTGATTAGCGGTGACACGAGGCTGGAGGAATTGGGGATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTgtgt                                                                                                                            <  4:470257/149‑1 (MQ=255)
                gtcggtcgTGATTAGCGGTGACACGAGGCTGGAGGAATTGGGGATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATg                                                                                                                         >  1:327282/1‑149 (MQ=255)
                     tcgTGATTAGCGGTGACACGAGGCTGGAGGAATTGGGGATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGaa                                                                                                                    >  2:40778/1‑149 (MQ=255)
                     tcgTGATTAGCGGTGACACGAGGCTGGAGGAATTGGGGATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGaa                                                                                                                    <  4:423411/149‑1 (MQ=255)
                        tGATTAGCGGTGACACGAGGCTGGAGGAATTGGGGATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAatat                                                                                                                 >  1:40753/1‑149 (MQ=255)
                           ttAGCGGTGACACGAGGCTGGAGGAATTGGGGATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGt                                                                                                              >  4:196833/1‑149 (MQ=255)
                                    acacGAGGCTGGAGGAATTGGGGATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGaa                                                                                                     >  4:305730/1‑149 (MQ=255)
                                      acGAGGCTGGAGGAATTGGGGATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATc                                                                                                   >  3:188822/1‑149 (MQ=255)
                                           gCTGGAGGAATTGGGGATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAg                                                                                              <  2:201173/149‑1 (MQ=255)
                                             tGGAGGAATTGGGGATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGc                                                                                            <  2:474523/149‑1 (MQ=255)
                                             tGGAGGAATTGGGGATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGc                                                                                            >  3:337242/1‑149 (MQ=255)
                                             tGGAGGAATTGGGGATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGc                                                                                            >  1:61261/1‑149 (MQ=255)
                                             tGGAGGAATTGGGGATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGc                                                                                            >  3:87852/1‑149 (MQ=255)
                                             tGGAGGAATTGGGGATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGc                                                                                            >  1:451987/1‑149 (MQ=255)
                                                aggaATTGGGGATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTg                                                                                         <  1:122242/149‑1 (MQ=255)
                                                 ggaATTGGGGATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTgg                                                                                        <  4:16925/149‑1 (MQ=255)
                                                     ttGGGGATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAg                                                                                              <  3:534089/139‑1 (MQ=255)
                                                     ttGGGGATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAg                                                                                              >  4:534089/1‑139 (MQ=255)
                                                      tGGGGATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTGGTTCtt                                                                                   >  2:474387/1‑149 (MQ=255)
                                                       ggggGTTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGTAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTGGTTCttt                                                                                  >  4:309272/1‑149 (MQ=255)
                                                       ggggATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGTAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTGGTTCttt                                                                                  >  4:347839/1‑149 (MQ=255)
                                                       ggggATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTGGTTCttt                                                                                  <  1:157868/149‑1 (MQ=255)
                                                         ggATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTGGTTCttttt                                                                                <  1:246346/149‑1 (MQ=255)
                                                         ggATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTGGTTCttttt                                                                                >  4:184303/1‑149 (MQ=255)
                                                          gATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTGGTTCtttttt                                                                               <  1:175460/149‑1 (MQ=255)
                                                           aTTAGCCAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGCAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAAGCGTGGTTCTTTTTTc                                                                              <  4:518860/149‑1 (MQ=255)
                                                            ttAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTGGTTCTTTTTTcc                                                                             <  2:40753/149‑1 (MQ=255)
                                                                            gcagGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCa                                                                                                                                                <  2:469531/66‑1 (MQ=255)
                                                                            gcagGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCa                                                                                                                                                >  1:469531/1‑66 (MQ=255)
                                                                              agGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTGGTTCTTTTTTCCGCCGCCGAGCGCGGGCGt                                                           >  1:240592/1‑149 (MQ=255)
                                                                                 aatgcaatgcCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGGGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTGGTTCTTTTTTCCGCCGCCGAGCGCGGGCGTTCa                                                        >  2:116223/1‑149 (MQ=255)
                                                                                   tgcaatgcCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTGGTTCTTTTTTCCGCCGCCGAGCGCGGGCGTTCAtt                                                      <  4:192908/149‑1 (MQ=255)
                                                                                       atgcCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTGGTTCTTTTTTCCGCCGCCGAGCGCGGGCGTTCATTGAGc                                                  <  3:305730/149‑1 (MQ=255)
                                                                                       atgcCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTGGTTCTTTTTTCCGCCGCCGAGCGCGGGCGTTCATTGAGc                                                  <  2:240592/149‑1 (MQ=255)
                                                                                        tgcCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTGGTTCTTTTTTCCGCCGCCGAGCGCGGGCGTTCATTGAGcc                                                 >  3:382246/1‑149 (MQ=255)
                                                                                        tgcCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTGTAACCGAATCGGGAGGCGTGGTTCTTTTTTCCGCCGCCGAGCGCGGGCGTTCATTGAGcc                                                 <  1:471257/149‑1 (MQ=255)
                                                                                              tgttTGGGATCAGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTGGTTCTTTTTTCCGCCGCCGAGCGCGGGCGTTCATTGAGCCTGCTca                                           <  3:196833/149‑2 (MQ=255)
                                                                                               gttTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGTAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTGGTTCTTTTTCCGCCGCCGAGCGCGGGCGTTCATTGAGCCTGCTCCCg                                         <  3:347839/149‑1 (MQ=255)
                                                                                                tttGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGTAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTGGTTCTTTTTTCCGCCGCCGAGCGCGGGCGTTCATTGAGCCTGCTCCCg                                         <  3:309272/149‑1 (MQ=255)
                                                                                                tttGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTGGTTCTTTTTTCCGCCGCCGAGCGCGGGCGTTCATTGAGCCTGCTCCCg                                         <  4:286191/149‑1 (MQ=255)
                                                                                                     gATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGCTGTTGTGTATGCTTAATATGGTTTAACCGAATCGGGAGGCGTGGTTCTTTTTTCCGCCGCCGAGCGCGGGCGTTCATTGAGCCAGCTCCCGGAata                                    <  4:480016/149‑1 (MQ=255)
                                                                                                          ggCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTGGTTCTTTTTTCCGCCGCCGAGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATtctc                               >  1:428165/1‑149 (MQ=255)
                                                                                                             gAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTGGTTCTTTTTTCCGCCGCCGAGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCttt                            <  4:382246/149‑1 (MQ=255)
                                                                                                                aGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTGGTTCTTTTTTCCGCCGCCGAGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTa                         <  1:298768/149‑1 (MQ=255)
                                                                                                                     ccGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTGGTTCTTTTTTCCGCCGCCGAGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTg                    >  1:208117/1‑149 (MQ=255)
                                                                                                                      cGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTGGTTCTTTTTTCCGCCGCCGAGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGc                   >  2:191549/1‑149 (MQ=255)
                                                                                                                      cGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTGGTTCTTTTTTCCGCCGCCGAGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGc                   >  2:205816/1‑149 (MQ=255)
                                                                                                                       gagaATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTGGTTCTTTTTTCCGCCGCCGAGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCg                  >  2:144992/1‑149 (MQ=255)
                                                                                                                       gagaATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTGGTTCTTTTTTCCGCCGCCAAGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCg                  >  3:513944/1‑149 (MQ=255)
                                                                                                                         gaATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTGGTTCTTTTTTCCGCCGCCGAGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTa                >  4:313723/1‑149 (MQ=255)
                                                                                                                          aaTATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTGGTTCTTTTTTCCGCCGCCGAGCGCGGGCGCTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTaa               >  2:84352/1‑149 (MQ=255)
                                                                                                                            tatCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTGGTTCTTTTTTCCGCCGCCGAGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTaaaa             <  4:188822/149‑1 (MQ=255)
                                                                                                                             atCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTGGTTCTTTTTTCCGCCGCCGAGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTAAAAg            <  4:107367/149‑1 (MQ=255)
                                                                                                                                tcGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTGGTTCTTTTTTCCGCCGCCGAGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTAAAAGCAc         >  3:248386/1‑149 (MQ=255)
                                                                                                                                  gcaggcagAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTaa                                                                                                          <  1:70377/50‑1 (MQ=255)
                                                                                                                                  gcaggcagAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTaa                                                                                                          >  2:70377/1‑50 (MQ=255)
                                                                                                                                    aggcagAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTGGTTCTTTTTTCCGCCGCCGAGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTAAAAGCACTATc     <  2:208117/149‑1 (MQ=255)
                                                                                                                                      gcagAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTGGTTCTTTTTTCCGCCGCCGAGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTAAAAGCACTATCCt   >  4:109053/1‑149 (MQ=255)
                                                                                                                                      gcagAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTGGTTCTTTTTTCCGCCGCCGAGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTAAAAGCACTATCCt   >  2:478793/1‑149 (MQ=255)
                                                                                                                                       cagAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTGGTTCTTTTTTCCGCCGCCGAGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTAAAAGCACTATCCTg  <  1:474387/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                               
AAGAGCTGGGAGATCAGTCGGTCGTGATTAGCGGTGACACGAGGCTGGAGGAATTGGGGATTAGCGAACGTTTGCAGCAGGAATGCAATGCCGTTGTTTGGGATCCGGCGAAAGGTGCCGAGAATATCTCGCAGGCAGAGCAGGCTAAAGTGGGTGTTGTGTATGCTGAATATGGTTTAACCGAATCGGGAGGCGTGGTTCTTTTTTCCGCCGCCGAGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTAAAAGCACTATCCTG  >  NC_000913/323993‑324276

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