Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0123061,091,9920TG37.7% 580.2 / 66.6 301intergenic (+38/‑145)narI/ychSnitrate reductase 1, gamma (cytochrome b(NR)) subunit/predicted protein
Reads supporting (aligned to +/- strand):  ref base T (54/133);  new base G (113/0);  total (168/133)
Fisher's exact test for biased strand distribution p-value = 3.44e-41
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 321 total aligned reads displayed.

GCTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGT  >  AP012306/1091843‑1092141
                                                                                                                                                      |                                                                                                                                                     
gCTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGt                                                                                                                                                       >  2:1093224/1‑150 (MQ=255)
 cTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGtt                                                                                                                                                      <  2:3461376/150‑1 (MQ=255)
 cTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGtt                                                                                                                                                      <  1:1319052/150‑1 (MQ=255)
   ggggATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGtttt                                                                                                                                                    <  1:2513757/150‑1 (MQ=255)
       aTGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTa                                                                                                                                                <  1:1922016/150‑1 (MQ=255)
       aTGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTAGCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTa                                                                                                                                                <  1:3518332/150‑1 (MQ=255)
         gACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATg                                                                                                                                              <  1:1564251/150‑1 (MQ=255)
          aCGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATgg                                                                                                                                             >  1:2987847/1‑150 (MQ=255)
           cGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTTTggg                                                                                                                                            >  2:2153366/1‑150 (MQ=255)
            gTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTTTTggg                                                                                                                                            >  2:3127926/1‑150 (MQ=255)
               aTTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTTTGGGCACg                                                                                                                                        >  2:28037/1‑150 (MQ=255)
                   ctgctgTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGc                                                                                                                                    <  2:910092/150‑1 (MQ=255)
                    tgctgTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCg                                                                                                                                   <  1:3632099/150‑1 (MQ=255)
                    tgctgTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCg                                                                                                                                   >  1:529859/1‑150 (MQ=255)
                    tgctgTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCg                                                                                                                                   >  2:3954180/1‑150 (MQ=255)
                      ctgTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCCCgggggggg                                                                                                                                 >  1:2044116/1‑150 (MQ=255)
                      ctgTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCggg                                                                                                                                 >  1:1037662/1‑150 (MQ=255)
                        gTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGt                                                                                                                               >  1:662500/1‑150 (MQ=255)
                            ccGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGt                                                                                                                           >  1:3577010/1‑150 (MQ=255)
                             cGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGtt                                                                                                                          <  1:3911282/150‑1 (MQ=255)
                               ttCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGt                                                                                                                        <  1:4081853/150‑1 (MQ=255)
                                tctcGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTc                                                                                                                       <  1:3753490/150‑1 (MQ=255)
                                    gcgTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGAt                                                                                                                   <  1:2344745/150‑1 (MQ=255)
                                       tCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGGGCGGGGGGAGTTGTCGGATgcg                                                                                                                >  2:502795/1‑150 (MQ=255)
                                            tACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCCCGGTGCGGGGTGAGGTGTCGGATGcgcttcg                                                                                                           >  1:623546/1‑150 (MQ=255)
                                             acacATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGcgcttcgc                                                                                                          >  1:1130407/1‑150 (MQ=255)
                                              cacaTCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGCTTTTTTATTGGCACGGTGCGGGGGGAGTTGTCGGGTGcgcttcgct                                                                                                         >  1:3213985/1‑150 (MQ=255)
                                                  tCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATTGGCACGGGGCGGGGGGAGTTGTCGGATGCGCTTCGCTTTTc                                                                                                     >  2:1717830/1‑150 (MQ=255)
                                                   cTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATcc                                                                                                    <  1:212182/150‑1 (MQ=255)
                                                    tGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCg                                                                                                   <  2:4242056/150‑1 (MQ=255)
                                                     ggAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGa                                                                                                  >  1:2594556/1‑150 (MQ=255)
                                                       aGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGAcc                                                                                                >  1:2322090/1‑150 (MQ=255)
                                                       aGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGGGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGAcc                                                                                                >  1:3368569/1‑150 (MQ=255)
                                                        gCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGGACGGGGGGGGGGGAGTTGTCGGGTGGGCTTCGCTTATCCGACCt                                                                                               >  2:2255812/1‑150 (MQ=255)
                                                        gCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCt                                                                                               <  1:291736/150‑1 (MQ=255)
                                                        gCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCt                                                                                               <  1:921057/150‑1 (MQ=255)
                                                        gCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTTTTGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCt                                                                                               >  2:2879370/1‑150 (MQ=255)
                                                         cGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTa                                                                                              >  2:3792181/1‑150 (MQ=255)
                                                         cGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTa                                                                                              >  1:3222594/1‑150 (MQ=255)
                                                          gTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTAc                                                                                             >  1:2273747/1‑150 (MQ=255)
                                                           tACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATTGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACa                                                                                            >  1:3016393/1‑150 (MQ=255)
                                                           tACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACa                                                                                            >  2:545622/1‑150 (MQ=255)
                                                              cGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAggg                                                                                         <  1:1463754/150‑1 (MQ=255)
                                                               ggTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTTTTGGCACGGGGCGGGGGGAGGTGTCGGGTGGGCTTTCCTTTTCCCACCCACAgggg                                                                                        >  2:1637145/1‑150 (MQ=255)
                                                                  ggAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGggagg                                                                                     <  1:833090/150‑1 (MQ=255)
                                                                        tCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTg                                                                               <  1:1175699/150‑1 (MQ=255)
                                                                           gACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAg                                                                            <  1:4240417/150‑1 (MQ=255)
                                                                            acacGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAgg                                                                           >  1:3297442/1‑150 (MQ=255)
                                                                               cGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGccc                                                                        <  1:2368968/150‑1 (MQ=255)
                                                                                 tAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCgg                                                                      >  1:865554/1‑150 (MQ=255)
                                                                                  aaGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGt                                                                     >  1:1418073/1‑150 (MQ=255)
                                                                                      aCCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGTTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAAc                                                                 <  2:3700348/150‑1 (MQ=255)
                                                                                      aCCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAAc                                                                 <  2:2204506/150‑1 (MQ=255)
                                                                                        cAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACgt                                                               <  1:994248/150‑1 (MQ=255)
                                                                                        cAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGGTATTGTAGGCCCGGGAAACgt                                                               >  1:2486746/1‑150 (MQ=255)
                                                                                           cTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACgtggt                                                            >  1:3804026/1‑150 (MQ=255)
                                                                                             ggTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGc                                                          >  1:472391/1‑150 (MQ=255)
                                                                                                 cgcgcTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATc                                                      <  2:2755423/150‑1 (MQ=255)
                                                                                                  gcgcTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCg                                                     <  1:3978694/150‑1 (MQ=255)
                                                                                                  gcgcTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGGTATTGTAGGCCCGGGAAACGTGGTGCCATCg                                                     >  1:3460026/1‑150 (MQ=255)
                                                                                                    gcTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTTTGGGGGCGGTGGGGGGTGGGTTGTCGGATGCGCTTCGCCTATACCACCTAAAGGGGGGGGTGTTGTGGGGCCGGTGAACGTGGGGGCCTCggg                                                   >  2:4270626/1‑150 (MQ=255)
                                                                                                     cTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGc                                                  <  2:342961/150‑1 (MQ=255)
                                                                                                     cTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGGTATTGTAGGCCCGGTAAACGTGGGGCCATCGGGc                                                  >  2:174066/1‑150 (MQ=255)
                                                                                                      tCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCa                                                 <  2:2489988/150‑1 (MQ=255)
                                                                                                       cGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCaa                                                <  2:3327984/150‑1 (MQ=255)
                                                                                                         tCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCaaaa                                              <  1:2698867/150‑1 (MQ=255)
                                                                                                         tCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCAGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCaaaa                                              >  1:2831362/1‑150 (MQ=255)
                                                                                                          cACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTTTCCGACCTACAGGGGGGGATATTTTAGGCCCGGTAAACGTGGTGCCATCGGGCCAAAt                                             >  2:188479/1‑150 (MQ=255)
                                                                                                           aCTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAAtt                                            <  1:1188163/150‑1 (MQ=255)
                                                                                                            cTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTa                                           <  2:3353677/150‑1 (MQ=255)
                                                                                                            cTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGGGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGGAAACGTGGTGCCATCGGGCAAAATTa                                           >  1:1081505/1‑150 (MQ=255)
                                                                                                               aGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTTTAGGCCCGGTAAACGTGGGGCCATCGGGCAAAATTAGCt                                        >  2:2364133/1‑150 (MQ=255)
                                                                                                                   aaTTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTTTGGGCACGGTGCGGGGGGAGTTGTCGGGTGCGCTTCGCTTATCCGACCTACAGGGGAGGGTATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCggg                                    >  2:4034887/1‑149 (MQ=255)
                                                                                                                    aTTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAg                                   >  1:766061/1‑150 (MQ=255)
                                                                                                                     ttttAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGt                                  <  1:505824/150‑1 (MQ=255)
                                                                                                                      tttAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTTTGGGCACGGGGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGGTATTGTAGGCCCGGGAAACGTGGTGCCCTCGGGCAAAATTAGCTCggggga                                 >  2:4134578/1‑146 (MQ=255)
                                                                                                                        tAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAAc                               <  2:1577231/150‑1 (MQ=255)
                                                                                                                           ttCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTTTGGGCACGGTGGGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGGTATTGTTGGCCCGGCAAAAGAGGTGCCATCGGGCAAAATTAGCTCGGGGTAACAgg                            >  2:3490449/1‑150 (MQ=255)
                                                                                                                            tCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTTTGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCCCTTATCCGACCTACAGGGGAGGGTATTGTAGGCCCGGTAAACGTGGTGCCATCGGGGAAAATTTGGCCGGAGGAAAAggg                           >  2:3995248/1‑149 (MQ=255)
                                                                                                                             cacaTAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGtt                          <  1:2744650/150‑1 (MQ=255)
                                                                                                                              acaTAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGGGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGATATTGTAGGCCCGGTAAACGTGGGGCCATCGGGCAAAATTAGCTCGGGGTAACAGGgtt                         >  2:1791617/1‑150 (MQ=255)
                                                                                                                              acaTAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGGTATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGctt                         >  1:439900/1‑150 (MQ=255)
                                                                                                                                aTAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGGGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTg                       >  1:3746479/1‑150 (MQ=255)
                                                                                                                                 tAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGa                      <  1:1752152/150‑1 (MQ=255)
                                                                                                                                 tAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGa                      <  1:2128923/150‑1 (MQ=255)
                                                                                                                                 tAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGGGCGGGGGGAGTTGTCGGGTGCGCTTCGCTTATCCGACCTACAGGGGGGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCCAAATTAGCTCGGGGGAACAGGGTTTGa                      >  2:2104169/1‑150 (MQ=255)
                                                                                                                                   gACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGGGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGGTTTGatt                    >  1:2237853/1‑150 (MQ=255)
                                                                                                                                     cccTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGattat                  >  1:151991/1‑150 (MQ=255)
                                                                                                                                     cccTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGattat                  <  2:1319182/150‑1 (MQ=255)
                                                                                                                                         gCTTCGGC‑GGGGGTTTTTTTTGGGCACGGTGCGGGGGGAGTTGTCGGGTGCGCTTCGCTTTTCCGACCTACAGGGGGGGGTATTGTAGGCCCGGTAAACGTGGGGCCATCGGGGCAAATTAGCTCGGGGGAACAGGGTTTTATTAtttta              >  2:4108761/1‑148 (MQ=255)
                                                                                                                                          cTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAg             <  1:4228817/150‑1 (MQ=255)
                                                                                                                                          cTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAg             <  1:3876889/150‑1 (MQ=255)
                                                                                                                                             ttggaggggGTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGa           <  2:118816/145‑1 (MQ=255)
                                                                                                                                              ggc‑ggGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAg         >  2:263577/1‑150 (MQ=255)
                                                                                                                                               gc‑ggGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGc        >  2:4042845/1‑150 (MQ=255)
                                                                                                                                                c‑ggGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCg       >  2:2511924/1‑150 (MQ=255)
                                                                                                                                                   gggTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGAt     >  2:869842/1‑150 (MQ=255)
                                                                                                                                                     gTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGAtgg   <  2:243959/150‑1 (MQ=255)
                                                                                                                                                     gTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGAtgg   <  2:3618689/150‑1 (MQ=255)
                                                                                                                                                      tttttttATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGAtggt  <  2:740166/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GCTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGT  >  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: ‑