Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0123061,091,9920TG38.3% 449.9 / 62.1 269intergenic (+38/‑145)narI/ychSnitrate reductase 1, gamma (cytochrome b(NR)) subunit/predicted protein
Reads supporting (aligned to +/- strand):  ref base T (60/106);  new base G (103/0);  total (163/106)
Fisher's exact test for biased strand distribution p-value = 9.14e-32
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 291 total aligned reads displayed.

GCTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGT  >  AP012306/1091843‑1092141
                                                                                                                                                      |                                                                                                                                                     
gCTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGt                                                                                                                                                       <  1:2062763/150‑1 (MQ=255)
gCTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGt                                                                                                                                                       <  1:1988935/150‑1 (MQ=255)
  tGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGttt                                                                                                                                                     <  1:2598503/150‑1 (MQ=255)
   ggggATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGtttt                                                                                                                                                    <  1:3233034/150‑1 (MQ=255)
    gggATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGgtttt                                                                                                                                                   >  1:3379237/1‑150 (MQ=255)
    gggATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCATGCTTCGGC‑GGGGgtttt                                                                                                                                                   >  2:1093637/1‑150 (MQ=255)
     ggATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGgttttt                                                                                                                                                  >  2:171440/1‑150 (MQ=255)
        tGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTAt                                                                                                                                               <  2:2770039/150‑1 (MQ=255)
         gACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGtttttttggg                                                                                                                                              <  1:3851524/150‑4 (MQ=255)
         gACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATg                                                                                                                                              >  1:3960616/1‑150 (MQ=255)
           cGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATtgg                                                                                                                                            >  1:1880821/1‑150 (MQ=255)
            gTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATTGGc                                                                                                                                           >  1:3684747/1‑150 (MQ=255)
             ttATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTTTTggga                                                                                                                                           >  2:4389309/1‑149 (MQ=255)
             ttATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCa                                                                                                                                          <  2:2219340/150‑1 (MQ=255)
               aTTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTATATGGGCACg                                                                                                                                        >  1:4372706/1‑150 (MQ=255)
                   ctgctgTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGc                                                                                                                                    <  1:1043556/150‑1 (MQ=255)
                   ctgctgTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTTTTGGGACggggg                                                                                                                                    >  2:886102/1‑147 (MQ=255)
                     gctgTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCgg                                                                                                                                  <  2:3173954/150‑1 (MQ=255)
                      ctgTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCggg                                                                                                                                 <  1:2386523/150‑1 (MQ=255)
                         ttCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTg                                                                                                                              <  2:3928584/150‑1 (MQ=255)
                         ttCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGTg                                                                                                                              >  1:13469/1‑150 (MQ=255)
                             cGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGTGAGtt                                                                                                                          >  1:445756/1‑150 (MQ=255)
                              gTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTg                                                                                                                         >  1:2553616/1‑150 (MQ=255)
                               ttCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTTTGGGGACGGTGCGGGGGGAGTTGt                                                                                                                        >  2:1523/1‑150 (MQ=255)
                                 ctcGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTTTTGGGACGGGGCGGGGGGGGGTGTCg                                                                                                                      >  2:726051/1‑150 (MQ=255)
                                     cgTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATg                                                                                                                  >  1:1353421/1‑150 (MQ=255)
                                      gTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATgc                                                                                                                 >  1:2696412/1‑150 (MQ=255)
                                      gTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTTTGGGCACGGTGCGGGGGGAGTTGTCGGATgg                                                                                                                 >  1:4156775/1‑149 (MQ=255)
                                       tCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATgcg                                                                                                                >  1:1582983/1‑150 (MQ=255)
                                          gATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATTGGCACGGTGCGGGGTGAGGTGTCGGATGcgctt                                                                                                             >  1:3406157/1‑150 (MQ=255)
                                           aTACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGcgcttc                                                                                                            >  2:1317182/1‑150 (MQ=255)
                                            tACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGcgcttcg                                                                                                           <  1:2764950/150‑1 (MQ=255)
                                              cacaTCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGcgcttcgct                                                                                                         <  1:478906/150‑1 (MQ=255)
                                               acaTCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGcgcttcgctt                                                                                                        >  1:4205179/1‑150 (MQ=255)
                                                caTCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTa                                                                                                       <  2:1824304/150‑1 (MQ=255)
                                                 aTCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTAt                                                                                                      >  1:3038811/1‑150 (MQ=255)
                                                  tCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATTGGGACGGTGCGGGGGGAGTTGTCGGGTGCGCTTCGCTTTTc                                                                                                     >  2:193335/1‑150 (MQ=255)
                                                   cTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATcc                                                                                                    >  2:2626965/1‑150 (MQ=255)
                                                    tGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCg                                                                                                   <  2:3899319/150‑1 (MQ=255)
                                                          gTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTAc                                                                                             <  1:1926623/150‑1 (MQ=255)
                                                            aCCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGGGCGGGGGGAGGTGTCGGATGCGCTTCGCTTATCCGACCTACAg                                                                                           >  2:3623472/1‑150 (MQ=255)
                                                              cGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAggg                                                                                         <  2:4245771/150‑1 (MQ=255)
                                                                 tGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCCCGGTGGGGGGGGAGGTGTCGGATGCGCTTCGCTTATCCCACCTACAGGggag                                                                                      >  2:1650779/1‑150 (MQ=255)
                                                                 tGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGggag                                                                                      <  1:2908759/150‑1 (MQ=255)
                                                                     gTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGata                                                                                  >  2:2408687/1‑150 (MQ=255)
                                                                      tATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGatat                                                                                 <  1:316153/150‑1 (MQ=255)
                                                                      tATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGatat                                                                                 >  1:11176/1‑150 (MQ=255)
                                                                       aTCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCCTCGCTTATCCGGCCTACAGGGGGGGATAtt                                                                                >  2:2792493/1‑150 (MQ=255)
                                                                        tCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTg                                                                               >  1:2874959/1‑150 (MQ=255)
                                                                          tGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGGTATTGTa                                                                             >  2:1562778/1‑150 (MQ=255)
                                                                            acacGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAgg                                                                           <  2:3335319/150‑1 (MQ=255)
                                                                             cacGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGGACGGTGCGGGGGGAGGTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGATATTTTAGGc                                                                          >  2:1791048/1‑150 (MQ=255)
                                                                             cacGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGGTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGATATTGTAGGc                                                                          >  2:211713/1‑150 (MQ=255)
                                                                                gTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCg                                                                       <  1:1645034/150‑1 (MQ=255)
                                                                                  aaGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGATTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGt                                                                     <  1:4348396/150‑1 (MQ=255)
                                                                                          gCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACgtgg                                                             <  2:2148881/150‑1 (MQ=255)
                                                                                           cTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACgtggt                                                            >  1:1472981/1‑150 (MQ=255)
                                                                                            tGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACgtggtg                                                           <  2:800920/150‑1 (MQ=255)
                                                                                             ggTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGATATTGTAGGCCCGGTAAACGTGGTGc                                                          >  2:3493909/1‑150 (MQ=255)
                                                                                              gTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGcc                                                         <  2:283195/150‑1 (MQ=255)
                                                                                               tGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCa                                                        <  2:2741956/150‑1 (MQ=255)
                                                                                                gcgcgcTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCAt                                                       >  1:1767650/1‑150 (MQ=255)
                                                                                                    gcTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCggg                                                   <  1:2499484/150‑1 (MQ=255)
                                                                                                     cTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGc                                                  >  1:2831691/1‑150 (MQ=255)
                                                                                                     cTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGc                                                  >  2:140063/1‑150 (MQ=255)
                                                                                                       cGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCaa                                                >  1:611777/1‑150 (MQ=255)
                                                                                                        gTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCaaa                                               >  1:1555648/1‑150 (MQ=255)
                                                                                                         tCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGGTATTGTAGGCCCGGGAAACGTGGGGCCATCGGGGaaaa                                              >  2:3515642/1‑150 (MQ=255)
                                                                                                          cACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAAt                                             >  1:1179659/1‑150 (MQ=255)
                                                                                                            cTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTa                                           <  2:619803/150‑1 (MQ=255)
                                                                                                            cTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTa                                           <  2:1212761/150‑1 (MQ=255)
                                                                                                             tAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAg                                          <  1:913700/150‑1 (MQ=255)
                                                                                                             tAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAg                                          >  1:3556074/1‑150 (MQ=255)
                                                                                                              aaGCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGc                                         <  1:2281863/150‑1 (MQ=255)
                                                                                                                gCGAATTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTc                                       <  2:4206262/150‑1 (MQ=255)
                                                                                                                   aaTTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGa                                    >  1:797472/1‑150 (MQ=255)
                                                                                                                   aaTTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGa                                    <  1:329563/150‑1 (MQ=255)
                                                                                                                    aTTTTAGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGGGCCATCGGGCAAAATTAGCTCGGAg                                   >  2:4165177/1‑150 (MQ=255)
                                                                                                                     ttttAGTTCACATAGACCCTGCTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCggggg                                  >  1:3730807/1‑147 (MQ=255)
                                                                                                                         aGTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACa                              <  1:3432281/150‑1 (MQ=255)
                                                                                                                          gTTCACATAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAg                             <  1:3361536/150‑1 (MQ=255)
                                                                                                                             cacaTAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGtt                          <  1:564924/150‑1 (MQ=255)
                                                                                                                              acaTAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGttt                         <  2:2076017/150‑1 (MQ=255)
                                                                                                                               caTAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCATGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGtttt                        <  1:3587843/150‑1 (MQ=255)
                                                                                                                               caTAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGgttt                        >  1:2929406/1‑150 (MQ=255)
                                                                                                                                 tAGACCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGa                      <  2:1003775/150‑1 (MQ=255)
                                                                                                                                  aGACCCTGCTTCGGC‑GGGGGTTTTTTTTGGGGACGGGGCGGGGGGAGGTGTCGGGTGCGCTTCGCTTATCCGACCTACAGGGGGGGGTATTGTAGGCCCGGGAAAAGTGGGGGCCTCGGGCAAAATTAGCTCGGGGTAACAGGGTTTTat                     >  2:2038244/1‑150 (MQ=255)
                                                                                                                                    aCCCTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGGGCGGGGGGAGTTGTCGGGTGCGCTTCGCTTTTCCGACCTACAGGGGAGGGTATTGTAGGCCCGGGAAACGGGGTGCCATCGGGCAAAATTAGCTCGGAGGAACAGGTTTttatta                   >  2:2025748/1‑150 (MQ=255)
                                                                                                                                     cccTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGattat                  >  1:735180/1‑150 (MQ=255)
                                                                                                                                       cTGCTTCGGC‑GGGGTTTTTTTATGGGCACGGGGCGGGGGGAGGTGTCGGGTGCGCTTCGCTTATCCGACCTACAGGGGAGGGTATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCCAAAATAGCTCGGGGGAAAAGGGTTTGAtttttt                >  2:3996559/1‑150 (MQ=255)
                                                                                                                                          cTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAg             <  2:1148345/150‑1 (MQ=255)
                                                                                                                                          cTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAg             <  1:4073573/150‑1 (MQ=255)
                                                                                                                                          cTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAg             <  1:2173274/150‑1 (MQ=255)
                                                                                                                                             cggc‑ggGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGaa          <  2:3399993/150‑1 (MQ=255)
                                                                                                                                             cggc‑ggGGGTTTTTTATTGGCACGGGGCGGGGGGGGTTGTCGGGTGCGCTTCGCTTATCCGACCTACAGGGGAGGGTATTTTAGGCCCGGGGAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGGAACAGGGTTTGATTATTTTAGGaa          >  2:963912/1‑150 (MQ=255)
                                                                                                                                              ggc‑ggGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAg         <  1:1275537/150‑1 (MQ=255)
                                                                                                                                               gc‑ggGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGc        <  2:2271991/150‑1 (MQ=255)
                                                                                                                                                      tttttttATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGAtggt  >  2:4059795/1‑150 (MQ=255)
                                                                                                                                                      tttttttATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGAtggt  <  2:3799474/150‑1 (MQ=255)
                                                                                                                                                      tttttttATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGAtggt  <  2:3750208/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 < 26 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑