Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0123063,412,0490TG26.5% 842.4 / 36.8 352intergenic (+25/‑122)proM/aslBtRNA‑Pro/predicted regulator of arylsulfatase activity
Reads supporting (aligned to +/- strand):  ref base T (97/161);  new base G (93/0);  total (191/161)
Fisher's exact test for biased strand distribution p-value = 1.38e-31
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 381 total aligned reads displayed.

GCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTT  >  AP012306/3411900‑3412198
                                                                                                                                                     |                                                                                                                                                      
     acGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGtttttt                                                                                                                                                   <  2:3067663/150‑1 (MQ=255)
       gACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGGTTTTTGt                                                                                                                                                 >  2:1041203/1‑150 (MQ=255)
        aCTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGtt                                                                                                                                                <  2:3553111/150‑1 (MQ=255)
         cTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGGTTTTTGttt                                                                                                                                               >  2:874369/1‑150 (MQ=255)
          tttAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTTGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGGTTTTTGtttt                                                                                                                                              >  1:4214187/1‑150 (MQ=255)
          tttAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCCCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGtttt                                                                                                                                              <  1:538708/150‑1 (MQ=255)
           ttAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTc                                                                                                                                             <  2:3664684/150‑1 (MQ=255)
           ttAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTc                                                                                                                                             <  1:1823335/150‑1 (MQ=255)
           ttAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGGTTTTTGTTTTc                                                                                                                                             >  2:3991876/1‑150 (MQ=255)
           ttAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGGTTTTTGTTTTc                                                                                                                                             >  1:4078839/1‑150 (MQ=255)
             aaaGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCtg                                                                                                                                           <  1:2793919/150‑1 (MQ=255)
             aaaGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCtg                                                                                                                                           >  2:493071/1‑150 (MQ=255)
             aaaGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCtg                                                                                                                                           >  1:3575585/1‑150 (MQ=255)
                    tGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCAttt                                                                                                                                    >  2:1930077/1‑150 (MQ=255)
                     gAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCGCTCGCCGACCAATTTTGAACCCCGCTTCGGTGGGGGTTTTTATTTTCTTTGCCTTTc                                                                                                                                   >  2:3498316/1‑150 (MQ=39)
                      aaCTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCg                                                                                                                                  <  1:3943697/150‑1 (MQ=255)
                       aCTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGt                                                                                                                                 <  1:1098709/150‑1 (MQ=255)
                        cTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTc                                                                                                                                >  2:590818/1‑150 (MQ=255)
                        cTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGGTTTTTGTTTTCTGTGCCTTTCGTc                                                                                                                                >  1:635623/1‑150 (MQ=255)
                         tAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCa                                                                                                                               <  1:756920/150‑1 (MQ=255)
                         tAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCa                                                                                                                               <  1:164113/150‑1 (MQ=255)
                         tAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCa                                                                                                                               <  1:1324595/150‑1 (MQ=255)
                            aaaTTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGGTTTTTGGTTTCTGTGCATTTTGTCAccc                                                                                                                            >  2:1211999/1‑150 (MQ=255)
                             aaTTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGCAGCGCAACTGGTTTGGGACCAGTGGATCGGAGGTTCGAATCCTCGCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTTTTTTCTTTGGATTTTCTCAcact                                                                                                                           >  2:751107/1‑150 (MQ=255)
                              aTTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCAccctc                                                                                                                          >  1:2605913/1‑150 (MQ=255)
                              aTTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGGTTTTTGTTTTCTGTGCATTTCGTCAccctc                                                                                                                          >  2:1404272/1‑150 (MQ=255)
                               ttCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTTTGTGCATTTCGTCAccctcc                                                                                                                         >  2:3638631/1‑150 (MQ=255)
                                    aaaGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCTGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCg                                                                                                                    >  1:2925723/1‑150 (MQ=255)
                                    aaaGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCg                                                                                                                    <  1:2121416/150‑1 (MQ=255)
                                      aGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCa                                                                                                                  <  2:1469329/150‑1 (MQ=255)
                                        cAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGGTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCCTCGCCAt                                                                                                                >  1:3155834/1‑150 (MQ=255)
                                         aGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGGTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATa                                                                                                               >  2:2803678/1‑150 (MQ=255)
                                           tATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGGTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATaaa                                                                                                             >  2:3913526/1‑150 (MQ=255)
                                               tCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGcc                                                                                                         <  2:1686692/150‑1 (MQ=255)
                                                cGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGccc                                                                                                        >  2:369099/1‑150 (MQ=255)
                                                 ggCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGGTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCg                                                                                                       >  2:2741803/1‑150 (MQ=255)
                                                 ggCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGGTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCg                                                                                                       >  2:2992071/1‑150 (MQ=255)
                                                  gCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGt                                                                                                      <  1:2363387/150‑1 (MQ=255)
                                                   cGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGGTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGta                                                                                                     >  2:2232504/1‑150 (MQ=255)
                                                         gcgcAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACt                                                                                               <  1:1200600/150‑1 (MQ=255)
                                                         gcgcAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACt                                                                                               <  2:4268880/150‑1 (MQ=255)
                                                         gcgcAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACt                                                                                               <  2:329180/150‑1 (MQ=255)
                                                            cAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCAt                                                                                           <  1:1736188/150‑1 (MQ=255)
                                                              gCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTg                                                                                         <  1:306799/150‑1 (MQ=255)
                                                               cTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGc                                                                                        <  2:3336106/150‑1 (MQ=255)
                                                               cTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGc                                                                                        <  1:3090832/150‑1 (MQ=255)
                                                                 tGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGccc                                                                                      >  1:3678454/1‑150 (MQ=255)
                                                                    tAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCAc                                                                                   >  1:3535979/1‑150 (MQ=255)
                                                                     aGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCCTTTCGTCACCCTCCCCTCGCCATAAACGCCCGTAATAACTCATTGGCCCAca                                                                                  >  1:2129698/1‑149 (MQ=255)
                                                                     aGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACg                                                                                  >  1:2462141/1‑150 (MQ=255)
                                                                         cAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTAt                                                                              <  1:1391353/150‑1 (MQ=255)
                                                                           aCTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGa                                                                            <  2:3364070/150‑1 (MQ=255)
                                                                              ggTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGAttt                                                                         >  1:1545771/1‑150 (MQ=255)
                                                                              ggTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGGTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGAttt                                                                         >  2:4154005/1‑150 (MQ=255)
                                                                               gTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTc                                                                        <  1:348076/150‑1 (MQ=255)
                                                                                tttGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCg                                                                       <  1:1019834/150‑1 (MQ=255)
                                                                                 ttGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGc                                                                      <  2:2466015/150‑1 (MQ=255)
                                                                                 ttGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGGTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCCCGGTATGATTTCGc                                                                      >  2:3247718/1‑150 (MQ=255)
                                                                                   gggACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTTATTGCCCCACGGTATGATTTCGccc                                                                    <  1:4166933/150‑1 (MQ=255)
                                                                                    ggACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGGTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTTCCCCACGGGATGATTTCGCCCt                                                                   >  2:930418/1‑150 (MQ=255)
                                                                                     gACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCtt                                                                  >  1:1469611/1‑150 (MQ=255)
                                                                                     gACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCCTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTTGCCCtt                                                                  >  2:4219317/1‑150 (MQ=255)
                                                                                      aCCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGTCCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTTTTTTATGTGCATTTTGTTACCCTCCCCTCGCAATAAACCCCCCTAATAACTCCATTCCCCACCGTATGAATTCGCCCTTa                                                                 >  2:2000951/1‑150 (MQ=255)
                                                                                      aCCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTa                                                                 <  2:1794015/150‑1 (MQ=255)
                                                                                      aCCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGGTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTa                                                                 >  1:1889146/1‑150 (MQ=255)
                                                                                       ccAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTaa                                                                >  1:3389289/1‑150 (MQ=255)
                                                                                        cAGTGGGTCGGGGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGGTTTTTTTTTTTTGTGCATTTTGTCACCCCCCCCTCGCAATAAAAGCCCGTAATAACTCACTTCCCCCCCGGATGGTTTTCCCCCTTaa                                                                >  2:3168258/1‑150 (MQ=255)
                                                                                         aGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACg                                                              <  1:2023093/150‑1 (MQ=255)
                                                                                           tGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTa                                                            >  1:229941/1‑150 (MQ=255)
                                                                                           tGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGGTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTa                                                            >  2:4108235/1‑150 (MQ=255)
                                                                                                 ggAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAg                                                      <  2:2960924/150‑1 (MQ=255)
                                                                                                 ggAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGGTTTTTGTTTTCTGTGCTTTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTTCCCCACGGTATGATTTCGCCCTTAACGTATTGAAg                                                      >  2:2091549/1‑150 (MQ=255)
                                                                                                   aGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGa                                                    >  1:2533755/1‑150 (MQ=255)
                                                                                                   aGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGGTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGa                                                    >  2:2596004/1‑150 (MQ=255)
                                                                                                    ggTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGAt                                                   <  1:2689491/150‑1 (MQ=255)
                                                                                                      ttCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGGTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGa                                                 >  2:1839327/1‑150 (MQ=255)
                                                                                                       tCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGAc                                                >  2:3781484/1‑150 (MQ=255)
                                                                                                        cGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACt                                               >  1:819498/1‑150 (MQ=255)
                                                                                                          aaTCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTTAAGGATGACTTc                                             >  2:111777/1‑150 (MQ=255)
                                                                                                            tCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAg                                           <  2:2670837/150‑1 (MQ=255)
                                                                                                            tCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGGTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAg                                           >  1:3006026/1‑150 (MQ=255)
                                                                                                             ccTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAgg                                          >  1:445342/1‑150 (MQ=255)
                                                                                                                   tcGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGa                                    <  2:942897/150‑1 (MQ=255)
                                                                                                                   tcGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCCCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATTATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGa                                    >  2:845457/1‑150 (MQ=255)
                                                                                                                       cGACCAATTTTGAACCCCGCTTCGGCGGGGGTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTTAAGGATGACTTCAGGCAAGGAGCGa                                >  2:1540481/1‑150 (MQ=255)
                                                                                                                        gACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGAc                               <  2:2372361/150‑1 (MQ=255)
                                                                                                                         aCCAATTTTGAACCCCGCTTCGGCGGGGGTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGGATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGAcc                              >  2:1761162/1‑150 (MQ=255)
                                                                                                                           cAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCAt                            >  1:159038/1‑150 (MQ=255)
                                                                                                                              ttttGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCAtgct                         >  1:1770074/1‑150 (MQ=255)
                                                                                                                                  gAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCaa                     >  2:1300335/1‑150 (MQ=255)
                                                                                                                                     ccccGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCCCGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAg                  >  2:2895138/1‑150 (MQ=255)
                                                                                                                                        cGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGtt               <  2:2047917/150‑1 (MQ=255)
                                                                                                                                        cGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCCCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGtt               >  2:2540623/1‑150 (MQ=255)
                                                                                                                                         gCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTc              <  2:523635/150‑1 (MQ=255)
                                                                                                                                             cggcggGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAAc          <  1:4001494/150‑1 (MQ=255)
                                                                                                                                             cggcggGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAAc          <  1:1694179/150‑1 (MQ=255)
                                                                                                                                               gcggGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGGAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAAcgc        <  1:3819967/150‑1 (MQ=255)
                                                                                                                                                cggGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAAcgcg       >  1:923407/1‑150 (MQ=255)
                                                                                                                                                 ggggTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGt      >  1:333417/1‑150 (MQ=255)
                                                                                                                                                 ggggTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGt      >  2:3016314/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTT  >  AP012306/3411900‑3412198

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: