Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0123061,091,9920TG41.1% 411.6 / 79.9 255intergenic (+38/‑145)narI/ychSnitrate reductase 1, gamma (cytochrome b(NR)) subunit/predicted protein
Reads supporting (aligned to +/- strand):  ref base T (48/101);  new base G (104/0);  total (154/101)
Fisher's exact test for biased strand distribution p-value = 8.32e-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 286 total aligned reads displayed.

GCTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGT  >  AP012306/1091843‑1092141
                                                                                                                                                     |                                                                                                                                                     
gCTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGt                                                                                                                                                       <  2:2053028/150‑1 (MQ=255)
   ggggATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGgttt                                                                                                                                                    >  2:887592/1‑150 (MQ=255)
   ggggATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGgttt                                                                                                                                                    >  2:1777324/1‑150 (MQ=255)
   ggggATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGgttt                                                                                                                                                    >  1:2170744/1‑150 (MQ=255)
    gggATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGgtttt                                                                                                                                                   >  2:4071290/1‑150 (MQ=255)
     ggATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGtttttt                                                                                                                                                  <  2:3792200/150‑1 (MQ=255)
       aTGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTa                                                                                                                                                <  1:1031859/150‑1 (MQ=255)
       aTGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTa                                                                                                                                                >  2:2707258/1‑150 (MQ=255)
          aCGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTAttg                                                                                                                                             >  1:3092672/1‑148 (MQ=255)
            gTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGCTTTTTTTTgggg                                                                                                                                           >  2:3024917/1‑149 (MQ=255)
             ttATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATtggga                                                                                                                                          >  2:3141444/1‑148 (MQ=255)
              tATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCAc                                                                                                                                         >  1:1797047/1‑150 (MQ=255)
               aTTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTTTAGGCACg                                                                                                                                        >  1:795579/1‑150 (MQ=255)
                ttCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTTTGGGCACgg                                                                                                                                       >  2:2621520/1‑150 (MQ=255)
                  ccTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTg                                                                                                                                     <  1:3246899/150‑1 (MQ=255)
                  ccTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTTTTGGGACCGTg                                                                                                                                     >  2:2083448/1‑150 (MQ=255)
                     gctgTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTTTGGGCACGGTGCgg                                                                                                                                  >  2:3501477/1‑150 (MQ=255)
                        gTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTTTGGGCGCCGTggggggg                                                                                                                               >  2:4248816/1‑149 (MQ=255)
                           cccGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAg                                                                                                                            <  1:2421038/150‑1 (MQ=255)
                              gTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGTGAGTTg                                                                                                                         >  1:890128/1‑150 (MQ=255)
                              gTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGGGCGGGGGGAGGTg                                                                                                                         >  1:1611416/1‑150 (MQ=255)
                                 ctcGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGGGAGGTGTCg                                                                                                                      >  1:2584414/1‑150 (MQ=255)
                                   cgcgTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGa                                                                                                                    >  1:439777/1‑150 (MQ=255)
                                    gcgTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGGACGGTGCGGGGGGGGTTGTCgggt                                                                                                                   >  2:357637/1‑148 (MQ=255)
                                      gTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATgc                                                                                                                 >  2:2874471/1‑150 (MQ=255)
                                      gTCTGATACACATCTGGAGCATACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATTGGCACGGTGCGGGGTGAGTTGTCGGATgc                                                                                                                 >  1:968605/1‑150 (MQ=255)
                                         tGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGcgct                                                                                                              <  1:708496/150‑1 (MQ=255)
                                          gATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGcgctt                                                                                                             >  1:14745/1‑150 (MQ=255)
                                           aTACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGGGCGGGGGGAGGTGTCGGATGcgcttc                                                                                                            >  1:2386781/1‑150 (MQ=255)
                                               acaTCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGcgcttcgctt                                                                                                        <  1:3283473/150‑1 (MQ=255)
                                                caTCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTa                                                                                                       >  2:1583613/1‑150 (MQ=255)
                                                 aTCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATTGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTAt                                                                                                      >  2:1803677/1‑150 (MQ=255)
                                                    tGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCg                                                                                                   <  1:3801770/150‑1 (MQ=255)
                                                     ggAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGGGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGa                                                                                                  >  2:501159/1‑150 (MQ=255)
                                                       aGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGAcc                                                                                                <  1:2824109/150‑1 (MQ=255)
                                                         cGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTa                                                                                              <  2:2306004/150‑1 (MQ=255)
                                                         cGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTa                                                                                              <  1:524796/150‑1 (MQ=255)
                                                         cGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTa                                                                                              <  1:772433/150‑1 (MQ=255)
                                                          gTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGGGGGGTGTCGGGTGCGGTTCGCCTATCCGACCTAc                                                                                             >  2:482661/1‑150 (MQ=255)
                                                            aCCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAg                                                                                           <  2:1231711/150‑1 (MQ=255)
                                                            aCCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAg                                                                                           <  2:138308/150‑1 (MQ=255)
                                                              cGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAggg                                                                                         <  2:2819684/150‑1 (MQ=255)
                                                               ggTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGGTGCGCTTCGCTTATCCGACCTACAgggg                                                                                        >  2:1979404/1‑150 (MQ=255)
                                                                gTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGgga                                                                                       >  2:2638164/1‑150 (MQ=255)
                                                                    aGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGGTGCGCTTCGCTTATCCGACCTACAGgggggggt                                                                                   >  2:1363992/1‑148 (MQ=255)
                                                                     gTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGata                                                                                  <  2:4290395/150‑1 (MQ=255)
                                                                      tATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGGTGCGCCTCGCTTATTCGACCTACAGGGGGGGatat                                                                                 >  2:600709/1‑150 (MQ=255)
                                                                        tCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTg                                                                               <  2:1507714/150‑1 (MQ=255)
                                                                             cacGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGc                                                                          <  1:1077819/150‑1 (MQ=255)
                                                                             cacGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGGGAGGTGTCGGGTGCGCTTCGCTTATTCGACCTACAGGGGAGGGGATTGTAggg                                                                          >  2:373336/1‑149 (MQ=255)
                                                                              acGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGcc                                                                         <  1:3901725/150‑1 (MQ=255)
                                                                                gTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCg                                                                       <  2:3656528/150‑1 (MQ=255)
                                                                                gTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCg                                                                       <  1:47635/150‑1 (MQ=255)
                                                                                  aaGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGt                                                                     <  2:793806/150‑1 (MQ=255)
                                                                                    gTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTTTGGGCACGGGGGGGGGGGAGTTGTCGGATGCGCTTCGCTTTTCCGACCTACAGGGGGGGATATTTTAGGCCCGGTaa                                                                   >  2:4130478/1‑150 (MQ=255)
                                                                                       ccAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGGTATTGTAGGCCCGGTAAACg                                                                >  1:3716573/1‑150 (MQ=255)
                                                                                       ccAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGGGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGATATTGTAGGCCCGGTAAACg                                                                >  2:1902460/1‑150 (MQ=255)
                                                                                        cAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACgt                                                               <  1:1229774/150‑1 (MQ=255)
                                                                                         aGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGGGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACgtg                                                              >  1:2330112/1‑150 (MQ=255)
                                                                                            tGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACgtggtg                                                           <  2:66377/150‑1 (MQ=255)
                                                                                              gTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGcc                                                         <  2:2997374/150‑1 (MQ=255)
                                                                                               tGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCa                                                        <  2:1856887/150‑1 (MQ=255)
                                                                                               tGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCa                                                        <  2:4221169/150‑1 (MQ=255)
                                                                                                gcgcgcTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCAt                                                       <  1:4022707/150‑1 (MQ=255)
                                                                                                  gcgcTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCg                                                     <  2:2114493/150‑1 (MQ=255)
                                                                                                    gcTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGGTATTGTAGGCCCGGTAAACGTGGTGCCATCggg                                                   >  2:2927598/1‑150 (MQ=255)
                                                                                                     cTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGc                                                  >  1:26979/1‑150 (MQ=255)
                                                                                                     cTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGATATTGTAGGCCCGGGAAACGTGGGGCCATCGGGc                                                  >  2:116192/1‑150 (MQ=255)
                                                                                                           aCTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAAtt                                            <  1:1703111/150‑1 (MQ=255)
                                                                                                            cTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGGTATTGTAGGCCCGGGAAACGTGGTGCCATCGGGCAAAATTa                                           >  2:1633535/1‑150 (MQ=255)
                                                                                                               aGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCt                                        <  1:1277413/150‑1 (MQ=255)
                                                                                                                  gAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCCACCTACAGGGGGGGATATTGTAGGCCCGGGAAACGTGGTGCCATCCGGCCAAATTAGCTCgg                                     >  2:2781106/1‑150 (MQ=255)
                                                                                                                   aaTTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGGTATTGTAGGCCCGGGAAACGTGGTGCCATCGGGCAAAATTAGCTCGGa                                    >  2:2555183/1‑150 (MQ=255)
                                                                                                                     ttttAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGt                                  <  1:398827/150‑1 (MQ=255)
                                                                                                                     ttttAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGt                                  <  1:64555/150‑1 (MQ=255)
                                                                                                                      tttAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTa                                 <  1:4201122/150‑1 (MQ=255)
                                                                                                                       ttAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTaa                                <  2:3221345/150‑1 (MQ=255)
                                                                                                                        tAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAAc                               <  1:583701/150‑1 (MQ=255)
                                                                                                                             cacaTAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGtt                          >  1:3995041/1‑150 (MQ=255)
                                                                                                                               caTAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGtttt                        <  2:2275531/150‑1 (MQ=255)
                                                                                                                               caTAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGtttt                        <  1:3327464/150‑1 (MQ=255)
                                                                                                                               caTAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGtttt                        <  1:1403399/150‑1 (MQ=255)
                                                                                                                               caTAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGGTATTGTAGGCCCGGTAAACGTGGGGCCATCGGGCAAAATTAGCTCGGAGTAACAGGgttt                        >  1:3416621/1‑150 (MQ=255)
                                                                                                                                aTAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGGGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTg                       >  1:224250/1‑150 (MQ=255)
                                                                                                                                  aGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGGGCCATCGGGCAAAATTAGCTCGGAGTAACAGGGTTTGat                     >  2:2439777/1‑150 (MQ=255)
                                                                                                                                   gACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGGTGCGCTTCGCTTATCCGACCTACAGGGGGGGATATTGTAGGCCCGGTAAAAGTGGTGCCATCGGGCAAAATTAGCTCGGAGGAACAGGGTTTGatt                    >  2:2470934/1‑150 (MQ=255)
                                                                                                                                     cccTGCTTCGGCGGGGTTTTTTTATGGGGACGGTGGGGGGTGAGTTGTCGGATGCGGTTCGCTTATCCGACCTACCGGGGGGGGTATTGTAGGCCCGGTAAACGTGGGGCCATCGGGCAAAAGTAGCTCGGGGTAACAGGTTTTGattat                  >  2:1270650/1‑150 (MQ=255)
                                                                                                                                     cccTGCTTCGGCGGGGGTTTTTTTTGGGCACGGGGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGGTATTGTAGGCCCGGTAAACGTGGGGCCATCGGGGAAAATTAGCTCGGGGGAAAAGGGTTTGatttt                  >  2:3411643/1‑148 (MQ=255)
                                                                                                                                      ccTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAAATAGCTCGGAGGAACAGGGTTTGattatt                 >  2:3062744/1‑150 (MQ=255)
                                                                                                                                      ccTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGGAACAGGTTTTGattatt                 >  1:2792837/1‑150 (MQ=255)
                                                                                                                                         gCTTCGGCGGGGGTTTTTTTTGGGCACGGGGCGGGGGGAGGTGTCGGATGCGCTTCGCTTTTCCGACCCACAGGGGGGGATATTGTAGGCCCGGGAAAAGTGGGGCCCTCGGGCAAAATTAGCTCGGAGTAAAAGGTTTTGATTATTTaa              >  2:90803/1‑150 (MQ=255)
                                                                                                                                          cTTCGGCGGGGTTTTTTTATGGGCACGGGGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGGTTTGATTATTTAAg             >  1:1312725/1‑150 (MQ=255)
                                                                                                                                               gcggGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGc        >  1:410121/1‑150 (MQ=255)
                                                                                                                                               gcggGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGc        >  2:2565928/1‑150 (MQ=255)
                                                                                                                                                  gggTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGAt     <  1:919561/150‑1 (MQ=255)
                                                                                                                                                  gggTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGAt     >  2:3975149/1‑150 (MQ=255)
                                                                                                                                                  gggTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGAt     >  2:1134923/1‑150 (MQ=255)
                                                                                                                                                   ggTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGAtg    <  2:3899996/150‑1 (MQ=255)
                                                                                                                                                   ggTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGAtg    >  2:2187378/1‑150 (MQ=255)
                                                                                                                                                     tttttttATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGAtggt  <  1:202564/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGT  >  AP012306/1091843‑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: