Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,492,273 T→C 5.9% intergenic (+27/+68) cysNC → / ← imm sulfate adenylyltransferase subunit 1/adenylyl‑sulfate kinase/colicin‑E7 immunity protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,492,2730TC 5.9% 643.5 / 18.1 204intergenic (+27/+68)cysNC/immsulfate adenylyltransferase subunit 1/adenylyl‑sulfate kinase/colicin‑E7 immunity protein
Reads supporting (aligned to +/- strand):  ref base T (99/93);  new base C (5/7);  total (104/100)
Fisher's exact test for biased strand distribution p-value = 5.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.19e-01

Only 100 of 209 total aligned reads displayed.

GTGAAAGCTTCCCGTTCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACTTTTTCAGGTTTACCATTT  >  NC_002947/1492125‑1492422
                                                                                                                                                    |                                                                                                                                                      
   aaaGCTTCCCGTTCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTcgg                                                                                                                                                    >  2:826607/1‑150 (MQ=255)
    aaGCTTCCCGTTCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCGGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTcggc                                                                                                                                                   >  1:2599027/1‑150 (MQ=255)
       cTTCCCGTTCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCggg                                                                                                                                                <  2:2628447/150‑1 (MQ=255)
       cTTCCCGTTCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCggg                                                                                                                                                >  2:3755999/1‑150 (MQ=255)
       cTTCCCGTTCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCggg                                                                                                                                                <  1:2404811/150‑1 (MQ=255)
          cccGTTCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCCGAAGCGGGGtt                                                                                                                                             <  2:2983272/150‑1 (MQ=255)
          cccGTTCGATGTGCCGCTGGATGCTGATATTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCtt                                                                                                                                                       >  2:5093247/1‑140 (MQ=255)
            cGTTCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGtttt                                                                                                                                           >  1:3596873/1‑150 (MQ=255)
            cGTTCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGtttt                                                                                                                                           >  2:3185369/1‑150 (MQ=255)
            cGTTCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCCGAAGCGGGGtttt                                                                                                                                           <  1:3923229/150‑1 (MQ=255)
             gTTCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGttttt                                                                                                                                          >  2:1406868/1‑150 (MQ=255)
              ttCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGAGGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTGCGGCGGGGttgtgt                                                                                                                                         >  2:4425556/1‑146 (MQ=255)
              ttCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGtttttt                                                                                                                                         <  2:5335213/150‑1 (MQ=255)
              ttCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGtttttt                                                                                                                                         <  2:2486133/150‑1 (MQ=255)
               tCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTc                                                                                                                                        <  2:4290947/150‑1 (MQ=255)
                 gATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTa                                                                                                                                      >  2:676947/1‑150 (MQ=255)
                  aTGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCCAGCCTGCAAGAAAAACCCCGCCGAAGCGGGGttttccttgc                                                                                                                                     <  2:3911499/150‑7 (MQ=255)
                        cGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCa                                                                                                                               <  1:3196889/150‑1 (MQ=255)
                             gATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCg                                                                                                                         >  1:4928548/1‑150 (MQ=255)
                              aTGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCgg                                                                                                                        <  1:634390/150‑1 (MQ=255)
                                gCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGt                                                                                                                      >  2:3486903/1‑150 (MQ=255)
                                  tGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGa                                                                                                                    >  2:2451430/1‑150 (MQ=255)
                                   gATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGaa                                                                                                                   <  2:64194/150‑1 (MQ=255)
                                     tCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGaaaa                                                                                                                 <  1:4521006/150‑1 (MQ=255)
                                       ttGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAgg                                                                                                               <  2:770946/150‑1 (MQ=255)
                                       ttGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGAAAAACCCCGCCGAAGCGGGGttt                                                                                                                                            <  1:379143/122‑1 (MQ=255)
                                       ttGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGAAAAACCCCGCCGAAGCGGGGttt                                                                                                                                            >  2:379143/1‑122 (MQ=255)
                                        tgtgATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGc                                                                                                              >  1:1899836/1‑150 (MQ=255)
                                           gATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCt                                                                                                           >  1:179417/1‑150 (MQ=255)
                                           gATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCt                                                                                                           >  2:2556560/1‑150 (MQ=255)
                                            aTCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTc                                                                                                          <  1:3831561/150‑1 (MQ=255)
                                            aTCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTc                                                                                                          >  2:476906/1‑150 (MQ=255)
                                             tCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCa                                                                                                         <  2:3154614/150‑1 (MQ=255)
                                             tCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCa                                                                                                         <  1:2789969/150‑1 (MQ=255)
                                              cGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCag                                                                                                        >  1:2486640/1‑150 (MQ=255)
                                                  acCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGtt                                                                                                    >  2:1272757/1‑150 (MQ=255)
                                                  acCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATCGCTCAGTTCGGGTGAAAAGGCCCTCAGAGtt                                                                                                    >  2:4884149/1‑150 (MQ=255)
                                                    ccAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCt                                                                                                                                       <  1:5259305/114‑1 (MQ=255)
                                                    ccAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCa                                                                                                  >  1:316932/1‑150 (MQ=255)
                                                     cAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAg                                                                                                 >  1:1253898/1‑150 (MQ=255)
                                                        gCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGc                                                                                              >  1:2018364/1‑150 (MQ=255)
                                                         ccaccaGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGcc                                                                                             <  2:1051244/150‑1 (MQ=255)
                                                           accaGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGGTGAAAAGGCCCTCAGAGTTCATGGCCt                                                                                            <  2:4129030/150‑1 (MQ=255)
                                                             caGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGc                                                                                         >  1:2345737/1‑150 (MQ=255)
                                                                 gTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCgt                                                                                     <  1:4470673/150‑1 (MQ=255)
                                                                   ggACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGTCCTTGCTCgtgt                                                                                   <  2:2746052/150‑1 (MQ=255)
                                                                   ggACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCgtgt                                                                                   <  1:1398597/150‑1 (MQ=255)
                                                                    gACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGtt                                                                                  >  1:925324/1‑150 (MQ=255)
                                                                      cGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAg                                                                                >  1:3578055/1‑150 (MQ=255)
                                                                       gAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGc                                                                               <  2:4182815/150‑1 (MQ=255)
                                                                        aaGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCt                                                                              <  1:4629453/150‑1 (MQ=255)
                                                                        aaGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCt                                                                              >  1:2674648/1‑150 (MQ=255)
                                                                         aGGCGTGAAACAGGTGCTGGACGTGATGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTg                                                                             <  1:890965/150‑1 (MQ=255)
                                                                          ggCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTgg                                                                            <  1:5016133/150‑1 (MQ=255)
                                                                            cGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGc                                                                          <  1:1617203/150‑1 (MQ=255)
                                                                             gTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCt                                                                         >  2:4069858/1‑150 (MQ=255)
                                                                              tGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCtt                                                                        <  1:1223041/150‑1 (MQ=255)
                                                                                aaaCAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAACAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGa                                                                      <  2:4218811/150‑1 (MQ=255)
                                                                                aaaCAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGa                                                                      >  2:4801375/1‑150 (MQ=255)
                                                                                  aCAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGaaa                                                                    >  1:2598212/1‑150 (MQ=255)
                                                                                  aCAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGaaa                                                                    <  1:4096278/150‑1 (MQ=255)
                                                                                    aGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAAcc                                                                  <  1:514165/150‑1 (MQ=255)
                                                                                      gTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGt                                                                >  2:2373378/1‑150 (MQ=255)
                                                                                       tGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTc                                                               <  2:524637/150‑1 (MQ=255)
                                                                                         cTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCtt                                                             <  1:4678756/150‑1 (MQ=255)
                                                                                          tGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCttt                                                            >  1:1542026/1‑150 (MQ=255)
                                                                                           ggACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTa                                                           <  1:2640952/150‑1 (MQ=255)
                                                                                             aCGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTAcc                                                         <  1:2985788/150‑1 (MQ=255)
                                                                                               gTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGt                                                       <  1:724406/150‑1 (MQ=255)
                                                                                                         aGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAgccg                                             >  1:1452098/1‑150 (MQ=255)
                                                                                                            gTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCa                                          >  1:4755799/1‑150 (MQ=255)
                                                                                                                cgcgATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGc                                      <  2:5299048/150‑1 (MQ=255)
                                                                                                                  cgATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGTGTTTTTTCTTTTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCtt                                    <  2:3235061/150‑1 (MQ=255)
                                                                                                                  cgATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCtt                                    <  2:4593342/150‑1 (MQ=255)
                                                                                                                  cgATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCtt                                    <  2:1198975/150‑1 (MQ=255)
                                                                                                                   gATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCtt                                    <  1:4898622/149‑1 (MQ=255)
                                                                                                                    aTCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGa                                  <  1:5057439/150‑1 (MQ=255)
                                                                                                                     tCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGAc                                 <  2:2010962/150‑1 (MQ=255)
                                                                                                                      cTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTc                                                                                                                                        <  1:2070501/47‑1 (MQ=255)
                                                                                                                      cTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTc                                                                                                                                        >  2:2070501/1‑47 (MQ=255)
                                                                                                                       tAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTa                               >  1:4039490/1‑150 (MQ=255)
                                                                                                                       tAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTa                               >  2:2191613/1‑150 (MQ=255)
                                                                                                                        aaGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTAt                              >  2:5156412/1‑150 (MQ=255)
                                                                                                                          gcgcTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTAttt                            <  1:104901/150‑1 (MQ=255)
                                                                                                                           cgcTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTc                           >  2:5414459/1‑150 (MQ=255)
                                                                                                                             cgaGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTa                         <  2:3580507/148‑1 (MQ=255)
                                                                                                                             cTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTa                         <  1:3521195/150‑1 (MQ=255)
                                                                                                                             cTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTa                         <  1:4425556/150‑1 (MQ=255)
                                                                                                                              tAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGc                                                                                              <  2:1838568/80‑1 (MQ=255)
                                                                                                                                    gCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACttt                  >  2:499349/1‑150 (MQ=255)
                                                                                                                                    gCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACttt                  >  1:1408021/1‑150 (MQ=255)
                                                                                                                                        ggAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACTTTTTCa              <  1:398282/150‑1 (MQ=255)
                                                                                                                                        ggAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACTTTTTCa              <  1:826607/150‑1 (MQ=255)
                                                                                                                                         gAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCt                                                                                        >  1:4795930/1‑75 (MQ=255)
                                                                                                                                         gAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACTTTTTCAg             >  1:2331587/1‑150 (MQ=255)
                                                                                                                                           aaaCCCCGCCGAAGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACTTTTTCAGGt           >  1:2134706/1‑150 (MQ=255)
                                                                                                                                              ccccGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAATAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACTTTTTCAGGTTTa        <  2:4648307/150‑1 (MQ=255)
                                                                                                                                              ccccGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACTTTTTCAGGTTTa        <  1:3486903/150‑1 (MQ=255)
                                                                                                                                               cccGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACTTTTTCAGGTTTAc       >  2:3072207/1‑150 (MQ=255)
                                                                                                                                                    ttCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACTTTTTCAGGTTTACCAttt  >  2:458490/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
GTGAAAGCTTCCCGTTCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACTTTTTCAGGTTTACCATTT  >  NC_002947/1492125‑1492422

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