Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0123061,091,9920TG36.7% 373.7 / 54.6 196intergenic (+38/‑145)narI/ychSnitrate reductase 1, gamma (cytochrome b(NR)) subunit/predicted protein
Reads supporting (aligned to +/- strand):  ref base T (35/89);  new base G (72/0);  total (107/89)
Fisher's exact test for biased strand distribution p-value = 3.64e-27
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 218 total aligned reads displayed.

TGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGT  >  AP012306/1091845‑1092141
                                                                                                                                                    |                                                                                                                                                     
tGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGttt                                                                                                                                                     >  1:2730478/1‑150 (MQ=255)
tGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGttt                                                                                                                                                     <  2:269139/150‑1 (MQ=255)
 ggggATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGtttt                                                                                                                                                    >  1:3268715/1‑150 (MQ=255)
 ggggATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGtttt                                                                                                                                                    <  1:1477783/150‑1 (MQ=255)
 ggggATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGgttt                                                                                                                                                    >  2:2888707/1‑150 (MQ=255)
  gggATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGttttt                                                                                                                                                   <  1:860101/150‑1 (MQ=255)
  gggATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGttttt                                                                                                                                                   <  2:3566079/150‑1 (MQ=255)
   ggATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGgttttt                                                                                                                                                  >  1:2314385/1‑150 (MQ=255)
   ggATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGgttttt                                                                                                                                                  >  1:402403/1‑150 (MQ=255)
    gATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGgtttttt                                                                                                                                                 >  1:3139955/1‑150 (MQ=255)
      tGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTAt                                                                                                                                               <  1:1137252/150‑1 (MQ=255)
       gACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTAtt                                                                                                                                              >  1:809420/1‑149 (MQ=255)
       gACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATg                                                                                                                                              >  2:1515565/1‑150 (MQ=255)
          gTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTTTTggg                                                                                                                                            >  2:1721765/1‑150 (MQ=255)
          gTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGc                                                                                                                                           <  2:768004/150‑1 (MQ=255)
           ttATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATTGGCa                                                                                                                                          >  1:1020185/1‑150 (MQ=255)
            tATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCAc                                                                                                                                         >  2:3977647/1‑150 (MQ=255)
            tATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACACTGCTTCGGC‑GGGGGTTTTTTATGGGCAc                                                                                                                                         >  1:2040598/1‑150 (MQ=255)
             aTTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTTTTGGGACg                                                                                                                                        >  2:95694/1‑150 (MQ=255)
             aTTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTTTTGGCACg                                                                                                                                        >  2:944377/1‑150 (MQ=255)
              ttCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACgg                                                                                                                                       <  1:3952742/150‑1 (MQ=255)
                ccTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGGACgggg                                                                                                                                     >  2:446886/1‑148 (MQ=255)
                 ctgctgTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGc                                                                                                                                    >  2:1077335/1‑150 (MQ=255)
                  tgctgTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCg                                                                                                                                   <  1:3243317/150‑1 (MQ=255)
                  tgctgTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCg                                                                                                                                   >  1:405990/1‑150 (MQ=255)
                    ctgTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCggg                                                                                                                                 <  2:209088/150‑1 (MQ=255)
                        tCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGa                                                                                                                             <  2:2072036/150‑1 (MQ=255)
                          ccGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGt                                                                                                                           <  1:1198530/150‑1 (MQ=255)
                             ttCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTTTGGGCACGGTGCGGGGGGAGTTGt                                                                                                                        >  2:2154391/1‑150 (MQ=255)
                               ctcGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCg                                                                                                                      >  1:3377782/1‑150 (MQ=255)
                                 cgcgTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTTTGGGCACGGGGCGGGGGGGGGTGTCGGa                                                                                                                    >  2:1965955/1‑150 (MQ=255)
                                 cgcgTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGGTGTCGGa                                                                                                                    >  1:356539/1‑150 (MQ=255)
                                  gcgTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGAt                                                                                                                   <  1:1590939/150‑1 (MQ=255)
                                       tGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGcgct                                                                                                              >  1:2808651/1‑150 (MQ=255)
                                        gATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGcgctt                                                                                                             <  2:1431419/150‑1 (MQ=255)
                                          tACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGcgcttcg                                                                                                           >  2:2176380/1‑150 (MQ=255)
                                             acaTCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGcgcttcgctt                                                                                                        <  1:631470/150‑1 (MQ=255)
                                               aTCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTAt                                                                                                      <  1:713077/150‑1 (MQ=255)
                                               aTCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTAt                                                                                                      <  2:1946356/150‑1 (MQ=255)
                                               aTCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGcgcttcgctttt                                                                                                      >  2:818120/1‑148 (MQ=255)
                                                    gAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTTTCCGAc                                                                                                 >  2:753104/1‑150 (MQ=255)
                                                       cGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTTTTGGCACGGTGCGGGGGGAGGTGTCGGGTGCGCTTCGCTTATTCGAccca                                                                                              >  2:3068451/1‑148 (MQ=255)
                                                        gTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTAc                                                                                             <  1:3859628/150‑1 (MQ=255)
                                                        gTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTAc                                                                                             <  2:64714/150‑1 (MQ=255)
                                                        gTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCCCGGTGCGGGGGGAGTTGTCGGATGCGCTTCCCTTATCCCACCTAc                                                                                             >  2:2981284/1‑150 (MQ=255)
                                                          aCCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAg                                                                                           >  2:2898946/1‑150 (MQ=255)
                                                               tGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGggag                                                                                      >  1:3720875/1‑150 (MQ=255)
                                                                ggAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGggagg                                                                                     <  1:2364319/150‑1 (MQ=255)
                                                                  aGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGat                                                                                   >  1:2942567/1‑150 (MQ=255)
                                                                   gTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGata                                                                                  <  1:148950/150‑1 (MQ=255)
                                                                   gTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGata                                                                                  <  1:181541/150‑1 (MQ=255)
                                                                         gACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAg                                                                            >  1:500899/1‑150 (MQ=255)
                                                                         gACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAg                                                                            >  1:2789773/1‑150 (MQ=255)
                                                                            acGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGGTATTGTAGGcc                                                                         >  2:1028468/1‑150 (MQ=255)
                                                                                aaGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGt                                                                     <  1:3817025/150‑1 (MQ=255)
                                                                                 aGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTa                                                                    <  2:2250160/150‑1 (MQ=255)
                                                                                 aGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGGGCGGGGGGGGATGTCGGATGCGCTTCGCTTATCCCACCTACAGGGGAGGGTATTGTAGGCCCGGTa                                                                    >  2:1353808/1‑150 (MQ=255)
                                                                                   tACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTaaa                                                                  <  1:1226967/150‑1 (MQ=255)
                                                                                      cAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACgt                                                               <  1:3409881/150‑1 (MQ=255)
                                                                                      cAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACgt                                                               <  2:2951236/150‑1 (MQ=255)
                                                                                       aGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGATATTGTAGGCCCGGTTAACgtg                                                              >  1:1853581/1‑150 (MQ=255)
                                                                                         cTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACgtggt                                                            <  1:153780/150‑1 (MQ=255)
                                                                                               cgcgcTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATc                                                      <  1:2007615/150‑1 (MQ=255)
                                                                                               cgcgcTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATc                                                      <  2:885268/150‑1 (MQ=255)
                                                                                                 cgcTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCgg                                                    <  1:109615/150‑1 (MQ=255)
                                                                                                  gcTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCggg                                                   >  1:1389487/1‑150 (MQ=255)
                                                                                                  gcTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCggg                                                   >  1:2134848/1‑150 (MQ=255)
                                                                                                    tCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCa                                                 <  1:3203160/150‑1 (MQ=255)
                                                                                                        cACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCCTCGGGCAAAAt                                             >  1:762633/1‑150 (MQ=255)
                                                                                                        cACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAAt                                             >  1:3779099/1‑150 (MQ=255)
                                                                                                          cTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTTTGGGCACGGGGCGGGGGGAGTTGTCGGATGCGCTTCGCTTTTCCGACCTACAGGGGAGGGTATTGTAGGCCCGGTAAACATGGTGCCATCGGGGAAAATTa                                           >  2:2789774/1‑150 (MQ=255)
                                                                                                            aaGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGGTATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGc                                         >  1:899177/1‑150 (MQ=255)
                                                                                                            aaGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGGTATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGc                                         >  2:2480471/1‑150 (MQ=255)
                                                                                                                gAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCgg                                     <  1:2862107/150‑1 (MQ=255)
                                                                                                                 aaTTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGa                                    >  1:1312268/1‑150 (MQ=255)
                                                                                                                    tttAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTTTGGGCCCGGGGGGGGGGGTGATGTCGGGTGCGCTTTGCTTATCCGACCTACAGGGGGGGAGATTGTAGGCCCGGTAAACGTGGGGCCATCCGGGAAAAATAGCTCGGAgga                                 >  2:2834283/1‑148 (MQ=255)
                                                                                                                     ttAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTaa                                <  1:2535628/150‑1 (MQ=255)
                                                                                                                      tAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAAc                               <  2:1985725/150‑1 (MQ=255)
                                                                                                                      tAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGGGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGATATTGTAGGCCCGGTAAACGTGGGGCCATCGGGCAAAATTAGCTCGGGGGAAc                               >  2:848313/1‑150 (MQ=255)
                                                                                                                       aGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACa                              <  2:1580589/150‑1 (MQ=255)
                                                                                                                       aGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTTTGGGCACGGGGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGGTATTGTAGGCCCGGGAAACGTGGGGCCATCGGGCAAAAATAGCTCGGGGGAACa                              >  2:3243943/1‑150 (MQ=255)
                                                                                                                        gTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAg                             <  1:1712629/150‑1 (MQ=255)
                                                                                                                        gTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAg                             <  2:3648794/150‑1 (MQ=255)
                                                                                                                         ttCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAgg                            <  2:2853306/150‑1 (MQ=255)
                                                                                                                           cacaTAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGtt                          >  1:1451934/1‑150 (MQ=255)
                                                                                                                              aTAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGGACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGATATTGTAGGCCCGGGAAACGTGGGGGCATCGGGCAAAATTAGCTCGGAGGAACAGGTTTTg                       >  1:2455828/1‑150 (MQ=255)
                                                                                                                                aGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGGAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGGAACAGGGTTTGat                     >  1:509743/1‑150 (MQ=255)
                                                                                                                                 gACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGGGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGATATTGTAGGCCCGGGAAACGTGGTGCCATCGGGCAAAATTAGCTCGGGGGAACAGGTTTTGatt                    >  2:2425279/1‑150 (MQ=255)
                                                                                                                                  aCCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGatta                   <  1:694225/150‑1 (MQ=255)
                                                                                                                                    ccTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGattatt                 <  2:3342389/150‑1 (MQ=255)
                                                                                                                                     cTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTAttt                <  1:111991/150‑1 (MQ=255)
                                                                                                                                      tGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTa               <  2:2808561/150‑1 (MQ=255)
                                                                                                                                        cTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAg             <  2:508537/150‑1 (MQ=255)
                                                                                                                                          tCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGa           <  2:3499586/150‑1 (MQ=255)
                                                                                                                                             gc‑ggGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGc        >  2:2349855/1‑150 (MQ=255)
                                                                                                                                              c‑ggGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCg       >  2:1452055/1‑150 (MQ=255)
                                                                                                                                              c‑ggGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCg       >  2:1175866/1‑150 (MQ=255)
                                                                                                                                               ccggggTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCg       <  1:3650864/148‑1 (MQ=255)
                                                                                                                                                ggggTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGa      >  2:3178005/1‑150 (MQ=255)
                                                                                                                                                    tttttttATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGAtggt  >  2:1250057/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGT  >  AP012306/1091845‑1092141

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑