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

TGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGT  >  AP012306/1091845‑1092141
                                                                                                                                                     |                                                                                                                                                     
tGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGgtt                                                                                                                                                     >  1:3450773/1‑150 (MQ=255)
  gggATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGttttt                                                                                                                                                   <  1:1917003/150‑1 (MQ=255)
  gggATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGttttt                                                                                                                                                   <  1:1453958/150‑1 (MQ=255)
  gggATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGttttt                                                                                                                                                   <  1:1682650/150‑1 (MQ=255)
  gggATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGgtttt                                                                                                                                                   >  1:2285261/1‑150 (MQ=255)
    gATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGttttttt                                                                                                                                                 <  2:2399042/150‑1 (MQ=255)
    gATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGctttttt                                                                                                                                                 >  1:1802232/1‑150 (MQ=255)
       gACGTTATTCCTGCTGTTCCCGTTCTCGCGTATGATACACAACTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTTTg                                                                                                                                              >  2:3633307/1‑150 (MQ=255)
        aCGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGAGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGGTTTTTTATgg                                                                                                                                             >  2:2772127/1‑150 (MQ=255)
        aCGTTATTCCTGCTGTTCCCGTTCCCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATgg                                                                                                                                             <  2:3397959/150‑1 (MQ=255)
         cGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGGTTTTTTTTggg                                                                                                                                            >  2:2918320/1‑150 (MQ=255)
          gTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGGTTTTTTTTGGGc                                                                                                                                           >  2:772739/1‑150 (MQ=255)
          gTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGAGGTTTTTTAttagc                                                                                                                                           >  1:877317/1‑146 (MQ=255)
          gTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGG‑GCGCGCCCGTCCCTAAGCGAATTTTAGTTCACCTAGACCCCGCCTTCGGC‑GGGGTTTTTTTTTGGGc                                                                                                                                           >  2:2103255/1‑150 (MQ=255)
           ttATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGGTTTTTTTTTGGCa                                                                                                                                          >  2:1516135/1‑150 (MQ=255)
               tCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGGTTTTTTATGGGCGCGGt                                                                                                                                      >  2:2279613/1‑150 (MQ=255)
                    ctgTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGGTTTTTTATGGGCACGGGGCggg                                                                                                                                 >  1:1723066/1‑150 (MQ=255)
                     tgTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCgggg                                                                                                                                <  1:1017770/150‑1 (MQ=255)
                       ttCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTg                                                                                                                              <  2:424585/150‑1 (MQ=255)
                       ttCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTg                                                                                                                              <  1:3888570/150‑1 (MQ=255)
                        tCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGGTTTTTTTTGGGCACGGTGCGGGGGGa                                                                                                                             >  2:1917127/1‑150 (MQ=255)
                          ccGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGt                                                                                                                           <  1:4071830/150‑1 (MQ=255)
                          ccGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGt                                                                                                                           <  1:3564010/150‑1 (MQ=255)
                            gTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGCGGGGGTTTTTTTTTGGGC‑CGGGGCGGGGGGAGTTg                                                                                                                         >  2:4144465/1‑150 (MQ=255)
                            gTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATTGGCACGGTGCGGGGGGAGTTg                                                                                                                         >  2:507449/1‑150 (MQ=255)
                            gTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGTACGGTGCGGGGTGAGGTg                                                                                                                         >  1:3196555/1‑150 (MQ=255)
                            gTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGGTg                                                                                                                         >  1:2597524/1‑150 (MQ=255)
                             ttCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGGTTTTTTTTTGGGACGGGGCGGGGGGAGgtgt                                                                                                                        >  2:899209/1‑150 (MQ=255)
                              tctcGCGTCTGATACACATCTGGAGCGTACTGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTc                                                                                                                       <  1:1112119/150‑1 (MQ=255)
                              tctcGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGGTTTTTTTTAGGCACCGTGCGGGGTGGGtttgc                                                                                                                       >  1:2391808/1‑147 (MQ=255)
                                    gTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATgc                                                                                                                 >  1:1280728/1‑150 (MQ=255)
                                     tCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGGTTTTTTTTGGGCACGGTGCGGGGGGAGGTGTCGGATggg                                                                                                                >  2:1794033/1‑148 (MQ=255)
                                      cTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATgcgc                                                                                                               >  2:688849/1‑150 (MQ=255)
                                       tGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGcgct                                                                                                              <  2:1866092/150‑1 (MQ=255)
                                       tGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGGTGcgct                                                                                                              >  1:476870/1‑150 (MQ=255)
                                          tACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGcgcttcg                                                                                                           <  1:1478482/150‑1 (MQ=255)
                                          tACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGGTTTTTTTTTGGCACGGTGCGGGGTGAGTTGTCGGATGcgcttcg                                                                                                           >  2:2220107/1‑150 (MQ=255)
                                           acacATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGcgcttcgc                                                                                                          >  1:3424320/1‑150 (MQ=255)
                                            cacaTCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGcgcttcgct                                                                                                         <  1:375243/150‑1 (MQ=255)
                                                      gCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCt                                                                                               <  1:3468427/150‑1 (MQ=255)
                                                      gCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCt                                                                                               <  2:4039974/150‑1 (MQ=255)
                                                       cGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTa                                                                                              <  1:958736/150‑1 (MQ=255)
                                                       cGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTa                                                                                              <  1:4010402/150‑1 (MQ=255)
                                                        gTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGGTTTTTTTTGGGCACGGTGCGGGGGGGGTTGTCGGATGCGCTTCGCTTTTCCGACCCAc                                                                                             >  2:3333999/1‑150 (MQ=255)
                                                          aCCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGGTTTTTTTTGGGCACGGTGCGGGGGGAGTTGTCGGATGGGCTTCGCTTATCCGACCCACAg                                                                                           >  2:3294388/1‑150 (MQ=255)
                                                           ccGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAgg                                                                                          <  1:856517/150‑1 (MQ=255)
                                                            cGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAggg                                                                                         <  1:2027652/150‑1 (MQ=255)
                                                            cGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAggg                                                                                         >  1:249529/1‑150 (MQ=255)
                                                            cGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAggg                                                                                         >  1:207760/1‑150 (MQ=255)
                                                              gTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGgga                                                                                       <  2:4067294/150‑1 (MQ=255)
                                                               tGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGGTGCGCTTCGCTTATCCGACCTACAGGggag                                                                                      >  2:3672853/1‑150 (MQ=255)
                                                                  aGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGat                                                                                   <  1:2523598/150‑1 (MQ=255)
                                                                  aGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGat                                                                                   <  2:1732995/150‑1 (MQ=255)
                                                                  aGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGGTTTTTTATGGGCACGGGGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGat                                                                                   >  2:260146/1‑150 (MQ=255)
                                                                      tCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGGTATTg                                                                               >  2:3887667/1‑150 (MQ=255)
                                                                      tCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGGTTTTTTTTGGGGACGGGGCGGGGGGGGGTGTCGGGTGGGCTTCGCTTTTCCGACCTACAGGGGGGGGGATTg                                                                               >  2:2016668/1‑150 (MQ=255)
                                                                      tCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGGTATTg                                                                               >  2:2672474/1‑150 (MQ=255)
                                                                      tCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGGTATTg                                                                               >  1:4035103/1‑150 (MQ=255)
                                                                         gACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAg                                                                            <  2:1574366/150‑1 (MQ=255)
                                                                          acacGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAgg                                                                           <  1:612301/150‑1 (MQ=255)
                                                                              gTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGGTATTGTAGGCCCg                                                                       >  1:2220882/1‑150 (MQ=255)
                                                                               tAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCgg                                                                      <  1:2641036/150‑1 (MQ=255)
                                                                                  gTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGGTTTTTTATGGGCACGGGGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGGGATTTTAGGGCCGGTaa                                                                   >  2:1673354/1‑150 (MQ=255)
                                                                                   tACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGGGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGATATTGTAGGCCCGGTaaa                                                                  >  1:3836984/1‑150 (MQ=255)
                                                                                      cAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATGGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACgt                                                               <  1:41343/150‑1 (MQ=255)
                                                                                      cAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGGTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGGTATTGTAGGCCCGGGAAACgt                                                               >  2:3028984/1‑150 (MQ=255)
                                                                                             tGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCa                                                        <  1:3137204/150‑1 (MQ=255)
                                                                                                  gcTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGATATTGTAGGCCCGGTAAACGTGGTGCCATCggg                                                   >  1:1483117/1‑150 (MQ=255)
                                                                                                      gTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCaaa                                               <  1:2162473/150‑1 (MQ=255)
                                                                                                         aCTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGGTATTGTAGGCCCGGTAAACGTGGGGCCATCGGGGAAAAtt                                            >  2:3254855/1‑150 (MQ=255)
                                                                                                          cTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTa                                           >  1:2938088/1‑150 (MQ=255)
                                                                                                            aaGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGGTATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGc                                         >  2:3794445/1‑150 (MQ=255)
                                                                                                             aGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCt                                        <  1:3906890/150‑1 (MQ=255)
                                                                                                              gCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTc                                       >  1:3004923/1‑150 (MQ=255)
                                                                                                                gAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCgg                                     >  1:3317178/1‑150 (MQ=255)
                                                                                                                 aaTTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGa                                    <  1:1950188/150‑1 (MQ=255)
                                                                                                                  aTTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGGTATTGTAGGCCCCGGAAACGTGGTGCCATCGGGCAAAATTAGCTCgggg                                   >  2:737479/1‑148 (MQ=255)
                                                                                                                   ttttAGTTCTCATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGt                                  <  1:2114889/150‑1 (MQ=255)
                                                                                                                      tAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAAc                               <  1:1569576/150‑1 (MQ=255)
                                                                                                                       aGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACa                              <  2:2829265/150‑1 (MQ=255)
                                                                                                                         ttCACATAGACCCTG‑CTTCGGC‑GGGGGTTTTTTATGGGCACGGGGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGGTATTGTAGGCCCGGTAAACGTGGTGCCCTCGGGCAAAATTAGCTCGGGGGAACAgg                            >  2:441091/1‑150 (MQ=255)
                                                                                                                          tCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGt                           <  1:2691182/150‑1 (MQ=255)
                                                                                                                           cacaTAGACCCTG‑CTTCGGC‑GGGGGTTTTTTATGGGCACGGGGCGGGGGGAGTTGTCGGGTGCGCTTCGCTTATCCGACCTACAGGGGGGGGTATTGTAGGCCCGGGAAACGTGGGGCCATCGGGCAAAATTAGCTCGGAGGAACAgggt                          >  2:1413111/1‑148 (MQ=255)
                                                                                                                           cacaTAGACCCTG‑CTTCGGC‑GGGGGTTTTTTATGGGCACGGGGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGGTATTGTAGGCCCGGTAAACGTGGTGCCATCGGGGAAAATTAGCTCCGAGGAAAAgggt                          >  2:3139056/1‑148 (MQ=255)
                                                                                                                             caTAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGGACGGTGCGGGGGGAGTTGTCGGGTGCGCTTCGCTTATTCGACCTAAAGGGGAGGATATTGTAGGCCCGGGAAACGTGGTGCCATCGGGCCAAATTAGCTCGGGGAAACAGGgttt                        >  2:1475512/1‑150 (MQ=255)
                                                                                                                             caTAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGtttt                        <  1:1182684/150‑1 (MQ=255)
                                                                                                                             caTAGACCCTG‑CTTCGGC‑GGGGGTTTTTTATGGGCACGGGGCGGGGGGGGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGGGGATATTGTAGGCCCGGGAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGgttt                        >  2:3481996/1‑150 (MQ=255)
                                                                                                                             caTAGACCCTG‑CTTCGGC‑GGGGGTTTTTTATGGGCACGGGGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGGTATTGTAGGCCCGGTAAACGTGGGGCCATCGGGCAAAATTAGCTCGGAGGAACAGGgttt                        >  2:853745/1‑150 (MQ=255)
                                                                                                                                aGACCCTG‑CTTCGGC‑GGGGGTTTTTTATGGGCACGGGGCGGGGGGAGTTGTCGGGTGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGGAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAAAAGGGTTTGat                     >  2:3623077/1‑150 (MQ=255)
                                                                                                                                 gACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGatt                    <  1:2755708/150‑1 (MQ=255)
                                                                                                                                   cccTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGattat                  <  1:1316614/150‑1 (MQ=255)
                                                                                                                                   cccTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGattat                  <  2:2330916/150‑1 (MQ=255)
                                                                                                                                      tG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTa               <  2:124173/150‑1 (MQ=255)
                                                                                                                                      tG‑CTTCGGC‑GGGGGTTTTTTATGGGCACGGGGCGGGGGGAGTTGTCGGGTGGGCTTCGCTTATCCGACCTACAGGGGGGGGTATTTTAGGCCCGGTAAACGTGGGGCCATCGGGGAAAATTAGCTCGGGGGAACAGGGTTTGATTAtttt               >  2:1057663/1‑149 (MQ=255)
                                                                                                                                          ttCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAgg            <  1:416040/150‑1 (MQ=255)
                                                                                                                                             ggc‑ggGTTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAg         <  2:2880047/150‑1 (MQ=255)
                                                                                                                                              gc‑ggGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGc        <  1:2126908/150‑1 (MQ=255)
                                                                                                                                                 ggggTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGa      <  1:78336/150‑1 (MQ=255)
                                                                                                                                                   ggTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGAtg    >  2:399179/1‑150 (MQ=255)
                                                                                                                                                     tttttttATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGAtggt  <  1:1015745/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTG‑CTTCGGC‑GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGT  >  AP012306/1091845‑1092141

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