Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0123061,074,6190TC21.7% 380.0 / 99.9 198S29G (AGC→GGC) ldrAtoxic polypeptide, small
Reads supporting (aligned to +/- strand):  ref base T (108/47);  new base C (14/29);  total (122/76)
Fisher's exact test for biased strand distribution p-value = 1.54e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.70e-03
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 221 total aligned reads displayed.

GTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCAAGATTAG  >  AP012306/1074470‑1074768
                                                                                                                                                     |                                                                                                                                                     
gTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGCTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCt                                                                                                                                                       >  1:1114087/1‑150 (MQ=255)
 tATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTg                                                                                                                                                      <  1:2271665/150‑1 (MQ=255)
 tATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCACCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTg                                                                                                                                                      >  2:1849856/1‑150 (MQ=255)
  aTTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGa                                                                                                                                                     >  2:3525589/1‑150 (MQ=255)
  aTTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGa                                                                                                                                                     <  2:383761/150‑1 (MQ=255)
   ttAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGAc                                                                                                                                                    <  2:3190390/150‑1 (MQ=255)
   ttAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGAc                                                                                                                                                    >  2:664552/1‑150 (MQ=255)
   ttAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGAc                                                                                                                                                    >  2:1220315/1‑150 (MQ=255)
    tAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACa                                                                                                                                                   >  1:1810805/1‑150 (MQ=255)
     aGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACaa                                                                                                                                                  >  1:2088551/1‑150 (MQ=255)
     aGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACaa                                                                                                                                                  >  1:1617888/1‑150 (MQ=255)
       cGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATc                                                                                                                                                >  1:3025323/1‑150 (MQ=255)
       cGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATc                                                                                                                                                >  1:4040045/1‑150 (MQ=255)
       cGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATc                                                                                                                                                >  1:1365566/1‑150 (MQ=255)
       cGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATc                                                                                                                                                <  2:3789544/150‑1 (MQ=255)
        gTTGTGCCTTAAGACTTTTTGTATGGGTGCCGAATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACTAGCCGACAATCg                                                                                                                                               <  2:4396709/150‑1 (MQ=14)
        gTTCTGCGTTAAGTCTTTTTTCATGGGTGCCAGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCg                                                                                                                                               <  1:3776644/150‑1 (MQ=255)
        gTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCg                                                                                                                                               <  2:1727132/150‑1 (MQ=255)
         ttCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGc                                                                                                                                              >  1:3268348/1‑150 (MQ=255)
          tCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCt                                                                                                                                             <  2:2680184/150‑1 (MQ=255)
          tCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCt                                                                                                                                             >  1:3136102/1‑150 (MQ=255)
          tCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTACGCCACCAGCTGACAATCGCTg                                                                                                                                            >  2:794845/1‑150 (MQ=255)
           cTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTg                                                                                                                                            >  2:896568/1‑150 (MQ=255)
           cTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTg                                                                                                                                            >  1:2117954/1‑150 (MQ=255)
            tGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGc                                                                                                                                           >  1:4349608/1‑150 (MQ=255)
              cGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCgg                                                                                                                                         >  1:2813229/1‑150 (MQ=255)
              cGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCgg                                                                                                                                         >  1:2015571/1‑150 (MQ=255)
              cGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCgg                                                                                                                                         >  1:3585243/1‑150 (MQ=255)
               gTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGt                                                                                                                                        >  2:1904881/1‑150 (MQ=255)
                ttAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGta                                                                                                                                       >  2:1679207/1‑150 (MQ=255)
                ttAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGta                                                                                                                                       <  1:179562/150‑1 (MQ=255)
                ttAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGta                                                                                                                                       >  1:3913080/1‑150 (MQ=255)
                ttAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGta                                                                                                                                       >  1:1442755/1‑150 (MQ=255)
                 tAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGtaa                                                                                                                                      >  2:2782491/1‑150 (MQ=255)
                  aaGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCCCCAGCTGACAATCGCTGCGGtaat                                                                                                                                     >  2:733827/1‑150 (MQ=255)
                  aaGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGtaat                                                                                                                                     >  1:1343510/1‑150 (MQ=255)
                  aaGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGtaat                                                                                                                                     <  2:4236411/150‑1 (MQ=255)
                    gACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGtaataa                                                                                                                                   <  2:4159765/150‑1 (MQ=255)
                    gACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGtaataa                                                                                                                                   >  1:457039/1‑150 (MQ=255)
                    gACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGtaataa                                                                                                                                   <  2:2446129/150‑1 (MQ=255)
                    gACGTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCGGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACTTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGtaataa                                                                                                                                   <  2:337716/150‑1 (MQ=255)
                     aCTTTTTTCATGGGTGCCAGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTaataat                                                                                                                                  <  2:2360169/150‑1 (MQ=255)
                       ttttttCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTc                                                                                                                                <  2:2900899/150‑1 (MQ=255)
                       ttttttCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTc                                                                                                                                >  1:169046/1‑150 (MQ=255)
                       ttttttCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTc                                                                                                                                <  2:1480610/150‑1 (MQ=255)
                        tttttCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTcc                                                                                                                               >  1:112919/1‑150 (MQ=255)
                        tttttCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTcc                                                                                                                               >  2:4378736/1‑150 (MQ=255)
                          tttCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCg                                                                                                                             <  2:749999/150‑1 (MQ=255)
                          tttCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCg                                                                                                                             <  2:1818253/150‑1 (MQ=255)
                           ttCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGc                                                                                                                            >  2:3973754/1‑150 (MQ=255)
                           ttCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGc                                                                                                                            >  2:386553/1‑150 (MQ=255)
                           ttCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGc                                                                                                                            >  2:4055496/1‑150 (MQ=255)
                           ttCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGc                                                                                                                            >  1:3348768/1‑150 (MQ=255)
                            tCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGcc                                                                                                                           >  1:623092/1‑150 (MQ=255)
                            tCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGcc                                                                                                                           >  2:4254375/1‑150 (MQ=255)
                            tCATGGGTGCCAAATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGcc                                                                                                                           <  1:2757674/150‑1 (MQ=255)
                             cATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCa                                                                                                                          >  1:981048/1‑150 (MQ=255)
                             cATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCa                                                                                                                          >  2:620550/1‑150 (MQ=255)
                              aTGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAg                                                                                                                         >  1:114507/1‑150 (MQ=255)
                               tGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTCCCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTCCTCCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAgg                                                                                                                        <  2:2801480/150‑1 (MQ=255)
                               tGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAgg                                                                                                                        >  2:2015083/1‑150 (MQ=255)
                                 ggTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGAt                                                                                                                      >  2:1280381/1‑150 (MQ=40)
                                 ggTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGAt                                                                                                                      >  1:3047794/1‑150 (MQ=38)
                                  gTGCCGGATACAAAAAAGGCCGCAGGCTGTTGCCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATc                                                                                                                     <  1:802321/150‑1 (MQ=255)
                                  gTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATc                                                                                                                     <  1:2911399/150‑1 (MQ=39)
                                  gTGCCGGATAAAAAAAAGGCCGCAGGCTGTTCCCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATc                                                                                                                     <  2:1570814/150‑1 (MQ=255)
                                     ccGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGt                                                                                                                  >  1:138580/1‑150 (MQ=38)
                                      cGGATACAAAAAAGCCCGCAGGCTGTTCCCCCTGCGGCGGGTTTCGGGCGAATATTGCCATCACGGCAGCCTGCCGCCCGTTTTCCCCTTCCTCCCGGTTCCCCCACCAGCTGAAAATCGCGGGGGTAATAATTCCCGCCAGGATCGGtg                                                                                                                 <  2:3210905/150‑1 (MQ=11)
                                       ggATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGtgc                                                                                                                >  2:1804168/1‑150 (MQ=35)
                                       ggATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGtgc                                                                                                                >  2:3254715/1‑150 (MQ=35)
                                       ggATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGtgc                                                                                                                >  2:3022909/1‑150 (MQ=17)
                                        gATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCCCCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGtgcc                                                                                                               <  2:3905861/150‑2 (MQ=14)
                                        gATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGtgcc                                                                                                               >  2:46441/1‑149 (MQ=14)
                                         aTACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGtgctg                                                                                                              >  1:870607/1‑150 (MQ=33)
                                         aTACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGtgctg                                                                                                              >  1:3490673/1‑150 (MQ=33)
                                             aaaaaaGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAg                                                                                                          >  2:2743709/1‑150 (MQ=32)
                                               aaaaGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGt                                                                                                        <  1:1972169/150‑1 (MQ=32)
                                               aaaaGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGt                                                                                                        >  1:656391/1‑150 (MQ=17)
                                                aaaGGCCGCAGGCTGTTCCCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTc                                                                                                       <  2:3881737/150‑1 (MQ=17)
                                                                                                                                  ttCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAt                     <  1:565072/150‑1 (MQ=11)
                                                                                                                                  ttCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAt                     <  1:2543546/150‑1 (MQ=11)
                                                                                                                                  ttCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAt                     >  1:892132/1‑150 (MQ=11)
                                                                                                                                    ccGGTTACCCCACCGTCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCAGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTg                   <  2:3637701/150‑1 (MQ=255)
                                                                                                                                     cGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGt                  >  2:1776661/1‑150 (MQ=32)
                                                                                                                                      ggTTACGCCCCCAGCTGACAATCGCTGGGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGtt                 <  2:1395293/150‑1 (MQ=17)
                                                                                                                                      ggTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGtt                 >  1:1038750/1‑150 (MQ=32)
                                                                                                                                         tACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTg              <  2:2816136/150‑1 (MQ=35)
                                                                                                                                           cGCCACCAGCTGACAATGGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGt            <  2:1099535/150‑1 (MQ=37)
                                                                                                                                           cGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGt            >  2:1889529/1‑150 (MQ=37)
                                                                                                                                             ccaccaGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTc          >  1:70622/1‑150 (MQ=38)
                                                                                                                                             cacccaGATGAAAATAGCGGGGGAAAAATTTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGTGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTc          <  2:2621699/147‑1 (MQ=11)
                                                                                                                                              caccaGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCa         >  2:2231853/1‑150 (MQ=38)
                                                                                                                                              caccaGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCa         <  1:4061134/150‑1 (MQ=38)
                                                                                                                                               accaGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCaa        >  1:3221456/1‑150 (MQ=39)
                                                                                                                                                 caGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCAAGa      <  2:2999884/150‑1 (MQ=255)
                                                                                                                                                 caGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCAAGa      <  2:1076260/150‑1 (MQ=255)
                                                                                                                                                  aGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCAAGAt     <  1:1330302/150‑1 (MQ=255)
                                                                                                                                                   gCTGACAATCGCTGCGGTAATAATCCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCAAGAtt    <  2:2273213/150‑1 (MQ=255)
                                                                                                                                                    cTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCAAGATTa   >  1:1624384/1‑150 (MQ=40)
                                                                                                                                                     tGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCAAGATTAg  <  1:107913/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCAAGATTAG  >  AP012306/1074470‑1074768

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