Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0123061,074,6190TC21.6% 492.6 / 121.2 255S29G (AGC→GGC) ldrAtoxic polypeptide, small
Reads supporting (aligned to +/- strand):  ref base T (141/59);  new base C (19/36);  total (160/95)
Fisher's exact test for biased strand distribution p-value = 2.88e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.49e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 281 total aligned reads displayed.

GTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCAAGATTAG  >  AP012306/1074470‑1074768
                                                                                                                                                     |                                                                                                                                                     
gTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCt                                                                                                                                                       >  1:2286299/1‑150 (MQ=255)
gTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCt                                                                                                                                                       <  1:9835/150‑1 (MQ=255)
 tATTGGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTg                                                                                                                                                      >  2:745867/1‑150 (MQ=255)
  aTTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGa                                                                                                                                                     >  2:1742455/1‑150 (MQ=255)
  aTTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCAGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGa                                                                                                                                                     >  1:3691480/1‑150 (MQ=255)
    tAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACa                                                                                                                                                   <  1:954819/150‑1 (MQ=255)
    tAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACa                                                                                                                                                   <  1:4653542/150‑1 (MQ=255)
    tAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACa                                                                                                                                                   >  1:3895877/1‑150 (MQ=255)
     aGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACaa                                                                                                                                                  >  2:2654216/1‑150 (MQ=255)
     aGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACaa                                                                                                                                                  >  2:350590/1‑150 (MQ=255)
      gCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAAt                                                                                                                                                 >  1:1594619/1‑150 (MQ=255)
      gCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAAt                                                                                                                                                 >  1:3729422/1‑150 (MQ=255)
      gCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAAt                                                                                                                                                 <  2:3559575/150‑1 (MQ=255)
       cGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATc                                                                                                                                                >  2:3596945/1‑150 (MQ=255)
       cGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATc                                                                                                                                                <  2:612123/150‑1 (MQ=255)
        gTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCTGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCg                                                                                                                                               >  2:4337155/1‑150 (MQ=255)
        gTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAAttg                                                                                                                                               >  1:1655338/1‑148 (MQ=255)
          tCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCt                                                                                                                                             <  2:1109930/150‑1 (MQ=255)
          tCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCt                                                                                                                                             <  2:189561/150‑1 (MQ=255)
          ggtGCGTTAAGTCTTTTTTCAGAGGTGCTAAATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGAAAATCGCt                                                                                                                                             <  1:4325058/148‑1 (MQ=11)
            tGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGc                                                                                                                                           <  1:342374/150‑1 (MQ=255)
            tGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGc                                                                                                                                           >  1:3524680/1‑150 (MQ=255)
             gCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCg                                                                                                                                          >  2:4659043/1‑150 (MQ=255)
              cGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCgg                                                                                                                                         >  2:2368149/1‑150 (MQ=255)
              cGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCgg                                                                                                                                         >  2:293618/1‑150 (MQ=255)
              cGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCgg                                                                                                                                         >  1:1722319/1‑150 (MQ=255)
               gTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGATAATCGCTGCGGt                                                                                                                                        >  1:4121506/1‑150 (MQ=255)
               gTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGt                                                                                                                                        <  2:921744/150‑1 (MQ=255)
                ttAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGta                                                                                                                                       <  2:3225667/150‑1 (MQ=255)
                ttAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGta                                                                                                                                       >  2:3351824/1‑150 (MQ=255)
                ttAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGta                                                                                                                                       >  1:21976/1‑150 (MQ=255)
                ttAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGta                                                                                                                                       >  1:1898263/1‑150 (MQ=255)
                 tAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGtaa                                                                                                                                      <  1:560879/150‑1 (MQ=255)
                 tAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGtaa                                                                                                                                      <  2:2466442/150‑1 (MQ=255)
                  aaGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTCCCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGtaat                                                                                                                                     <  2:4086331/150‑1 (MQ=255)
                  aaGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGtaat                                                                                                                                     >  1:5157728/1‑150 (MQ=255)
                  aaGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGtaat                                                                                                                                     >  1:441343/1‑150 (MQ=255)
                   aGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGtaata                                                                                                                                    >  1:3192215/1‑150 (MQ=255)
                   aGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGtaata                                                                                                                                    >  1:1451689/1‑150 (MQ=255)
                    gACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGtaataa                                                                                                                                   >  1:736782/1‑150 (MQ=255)
                    gACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGtaataa                                                                                                                                   <  1:4372455/150‑1 (MQ=255)
                       ttttttTGTGGGTGACGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTc                                                                                                                                <  1:2888831/150‑1 (MQ=255)
                       ttttttCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTc                                                                                                                                >  1:4993376/1‑150 (MQ=255)
                       ttttttCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTc                                                                                                                                >  1:1803446/1‑150 (MQ=255)
                       ttttttCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTc                                                                                                                                <  2:2999706/150‑1 (MQ=255)
                        tttttCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTcc                                                                                                                               >  2:360880/1‑150 (MQ=255)
                        tttttCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTcc                                                                                                                               <  2:3738223/150‑1 (MQ=255)
                        tttttCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTcc                                                                                                                               <  2:2804199/150‑1 (MQ=255)
                         ttttCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTccc                                                                                                                              >  1:2687084/1‑150 (MQ=255)
                          tttCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCg                                                                                                                             <  1:658986/150‑1 (MQ=255)
                          tttCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCg                                                                                                                             >  2:3234476/1‑150 (MQ=255)
                           ttCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGc                                                                                                                            >  1:879320/1‑150 (MQ=255)
                           ttCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGc                                                                                                                            >  2:2087843/1‑150 (MQ=255)
                           ttCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGc                                                                                                                            >  2:526339/1‑150 (MQ=255)
                            tCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGcc                                                                                                                           <  2:3665524/150‑1 (MQ=255)
                             cATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCa                                                                                                                          >  1:2113430/1‑150 (MQ=255)
                             cATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCa                                                                                                                          >  1:4946621/1‑150 (MQ=255)
                             cATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCa                                                                                                                          >  2:4017229/1‑150 (MQ=255)
                               tGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAgg                                                                                                                        >  1:108450/1‑150 (MQ=255)
                               tGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAgg                                                                                                                        >  1:2067558/1‑150 (MQ=255)
                               ggggTGCCGAATAAAAAAAGGCCCGCAGGCTGTTCCCCCTGCGGCCGGTTTCGGGCGCATATTGCCACCACGGCAGCCTGACGCCCTTTTTCCCCTTCCTTCCGGTTACCCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAgg                                                                                                                        <  2:4351835/149‑1 (MQ=11)
                                 ggTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGAt                                                                                                                      >  2:2578548/1‑150 (MQ=40)
                                 ggTGCCGAATAAAAAAAAGGCCGCAGGCTGTTCCCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACCCCCCCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGAt                                                                                                                      <  2:177346/150‑1 (MQ=11)
                                  gTGCCGGATACAAAAAAGGCCGCAGGCTGTTCCCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATc                                                                                                                     <  2:4766549/150‑1 (MQ=255)
                                  gTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATc                                                                                                                     >  1:2727177/1‑150 (MQ=39)
                                  gTGCCGGATAAAAAAAAGGCCGCAGGCTTTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATc                                                                                                                     <  2:4593315/150‑1 (MQ=255)
                                   tGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCCCCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCg                                                                                                                    <  1:3633164/150‑1 (MQ=25)
                                   tGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCg                                                                                                                    >  2:1358914/1‑150 (MQ=39)
                                    gCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCgg                                                                                                                   >  1:1237165/1‑150 (MQ=38)
                                    gCCGGATACAAAAAAGCCCGCAGCCTTTTCCCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCCCGGCAGCCTGACGCCCGTTTTCACCTTCCTCCCGGTTACGCCCCCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCgg                                                                                                                   <  2:3953804/150‑1 (MQ=11)
                                     ccGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGt                                                                                                                  >  1:4103192/1‑150 (MQ=38)
                                      cGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGtg                                                                                                                 >  2:3836105/1‑150 (MQ=35)
                                      cGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGtg                                                                                                                 >  1:259145/1‑150 (MQ=35)
                                      cGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGtg                                                                                                                 >  2:2391578/1‑150 (MQ=18)
                                       ggATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGGAATTATTACCGCCAGGATCGGtgc                                                                                                                >  2:1440961/1‑150 (MQ=255)
                                       ggATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGtgc                                                                                                                <  1:4528698/150‑1 (MQ=17)
                                         aTACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGtgccg                                                                                                              >  2:2112660/1‑148 (MQ=18)
                                           aCAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTgccgcc                                                                                                            >  2:2055755/1‑150 (MQ=17)
                                            cAAAAAAGCCCGCAGGCTGTTCCCCCTGCGGCCGGTTTCGGGCGCATATTGCCACCCCGGCACCCTGCCGCCCGTTTTCCCCTTCCTTCCGGTTCCGCCCCCAGCTGCCAATCGCTGCGGTAATAATCCCCGCCAGGATCGGTGCTGCCa                                                                                                           <  2:519360/150‑1 (MQ=11)
                                                 aaGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCg                                                                                                      >  1:5200718/1‑150 (MQ=31)
                                                  aGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCCCCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGt                                                                                                     <  1:603422/150‑1 (MQ=11)
                                                                                                                                   tCCGGTTCCGCCACCGGCTGACAATCGCTGGGGAAATAATCCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAtt                    <  2:2423258/150‑1 (MQ=11)
                                                                                                                                    ccGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTg                   <  1:1030448/150‑1 (MQ=32)
                                                                                                                                    ccGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTg                   >  2:4951844/1‑150 (MQ=32)
                                                                                                                                     gggTTACGCCACCAGCTGAAAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGt                  <  1:3600610/149‑1 (MQ=17)
                                                                                                                                      ggTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGTCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGtt                 >  1:3110624/1‑150 (MQ=17)
                                                                                                                                      ggTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGtt                 <  1:2254821/150‑1 (MQ=32)
                                                                                                                                       gTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTc                >  1:4264054/1‑150 (MQ=18)
                                                                                                                                        ttACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCt               >  1:1143923/1‑150 (MQ=34)
                                                                                                                                          ccGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTgg             <  2:3144979/149‑1 (MQ=25)
                                                                                                                                           cGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGt            >  1:1529281/1‑150 (MQ=37)
                                                                                                                                           cGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGt            >  1:3477031/1‑150 (MQ=37)
                                                                                                                                            gCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGtt           <  2:3933837/150‑1 (MQ=37)
                                                                                                                                              caccaGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCa         <  1:401964/150‑1 (MQ=38)
                                                                                                                                                 caGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCAAGa      >  1:4723152/1‑150 (MQ=40)
                                                                                                                                                   gCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCAAGAtt    <  2:3696137/150‑1 (MQ=40)
                                                                                                                                                   gCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCAAGAtt    <  1:1988670/150‑1 (MQ=40)
                                                                                                                                                    cTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCAAGATTa   >  1:335974/1‑150 (MQ=40)
                                                                                                                                                    cTGACAATCGCTGCGGTAAGAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCAAGATTa   <  1:3940217/150‑1 (MQ=255)
                                                                                                                                                    cTGACAATCGCGGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCAAGATTa   <  2:4179149/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCAAGATTAG  >  AP012306/1074470‑1074768

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