Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0123063,412,0490TG23.6% 1641.8 / 68.9 646intergenic (+25/‑122)proM/aslBtRNA‑Pro/predicted regulator of arylsulfatase activity
Reads supporting (aligned to +/- strand):  ref base T (176/315);  new base G (152/0);  total (331/315)
Fisher's exact test for biased strand distribution p-value = 0.00e+00
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 700 total aligned reads displayed.

GCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTT  >  AP012306/3411900‑3412198
                                                                                                                                                      |                                                                                                                                                     
  accacGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑Gttt                                                                                                                                                     >  1:413867/1‑150 (MQ=255)
  accacGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑Gttt                                                                                                                                                     <  2:145812/150‑1 (MQ=255)
   ccacGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑Gattt                                                                                                                                                    >  1:3142435/1‑150 (MQ=255)
     acGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑Gtttttt                                                                                                                                                  >  1:2349722/1‑150 (MQ=255)
     acGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑Ggttttt                                                                                                                                                  >  1:901109/1‑150 (MQ=255)
       gACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGt                                                                                                                                                >  1:3889966/1‑150 (MQ=255)
       gACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGt                                                                                                                                                <  2:3695953/150‑1 (MQ=255)
        aCTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGtt                                                                                                                                               <  1:212771/150‑1 (MQ=255)
         cTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGttt                                                                                                                                              <  1:4396319/150‑1 (MQ=255)
          tttAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGtttt                                                                                                                                             <  2:3156387/150‑1 (MQ=255)
           ttAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GGTTTTTGTTTTc                                                                                                                                            >  1:1079553/1‑150 (MQ=255)
             aaaGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GGTTTTTGTTTTCtg                                                                                                                                          >  1:346950/1‑150 (MQ=255)
              aaGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GGTTTTTGTTTTCtgt                                                                                                                                         >  2:3100658/1‑150 (MQ=255)
              aaGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTCGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGAGTttttttttttttt                                                                                                                                          >  2:3609159/1‑147 (MQ=39)
               aGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GGTTTTTGTTTTCtgtg                                                                                                                                        >  1:3711238/1‑150 (MQ=255)
                 aaTTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GGTTTTTTTTTTCTTTGCa                                                                                                                                      >  2:323819/1‑150 (MQ=255)
                  aTTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCAt                                                                                                                                     <  2:840244/150‑1 (MQ=255)
                    tGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GGTTTTTGTTTTCTGTGGAttt                                                                                                                                   >  2:1252284/1‑150 (MQ=255)
                    tGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GGTTTTTGTTTTCTGTGCAttt                                                                                                                                   >  2:139058/1‑150 (MQ=255)
                     gAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GGTTTTTGTTTTCTGTGCATTTc                                                                                                                                  >  1:1255784/1‑150 (MQ=255)
                      aaCTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GGTTTTTGTTTTCTGTGCATTTCg                                                                                                                                 >  2:2434341/1‑150 (MQ=255)
                        cTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCGATTTTGAACCCCGCTTCGGCGGG‑GGTTTTTGTTTTCTGTGCATTTCGTc                                                                                                                               >  2:209944/1‑150 (MQ=255)
                         tAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GGTTTTTGTTTTCTGTGCATTTCGTCa                                                                                                                              >  1:2556311/1‑150 (MQ=255)
                          aaaaaTTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GGTTTTTGTTTTCTGTGCATTTCGCCAc                                                                                                                             >  1:2013873/1‑150 (MQ=255)
                             aaTTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCAccct                                                                                                                          >  2:1320458/1‑150 (MQ=255)
                             aaTTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCAccct                                                                                                                          <  2:151466/150‑1 (MQ=255)
                              aTTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GGTTTTTGTTTTCTGTGCATTTCGTCAccctc                                                                                                                         >  2:1122837/1‑150 (MQ=255)
                               ttCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTTTGTGCATTTCGTCAccctcc                                                                                                                        >  2:3537060/1‑150 (MQ=255)
                                tCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCAccctccc                                                                                                                       >  1:2702591/1‑150 (MQ=255)
                                tCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCAccctccc                                                                                                                       <  1:1658026/150‑1 (MQ=255)
                                 cAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCAccctccct                                                                                                                      >  1:497765/1‑150 (MQ=255)
                                  aaaaaGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCtt                                                                                                                     <  2:917192/150‑1 (MQ=255)
                                  aaaaaGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GGTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCtt                                                                                                                     >  2:448206/1‑150 (MQ=255)
                                   aaaaGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTc                                                                                                                    <  2:3464808/150‑1 (MQ=255)
                                    aaaGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCg                                                                                                                   <  1:150418/150‑1 (MQ=255)
                                    aaaGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCg                                                                                                                   >  1:3032717/1‑150 (MQ=255)
                                       gCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GGTTTTTGGTTTCTGTGCATTTCGTCACCCTCCCTTCGCaa                                                                                                                >  1:2944447/1‑150 (MQ=255)
                                        cAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAAt                                                                                                               >  1:3633031/1‑150 (MQ=255)
                                         aGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATa                                                                                                              <  2:1980652/150‑1 (MQ=255)
                                          gTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATaa                                                                                                             >  2:735466/1‑150 (MQ=255)
                                          gTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATaa                                                                                                             <  2:255290/150‑1 (MQ=255)
                                          gTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GGTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATaa                                                                                                             >  2:3955835/1‑150 (MQ=255)
                                           tATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTTTGCATTTCGTCACCCTCCCTTCGCAATtaa                                                                                                            >  1:359299/1‑150 (MQ=255)
                                           tATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATaaa                                                                                                            >  1:2514958/1‑150 (MQ=255)
                                           tATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATaaa                                                                                                            >  1:275221/1‑150 (MQ=255)
                                           tATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATaaa                                                                                                            <  2:3344386/150‑1 (MQ=255)
                                                 ggCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GGTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAAAGCCCg                                                                                                      >  2:1036253/1‑150 (MQ=255)
                                                  gCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGt                                                                                                     <  1:427254/150‑1 (MQ=255)
                                                  gCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGt                                                                                                     <  2:289089/150‑1 (MQ=255)
                                                   cGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GGTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGta                                                                                                    >  2:218089/1‑150 (MQ=255)
                                                    gAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGtaa                                                                                                   <  1:2263433/150‑1 (MQ=255)
                                                       tAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTTGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGtaataa                                                                                                >  1:4231447/1‑150 (MQ=255)
                                                         gcgcAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACt                                                                                              <  2:3765517/150‑1 (MQ=255)
                                                           gcAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GGTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAAATCa                                                                                            >  2:2339108/1‑150 (MQ=255)
                                                            cAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCAt                                                                                           >  1:3798172/1‑150 (MQ=255)
                                                            cAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCAt                                                                                           >  1:1577176/1‑150 (MQ=255)
                                                             aGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCAtt                                                                                          >  1:1922170/1‑150 (MQ=255)
                                                               cTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGc                                                                                        <  1:2783944/150‑1 (MQ=255)
                                                                  ggTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGcccc                                                                                     <  2:3002663/150‑1 (MQ=255)
                                                                  ggTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGcccc                                                                                     <  1:759275/150‑1 (MQ=255)
                                                                   gTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GGTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTTCCCCa                                                                                    >  2:620900/1‑150 (MQ=255)
                                                                    tAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCAc                                                                                   >  1:1216819/1‑150 (MQ=255)
                                                                       cgcAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGt                                                                                <  1:1831665/150‑1 (MQ=255)
                                                                           aCTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGa                                                                            <  1:636045/150‑1 (MQ=255)
                                                                           aCTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGa                                                                            <  2:225915/150‑1 (MQ=255)
                                                                           aCTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GGTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGa                                                                            >  2:3901889/1‑150 (MQ=255)
                                                                                tttGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAAATCATTGCCCCACGGTATGATTTCg                                                                       >  1:3252734/1‑150 (MQ=255)
                                                                                  tGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGcc                                                                     >  1:4364831/1‑150 (MQ=255)
                                                                                    ggACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCt                                                                   <  1:3334678/150‑1 (MQ=255)
                                                                                    ggACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCt                                                                   <  1:3755088/150‑1 (MQ=255)
                                                                                     gACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCtt                                                                  >  1:1965814/1‑150 (MQ=255)
                                                                                     gACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCtt                                                                  <  2:1175293/150‑1 (MQ=255)
                                                                                       ccAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTaa                                                                <  2:1819607/150‑1 (MQ=255)
                                                                                        cAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAAc                                                               <  1:349571/150‑1 (MQ=255)
                                                                                           tGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTa                                                            <  2:3284577/150‑1 (MQ=255)
                                                                                             ggTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTAtt                                                          <  2:3853991/150‑1 (MQ=255)
                                                                                             ggTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GGTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTAtt                                                          >  2:189862/1‑150 (MQ=255)
                                                                                                  gAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAgg                                                     >  1:599488/1‑150 (MQ=255)
                                                                                                    ggTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGAt                                                   >  1:2192234/1‑150 (MQ=255)
                                                                                                    ggTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGAt                                                   >  1:1413146/1‑150 (MQ=255)
                                                                                                    ggTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGAt                                                   <  1:286892/150‑1 (MQ=255)
                                                                                                     gTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCACTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATg                                                  >  2:4133643/1‑150 (MQ=255)
                                                                                                      ttCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGa                                                 <  1:710109/150‑1 (MQ=255)
                                                                                                      ttCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGATTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGa                                                 >  2:2658529/1‑150 (MQ=255)
                                                                                                       tCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGAc                                                <  1:826459/150‑1 (MQ=255)
                                                                                                        cGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACt                                               >  1:3071974/1‑150 (MQ=255)
                                                                                                         gAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACtt                                              <  1:100664/150‑1 (MQ=255)
                                                                                                          aaTCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTc                                             >  1:1731275/1‑150 (MQ=255)
                                                                                                          aaTCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTc                                             <  2:522619/150‑1 (MQ=255)
                                                                                                           aTCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCa                                            <  2:4252313/150‑1 (MQ=255)
                                                                                                            tCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTCTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAg                                           <  2:4031098/150‑1 (MQ=255)
                                                                                                                 tctcGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAg                                      >  1:1314259/1‑150 (MQ=255)
                                                                                                                     gCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGc                                  >  2:4176283/1‑150 (MQ=255)
                                                                                                                           cAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCAt                            <  2:1682926/150‑1 (MQ=255)
                                                                                                                             aTTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCAtgc                          <  2:272164/150‑1 (MQ=255)
                                                                                                                              ttttGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCAtgct                         <  1:96369/150‑1 (MQ=255)
                                                                                                                                         gCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTc              <  1:4058837/150‑1 (MQ=255)
                                                                                                                                                cggG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAAcgcg       >  2:1610364/1‑150 (MQ=255)
                                                                                                                                                cggG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAAcgcg       >  1:3947354/1‑150 (MQ=255)
                                                                                                                                                   g‑gTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGc    >  1:3219320/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGG‑GTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTT  >  AP012306/3411900‑3412198

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