Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 242,344 C→G 9.7% P320P (CCG→CCC)  fadE ← acyl coenzyme A dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913242,3440CG 9.7% 235.7 / 15.6 73P320P (CCG→CCC) fadEacyl coenzyme A dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (31/34);  new base G (3/4);  total (34/39)
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.41e-01

GCACTCCACCAGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGC  >  NC_000913/242197‑242491
                                                                                                                                                   |                                                                                                                                                   
gCACTCCACCAGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCg                                                                                                                                                    >  2:333871/1‑149 (MQ=255)
 cACTCCACCAGCATCCGCCAGCCTTGCCCGGCCGTTTTCGGCCCGCAGATGATGTAATAGATCGGCACGAAGACAACTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCgg                                                                                                                                                   >  4:1055/1‑149 (MQ=255)
    tCCACCAGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCgt                                                                                                                                                >  1:218524/1‑149 (MQ=255)
    tCCACCAGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCgt                                                                                                                                                >  3:77706/1‑149 (MQ=255)
    tCCACCAGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCgt                                                                                                                                                <  2:199949/149‑1 (MQ=255)
     ccaccaGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCgtg                                                                                                                                               <  2:283686/149‑1 (MQ=255)
      caccaGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCgtgg                                                                                                                                              >  3:533362/1‑149 (MQ=255)
         caGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCGGGCgtggaaa                                                                                                                                           <  1:442130/149‑4 (MQ=255)
           gCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCTGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGtt                                                                                                                                         <  4:479025/149‑1 (MQ=255)
           gCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTAACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGtt                                                                                                                                         <  1:56060/149‑1 (MQ=255)
            cATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTg                                                                                                                                        <  1:379800/149‑1 (MQ=255)
            cATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTg                                                                                                                                        <  4:18870/149‑1 (MQ=255)
             aTCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCGGGCgtggaaattgg                                                                                                                                       <  1:439505/149‑8 (MQ=255)
             aTCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTgg                                                                                                                                       <  1:571295/149‑1 (MQ=255)
                cGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGAt                                                                                                                                    <  2:480962/149‑1 (MQ=255)
                 gCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATc                                                                                                                                   >  4:174492/1‑149 (MQ=255)
                  ccAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCa                                                                                                                                  <  2:104679/149‑1 (MQ=255)
                          gCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAg                                                                                                                          <  3:500928/149‑1 (MQ=255)
                              ggCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTaa                                                                                                                      >  1:533068/1‑149 (MQ=255)
                              ggCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTaa                                                                                                                      <  1:435246/149‑1 (MQ=255)
                               gCCATTCTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAAt                                                                                                                     >  4:484418/1‑149 (MQ=255)
                                ccATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATg                                                                                                                    >  4:619968/1‑149 (MQ=255)
                                  aTTTTCGGCCCTCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGcc                                                                                                                  >  1:74034/1‑149 (MQ=255)
                                   ttttCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCt                                                                                                                 >  2:84742/1‑149 (MQ=255)
                                        ggCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAAtc                                                                                                            <  1:586145/149‑1 (MQ=255)
                                           cctccgATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAAtcttc                                                                                                         <  2:74034/146‑1 (MQ=255)
                                           ccgccgATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAAtcttc                                                                                                         <  3:139811/149‑1 (MQ=255)
                                            cgccgATGATGTAATCGATCGGCGCGAAGACAGCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAAtcttct                                                                                                        >  3:391549/1‑149 (MQ=255)
                                            cgccgATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAAtcttct                                                                                                        <  1:324409/149‑1 (MQ=255)
                                                        aaTCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGt                                                                                            >  4:642114/1‑149 (MQ=255)
                                                         atcgatcgGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTa                                                                                           >  1:579741/1‑149 (MQ=255)
                                                           cgatcgGCACGAAGACCTCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCATTTTCCACCCCAGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAAt                                                                                         <  1:447343/149‑1 (MQ=255)
                                                              tcgGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAAtttt                                                                                      >  1:508422/1‑149 (MQ=255)
                                                                 gCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTcc                                                                                   <  3:642114/149‑1 (MQ=255)
                                                                 gCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTcc                                                                                   <  4:360461/149‑1 (MQ=255)
                                                                   aCGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCgg                                                                                 <  3:326735/149‑1 (MQ=255)
                                                                       aGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCg                                                                             <  4:172714/149‑1 (MQ=255)
                                                                         aCATCTTTACCGCGCGTCAGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGga                                                                           <  4:480602/149‑1 (MQ=255)
                                                                           aTCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGgaga                                                                         <  3:174492/149‑1 (MQ=255)
                                                                            tCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCAccgc                                                                                                 >  2:47972/1‑124 (MQ=255)
                                                                            tCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCAccgc                                                                                                 <  1:47972/124‑1 (MQ=255)
                                                                                  ccGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTaaa                                                                  <  3:484418/149‑1 (MQ=255)
                                                                                                gTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCa                                                    >  2:279664/1‑149 (MQ=255)
                                                                                                 ttCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCAc                                                   >  2:260462/1‑149 (MQ=255)
                                                                                                        aCGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCa                                            >  3:213913/1‑149 (MQ=255)
                                                                                                          ggTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAAt                                          >  4:152290/1‑149 (MQ=255)
                                                                                                           gTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATc                                         <  1:419345/149‑1 (MQ=255)
                                                                                                           gTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGCTTGGGATCAGATCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTGTGCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATc                                         <  4:420091/149‑1 (MQ=255)
                                                                                                            tACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCg                                        >  1:595386/1‑149 (MQ=255)
                                                                                                               gTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCa                                                        <  3:198418/130‑1 (MQ=255)
                                                                                                               gTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCa                                                        >  4:198418/1‑130 (MQ=255)
                                                                                                               gTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTg                                     >  4:78493/1‑149 (MQ=255)
                                                                                                                  cAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCa                                  >  2:601728/1‑149 (MQ=255)
                                                                                                                    gCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGc                                <  4:533362/149‑1 (MQ=255)
                                                                                                                              gCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGc                      <  1:404461/149‑1 (MQ=255)
                                                                                                                               cgacgaCCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCg                     <  2:511354/149‑1 (MQ=255)
                                                                                                                                gacgacCGATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGt                    >  2:377970/1‑149 (MQ=255)
                                                                                                                                    acCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAGTGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGt                >  1:650947/1‑149 (MQ=255)
                                                                                                                                    acCAATTCCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTGCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGt                <  2:422037/149‑1 (MQ=255)
                                                                                                                                      cAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTc              >  3:376246/1‑149 (MQ=255)
                                                                                                                                       aaccaccACGCCGGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTcc             >  3:510099/6‑149 (MQ=255)
                                                                                                                                       aaTTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTcc             >  4:218002/1‑149 (MQ=255)
                                                                                                                                       aaTTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTcc             <  2:365140/149‑1 (MQ=255)
                                                                                                                                        aTTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCa            <  4:627435/149‑1 (MQ=255)
                                                                                                                                         tttCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTaaa                                                                  >  4:313777/1‑94 (MQ=255)
                                                                                                                                         tttCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTaaa                                                                  <  3:313777/94‑1 (MQ=255)
                                                                                                                                         tttCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAg           >  3:264708/1‑149 (MQ=255)
                                                                                                                                         tttCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAg           <  1:260462/149‑1 (MQ=255)
                                                                                                                                         ccaccACGCCGGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTg                 >  3:120818/4‑143 (MQ=255)
                                                                                                                                         ccaccACGCCGGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTg                 <  4:120818/140‑1 (MQ=255)
                                                                                                                                             cACGCCGGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGgaga                                                                         >  4:491228/1‑83 (MQ=255)
                                                                                                                                             cACGCCGGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGgaga                                                                         <  3:491228/83‑1 (MQ=255)
                                                                                                                                             cACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCa       >  3:638486/1‑149 (MQ=255)
                                                                                                                                                  ccGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAAATTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCTAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGc  >  4:649991/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCACTCCACCAGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGC  >  NC_000913/242197‑242491

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