Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0123061,091,9920TG45.1% 306.8 / 80.0 226intergenic (+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 (102/0);  total (137/89)
Fisher's exact test for biased strand distribution p-value = 2.83e-34
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 257 total aligned reads displayed.

GCTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGT  >  AP012306/1091843‑1092141
                                                                                                                                                     |                                                                                                                                                      
gCTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGt                                                                                                                                                        <  2:2458064/150‑1 (MQ=255)
 cTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGtt                                                                                                                                                       <  1:3185757/150‑1 (MQ=255)
 cTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGtt                                                                                                                                                       <  1:3332901/150‑1 (MQ=255)
   ggggATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGgttt                                                                                                                                                     >  2:3191009/1‑150 (MQ=255)
     ggATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGgttttt                                                                                                                                                   >  1:3855156/1‑150 (MQ=255)
     ggATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGgttttt                                                                                                                                                   >  1:1973249/1‑150 (MQ=255)
      gATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGttttttt                                                                                                                                                  <  1:2258621/150‑1 (MQ=255)
      gATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGgtttttt                                                                                                                                                  >  2:2650807/1‑150 (MQ=255)
        tGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGgtttttttt                                                                                                                                                >  2:3775658/1‑148 (MQ=255)
         gACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGgttttttttt                                                                                                                                               >  2:2106934/1‑147 (MQ=255)
         gACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATg                                                                                                                                               >  2:4143567/1‑150 (MQ=255)
          aCGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATgg                                                                                                                                              <  1:185573/150‑1 (MQ=255)
            gTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTTTGGGc                                                                                                                                            >  2:4277129/1‑150 (MQ=255)
             ttATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCa                                                                                                                                           <  1:2394297/150‑1 (MQ=255)
             ttATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATgggga                                                                                                                                           >  2:1457655/1‑148 (MQ=255)
             ttATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCa                                                                                                                                           >  2:2294369/1‑150 (MQ=255)
              tATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCAc                                                                                                                                          >  1:4003860/1‑150 (MQ=255)
                ttCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACgg                                                                                                                                        >  2:3847303/1‑150 (MQ=255)
                 tCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTTTTGGCACGGt                                                                                                                                       >  2:3555589/1‑150 (MQ=255)
                    tgctgTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCg                                                                                                                                    >  1:3914598/1‑150 (MQ=255)
                     gctgTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCgg                                                                                                                                   >  1:2061294/1‑150 (MQ=255)
                        gTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATTGGCACGGTGCggggg                                                                                                                                >  2:158196/1‑149 (MQ=255)
                           cccGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGTGAg                                                                                                                             >  1:217579/1‑150 (MQ=255)
                            ccGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTCCGGGGTGAgc                                                                                                                            >  1:4365004/1‑149 (MQ=255)
                            ccGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGTGAgg                                                                                                                            >  1:1311503/1‑149 (MQ=255)
                            ccGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGTGAGt                                                                                                                            >  1:2645135/1‑150 (MQ=255)
                            ccGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGTGAGt                                                                                                                            >  1:54776/1‑150 (MQ=255)
                              gTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGTGAGTTg                                                                                                                          >  1:1213535/1‑150 (MQ=255)
                                tctcGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTc                                                                                                                        <  1:4098497/150‑1 (MQ=255)
                                tctcGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTTTTGGCACGGGGCGGGGGGGGTTGTc                                                                                                                        >  2:4381795/1‑150 (MQ=255)
                                 ctcGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCg                                                                                                                       <  1:3541364/150‑1 (MQ=255)
                                  tcGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGGGCGGGGGGAGGTGTCgg                                                                                                                      >  2:2636761/1‑150 (MQ=255)
                                    gcgTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATTGGCACGGTGCGGGGTGAGTTGTCGGAt                                                                                                                    >  1:2592586/1‑150 (MQ=255)
                                     cgTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATg                                                                                                                   <  1:4283940/150‑1 (MQ=255)
                                            tACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGcgcttcg                                                                                                            <  2:3580651/150‑1 (MQ=255)
                                            tACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGCACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGcgcttcg                                                                                                            >  2:591280/1‑150 (MQ=255)
                                             acacATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGcgcttcgc                                                                                                           >  1:3444122/1‑150 (MQ=255)
                                                 aTCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTAt                                                                                                       <  1:2199314/150‑1 (MQ=255)
                                                    tGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCg                                                                                                    <  1:4130980/150‑1 (MQ=255)
                                                      gAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGAc                                                                                                  >  1:600601/1‑150 (MQ=255)
                                                         cGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTa                                                                                               >  1:921311/1‑150 (MQ=255)
                                                          gTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTAc                                                                                              <  1:1354942/150‑1 (MQ=255)
                                                                gTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGAGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGgga                                                                                        >  2:1142797/1‑150 (MQ=255)
                                                                     gTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGGACGGGGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGata                                                                                   >  2:1827561/1‑150 (MQ=255)
                                                                      tATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGatat                                                                                  >  1:2975280/1‑150 (MQ=255)
                                                                       aTCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGGGCGGGGGGAGGTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGGTAtt                                                                                 >  2:1354271/1‑150 (MQ=255)
                                                                          tGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑a                                                                             <  2:894046/150‑1 (MQ=255)
                                                                           gACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACGGGGGAGGATATTGT‑Ag                                                                            <  1:1597256/150‑1 (MQ=255)
                                                                            acacGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGGGCGGGGGGGGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGATATTGT‑Agg                                                                           >  2:4223367/1‑150 (MQ=255)
                                                                               cGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGccc                                                                        <  2:3695931/150‑1 (MQ=255)
                                                                               cGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTTTGGGCACGGTGCGGGGGGGGTTGTCGGGTGCGCTTCGCTTATCCGACCTACAGGGGAGGGTATTTT‑AGGgcc                                                                        >  2:2023112/1‑150 (MQ=255)
                                                                                 tAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAGTTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGGGCGGGGGGAGTTGTCGGATGAGCTTCTCTTATACGACCTACAGGGGGGGGTATTGT‑AGGCCCgg                                                                      >  2:2894564/1‑150 (MQ=255)
                                                                                 tAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCgg                                                                      <  1:4422527/150‑1 (MQ=255)
                                                                                 tAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCgg                                                                      <  2:2044650/150‑1 (MQ=255)
                                                                                  aaGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGGGCTTCGCTTATCCGACCTACAGGGGGGGGTATTGT‑AGGCCCggg                                                                     >  2:547548/1‑149 (MQ=255)
                                                                                     tACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATTGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGGTATTTT‑AGGCCCGGTaaa                                                                  >  2:2367142/1‑150 (MQ=255)
                                                                                      aCCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCGGGAAAc                                                                 >  1:971667/1‑150 (MQ=255)
                                                                                         aGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCGGTAAACgtg                                                              <  1:2732451/150‑1 (MQ=255)
                                                                                          gCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGGTATTGT‑AGGCCCGGTAAACgtgg                                                             >  2:1670618/1‑150 (MQ=255)
                                                                                           cTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCGGGAAACgtggt                                                            >  2:196540/1‑150 (MQ=255)
                                                                                            tGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCGGTAAACgtggtg                                                           <  1:2814803/150‑1 (MQ=255)
                                                                                              gTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGGTATTGT‑AGGGCCGGGAAACGTGGGGcc                                                         >  2:3074477/1‑150 (MQ=255)
                                                                                                gcgcgcTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCAACCTACAGGGGAGGATATTGT‑AGGCCCGGTAAACGTGGTGCCAt                                                       <  1:77956/150‑1 (MQ=255)
                                                                                                  gcgcTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCGGTAAACGTGGTGCCATCg                                                     <  2:710867/150‑1 (MQ=255)
                                                                                                   cgcTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCGGTAAACGTGGTGCCATCgg                                                    <  1:4037364/150‑1 (MQ=255)
                                                                                                   cgcTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCGGGAAACGTGGTGCCATCgg                                                    >  2:505095/1‑150 (MQ=255)
                                                                                                    gcTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCGGTAAACGTGGTGCCATCggg                                                   <  2:3677771/150‑1 (MQ=255)
                                                                                                    gcTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTTTGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGGTATTGT‑AGGCCCGGTAAACGTGGTGCCATCggg                                                   >  2:2964889/1‑150 (MQ=255)
                                                                                                    gcTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCGGTAAACGTGGGGCCATCggg                                                   >  1:3674466/1‑150 (MQ=255)
                                                                                                       cGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCGGGAAACGTGGTGCCATCGGGCaa                                                >  2:1700633/1‑150 (MQ=255)
                                                                                                        gTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTTTGGGGACGGGGCGGGGGGAGGTTTCGGGTGGGCTTCGCTTTTCCCACCTACAGGGGGGGGTATTGTAAGGCCCGGAAAAGGGGGGGCC‑TCGGGGaaa                                               >  2:1013711/1‑150 (MQ=255)
                                                                                                         tCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGGTATTGT‑AGGCCCGGTAAACGTGGGGCCATCGGGCaaaa                                              >  2:870368/1‑150 (MQ=255)
                                                                                                           aCTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCGGTAAACGTGGTGCCATCGGGCAAAAtt                                            <  2:1752592/150‑1 (MQ=255)
                                                                                                            cTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTa                                           <  2:1510514/150‑1 (MQ=255)
                                                                                                               aGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCt                                        <  2:3152702/150‑1 (MQ=255)
                                                                                                               aGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTTTGGGGACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGGTATTGT‑AGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCt                                        >  2:2525605/1‑150 (MQ=255)
                                                                                                                   aaTTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTTTGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGATATTGT‑AGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGa                                    >  2:2270554/1‑150 (MQ=255)
                                                                                                                   aaTTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGa                                    >  1:1163430/1‑150 (MQ=255)
                                                                                                                    aTTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAg                                   <  1:1535748/150‑1 (MQ=255)
                                                                                                                     ttttAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCGGTAAACGTGGGGCCATCGGGCAAAATTAGCTCGGAGt                                  >  1:1407626/1‑150 (MQ=255)
                                                                                                                      tttAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAgga                                 >  1:1723684/1‑148 (MQ=255)
                                                                                                                      tttAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGGTATTGT‑AGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTa                                 >  2:214763/1‑150 (MQ=255)
                                                                                                                          gTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAg                             >  1:6696/1‑150 (MQ=255)
                                                                                                                           ttCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAgg                            <  2:328380/150‑1 (MQ=255)
                                                                                                                           ttCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCCACCTACAGGGGAGGATATTGT‑AGGGCCGGGAAACGTGGGGCCATCGGGCAAAATTAGCTCGGGGGAACAgg                            >  2:961445/1‑150 (MQ=255)
                                                                                                                              acaTAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGttt                         <  1:3796083/150‑1 (MQ=255)
                                                                                                                               caTAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGtttt                        <  1:4247448/150‑1 (MQ=255)
                                                                                                                               caTAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGtttt                        <  1:3118554/150‑1 (MQ=255)
                                                                                                                               caTAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGtttt                        <  2:3064218/150‑1 (MQ=255)
                                                                                                                                aTAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGGTTTg                       >  1:2873400/1‑150 (MQ=255)
                                                                                                                                  aGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGat                     <  1:1799518/150‑1 (MQ=255)
                                                                                                                                      ccTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGattatt                 >  1:3811159/1‑150 (MQ=255)
                                                                                                                                        tGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTa               <  1:1001108/150‑1 (MQ=255)
                                                                                                                                           ttCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAgg            <  1:3483109/150‑1 (MQ=255)
                                                                                                                                                cggGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCg       <  2:3931195/150‑1 (MQ=255)
                                                                                                                                                cggGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCg       >  2:3367072/1‑150 (MQ=255)
                                                                                                                                                 ggggTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGa      <  1:618808/150‑1 (MQ=255)
                                                                                                                                                  gggTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGAt     <  2:4336157/150‑1 (MQ=255)
                                                                                                                                                  gggTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGAt     >  2:3429786/1‑150 (MQ=255)
                                                                                                                                                     tttttttATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGAtggt  >  2:3485826/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GCTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGT‑AGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGT  >  AP012306/1091843‑1092141

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