Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0123061,128,4120TG23.6% 633.1 / 38.2 255V29G (GTG→GGG) yciOconserved protein
Reads supporting (aligned to +/- strand):  ref base T (55/139);  new base G (60/0);  total (116/139)
Fisher's exact test for biased strand distribution p-value = 2.33e-26
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 320 total aligned reads displayed.

CCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCAGCTGCCGGACGGTCACAACTTTACCCTGATGTGTCGCGATCTTTCTGAACTCTCGACCTATT  >  AP012306/1128263‑1128561
                                                                                                                                                     |                                                                                                                                                      
   cAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGA‑t                                                                                                                                                    >  2:3448140/1‑150 (MQ=255)
    aGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGA‑tt                                                                                                                                                   >  1:574061/1‑150 (MQ=255)
        gTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGttt                                                                                                                                               <  1:1535637/150‑1 (MQ=255)
        gTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGgg‑ttggtt                                                                                                                                               >  2:3726126/1‑150 (MQ=255)
        gTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGG‑ATGttt                                                                                                                                               >  2:3992392/1‑150 (MQ=255)
        gTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGAGttgtt                                                                                                                                                >  1:1134328/1‑150 (MQ=255)
         ttGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTa                                                                                                                                              >  1:1716828/1‑150 (MQ=255)
          tGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGA‑TTGTTTAt                                                                                                                                             >  1:1891344/1‑150 (MQ=255)
           gAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATc                                                                                                                                            <  2:3959663/150‑1 (MQ=255)
                  tCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGa                                                                                                                                     <  1:3640053/150‑1 (MQ=255)
                         gCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGA‑TTGGTTGTCCAACACATTCCCGc                                                                                                                              >  2:4146914/1‑150 (MQ=255)
                            aTGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGA‑TTGTTTATCCAACTGATTCCGGCTAt                                                                                                                           >  2:2571737/1‑150 (MQ=255)
                             tgtgAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATg                                                                                                                          >  1:2716459/4‑150 (MQ=255)
                                   cAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGG‑TTGTTTATCCAAATGGTTTCGGGTATGCGCTCg                                                                                                                    >  2:174520/1‑150 (MQ=255)
                                     gCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGc                                                                                                                  >  2:3886226/1‑150 (MQ=255)
                                      ccGCAGAACACCACAGAGAGGGAATTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCt                                                                                                                 <  2:3866272/150‑1 (MQ=255)
                                        gCAGAACACCACAGAGAGGGAATTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGt                                                                                                               >  1:1984516/1‑150 (MQ=255)
                                           gAACACCACAGAGAGGGAATTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTTCGGCTATTCGCTCGGCTGGaaa                                                                                                            >  2:1856453/1‑150 (MQ=255)
                                              caccacAGAGAGGGAATTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAAtt                                                                                                         >  1:3738265/1‑150 (MQ=255)
                                               accacAGAGAGGGAATTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTg                                                                                                        >  1:2851710/1‑150 (MQ=255)
                                                ccacAGAGAGGGAATTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGA‑TTGTTTATACTACATTTTTTGGGTTAGCGCTCTGCTGTAGAAttta                                                                                                       >  2:2749238/1‑148 (MQ=255)
                                                 cacaGAGAGGGAATTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGaa                                                                                                      >  2:3405652/1‑150 (MQ=255)
                                                  acaGAGAGGGAATTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAg                                                                                                     >  1:2189873/1‑150 (MQ=255)
                                                   caGAGAGGGAATTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGA‑TTGTTTTTCCAACTGATTCCCGCTATGCGCTCGGCTGTAAAATTGAAGa                                                                                                    >  2:2485003/1‑150 (MQ=255)
                                                    agagagGGAATTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGA‑TTGTTTATCCAACCGATTCCGGGTATGCGCTCGGCTGTAAAATTGAAGAt                                                                                                   >  2:4110269/1‑150 (MQ=255)
                                                     gagagGGAATTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGGAGATa                                                                                                  >  1:2990027/1‑150 (MQ=255)
                                                     gagagGGAATTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATa                                                                                                  <  1:883792/150‑1 (MQ=255)
                                                      agagGGAATTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATaa                                                                                                 >  1:3175744/1‑150 (MQ=255)
                                                       gagGGAATTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATaaa                                                                                                <  2:4247295/150‑1 (MQ=255)
                                                        agGGAATTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATaaaa                                                                                               >  1:648587/1‑150 (MQ=255)
                                                         gggAATTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGA‑TTTTTTATCCCACTTATTCCGGCTATGCGCTCGGCTGTAAAATTTAAGAaaaaaa                                                                                              >  2:1301333/1‑150 (MQ=255)
                                                          ggAATTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAAc                                                                                             <  1:3466671/150‑1 (MQ=255)
                                                          ggAATTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGA‑TTTTTTTTCCAACTGTTTTCGGCTATGCGCTCGGCTGGTAAATTTAAGATaaaaaa                                                                                             >  2:2665244/1‑149 (MQ=255)
                                                            aaTTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGA‑TTGTTTATCCAACTGATTCCGGGTATGCGCTCGGCTGTAAAATTTGAGATAAAAAagc                                                                                           >  2:33725/1‑150 (MQ=255)
                                                             aTTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAAcgcg                                                                                          <  2:1914061/150‑1 (MQ=255)
                                                             aTTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGG‑TGGTTTATGCGACCGAGTCCGGCTATGCGGTCGGGGGTAAAAATTAAGGTAAAAAggcg                                                                                          >  2:460766/1‑150 (MQ=255)
                                                             aTTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGA‑TTGTTTTTCCAACTGATTCCGGCTATGCGCTCGGCTGGAAAATTTAAGATAAAAAagcg                                                                                          >  2:2228243/1‑150 (MQ=255)
                                                             aTTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAAcgcg                                                                                          >  1:732417/1‑150 (MQ=255)
                                                              ttATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGa                                                                                         <  1:938722/150‑1 (MQ=255)
                                                               tATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGGAAAATTGAAGATAAAAACGCGAt                                                                                        >  1:98603/1‑150 (MQ=255)
                                                                aTGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATg                                                                                       <  1:811877/150‑1 (MQ=255)
                                                                aTGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGG‑TTGTTTATCCAACTGATTCCGGCTATGGGCTCGGGTGTAAAATTGAAGAAAAAAAAGCGaat                                                                                       >  2:898615/1‑148 (MQ=255)
                                                                 tGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATgg                                                                                      <  2:4197153/150‑1 (MQ=255)
                                                                  gAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGa                                                                                     <  2:3820895/150‑1 (MQ=255)
                                                                         tttttttATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGGTGGAGCGTAtt                                                                              >  1:2394137/1‑150 (MQ=255)
                                                                         tttttttATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTAtt                                                                              <  2:2983715/150‑1 (MQ=255)
                                                                         tttttttATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTAtt                                                                              <  2:3599288/150‑1 (MQ=255)
                                                                          ttttttATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTAttt                                                                             <  2:3583410/150‑1 (MQ=255)
                                                                           tttttATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGA‑TTGTTTTTCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGTAGATAAAAAAGCGATGGGGCGGAtttt                                                                            >  2:691327/1‑149 (MQ=255)
                                                                            ttttATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGt                                                                           >  1:3034578/1‑150 (MQ=255)
                                                                             tttATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTc                                                                          <  1:1213227/150‑1 (MQ=255)
                                                                             tttATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGG‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTTAAGATAAAAACGCGATGGGGCGTATTTGTc                                                                          >  2:2051674/1‑150 (MQ=255)
                                                                             tttATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGA‑TTGTTTATCCAACTGATTCCGGCCATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTc                                                                          >  1:446142/1‑150 (MQ=255)
                                                                                atatTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTa                                                                       <  1:1656042/150‑1 (MQ=255)
                                                                                atatTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTa                                                                       <  2:2009025/150‑1 (MQ=255)
                                                                                 tatTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGG‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGGAAAATTGAAGATAAAAACGCGATGGAGCGTATTTTTCGTAt                                                                      >  2:1715494/1‑150 (MQ=255)
                                                                                   ttCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTc                                                                    >  1:3578758/1‑150 (MQ=255)
                                                                                     cATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGc                                                                  <  1:3375720/150‑1 (MQ=255)
                                                                                     cATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGA‑TTGTTTATCCAACTGATTCCGGCTATTCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTTTCGGATTTGc                                                                  >  2:1233156/1‑150 (MQ=255)
                                                                                        ccTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCAg                                                               >  1:497692/1‑150 (MQ=255)
                                                                                        ccTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCCGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTTTCGGATTCGCCac                                                               >  2:2822094/1‑149 (MQ=255)
                                                                                           gATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCAGCTg                                                            <  1:2289452/150‑1 (MQ=255)
                                                                                           gATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCAGCTg                                                            <  1:1365823/150‑1 (MQ=255)
                                                                                            aTAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGGGCGTATTTGTCGGATTCGCCAGCTGc                                                           >  2:1205242/1‑150 (MQ=255)
                                                                                                   acaGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCAGCTGCCGGACgg                                                    <  1:1090035/150‑1 (MQ=255)
                                                                                                   acaGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGA‑TTGTTTATCCAACTGATTCCGGCTATTCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTTTCGTATTCGCCATCTGCCggcagg                                                    >  2:574006/1‑146 (MQ=255)
                                                                                                      gCAACTTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCAGCTGCCGGACGGTca                                                 >  1:3249050/1‑150 (MQ=255)
                                                                                                       cAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTTTCCAACTGATTCCGGGTATGCGCTCGGCTGTAAAAATGAAGATAAAAACGCGATGGAGCGTATTTTTCGTTTTTGCCAGCTGCCGGACGGGcac                                                >  2:238825/1‑150 (MQ=255)
                                                                                                         aCGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCAGCTGCCGGACGGTCACaa                                              <  2:2321417/150‑1 (MQ=255)
                                                                                                         aCGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCAGCTGCCGGACGGTCACaa                                              >  1:2913160/1‑150 (MQ=255)
                                                                                                          cGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCAGCTGCCGGACGGTCACAAc                                             >  2:1457861/1‑150 (MQ=255)
                                                                                                           gTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCAGCTGCCGGACGGTCAGAACt                                            <  2:3147035/150‑1 (MQ=255)
                                                                                                           gTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCAGCTGCCGGACGGTCACAACt                                            <  2:800066/150‑1 (MQ=255)
                                                                                                           gTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCAGCTGCCGGACGGTCACAACt                                            <  2:1942189/150‑1 (MQ=255)
                                                                                                           gTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGG‑TTGTTTATCCCACTGATTTCGGCTATGCGCTCGGCTGTAAAATTTGAGATAAAAACGCGATGGAGCGTATTTTTCGTATTTGCCCGCTGCCGGGCGGGCACAACt                                            >  2:2860113/1‑150 (MQ=255)
                                                                                                               gATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCAGCTGCCGGACGGTCACAACTTTAc                                        <  1:2439591/150‑1 (MQ=255)
                                                                                                                   aaCCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGA‑TTGTTTTTCCAAATGATTCCGGCTATGGGCTCGGCTGTAAAATTGGAGATAAAAAAGCGATGGAGCGGATTTGTCGGATTCCCCAGCTGCCGGACGGTCACAACTTTAACCtc                                    >  2:2887815/1‑149 (MQ=255)
                                                                                                                    aCCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGG‑TTGTTTATCCAACTGCTTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTTTTTGTCGTTTTCGCCAGCTGCCGGACGGGCACAACTTTACCCTGa                                   >  2:2646649/1‑150 (MQ=255)
                                                                                                                       aGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCAGCTGCCGGACGGTCACAACTTTACCCTGAtgt                                <  1:390905/150‑1 (MQ=255)
                                                                                                                       aGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGA‑TTGTTTATCCCACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCCGCTGCCGGACGGCCACCACTTTACCCTGAtgg                                >  2:2170039/1‑149 (MQ=255)
                                                                                                                        ggCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCAGCTGCCGGACGGTCACAACTTTACCCTGAtgtg                               <  1:2812108/150‑1 (MQ=255)
                                                                                                                           ggTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCAGCTGCCGGACGGTCACAACTTTACCCTGATGTGTcg                            >  1:2734336/1‑150 (MQ=255)
                                                                                                                            gTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCAGCTGCCGGACGGTCACAACTTTACCCTGATGTGTcgc                           <  1:4200380/150‑1 (MQ=255)
                                                                                                                              ggAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCAGCTGCCGGACGGTCACAACTTTACCCTGATGTGTCGCGa                         <  2:3196286/150‑1 (MQ=255)
                                                                                                                               gagaTCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCAGCTGCCGGACGGTCACAACTTTACCCTGATGTGTCGCGAt                        <  1:1352311/150‑1 (MQ=255)
                                                                                                                                 gaTCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCAGCTGCCGGACGGTCACAACTTTACCCTGATGTGTCGCGATCt                      <  1:2077488/150‑1 (MQ=255)
                                                                                                                                   tCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCAGCTGCCGGACGGTCACAACTTTACCCTGATGTGTCGCGATCttt                    <  2:3762984/150‑1 (MQ=255)
                                                                                                                                   tCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCAGCTGCCGGACGGTCACAACTTTACCCTGATGTGTCGCGATCttt                    <  2:2399745/150‑1 (MQ=255)
                                                                                                                                    cGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCAGCTGCCGGACGGTCACAACTTTACCCTGATGTGTCGCGATCTTTc                   >  1:3429762/1‑150 (MQ=255)
                                                                                                                                    cGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCAGCTGCCGGACGGTCACAACTTTACCCTGATGTGTCGCGATCTTTc                   <  1:405933/150‑1 (MQ=255)
                                                                                                                                     gTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGGGTATTTGTCGTATTCGCCAGCTGCCGGACGGTCACAACTTTACCCTGATGTGTCGCGATCTTTCt                  >  2:3541065/1‑150 (MQ=255)
                                                                                                                                      tGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCAGCTGCCGGACGGTCACAACTTTACCCTGATGTGTCGCGATCTTTCTg                 <  2:739852/150‑1 (MQ=255)
                                                                                                                                         gTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCAGCTGCCGGACGGTCACAACTTTACCCTGATGTGTCGCGATCTTTCTGAAc              >  2:1538643/1‑150 (MQ=255)
                                                                                                                                         gTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCAGCTGCCGGACGGTCACAACTTTACCCTGATGTGTCGCGATCTTTCTGAAc              <  1:3839796/150‑1 (MQ=255)
                                                                                                                                            aaGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCAGCTGCCGGACGGTCACAACTTTACCCTGATGTGTCGCGATCTTTCTGAActct           >  1:2227828/1‑150 (MQ=255)
                                                                                                                                             aGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCAGCTGCCGGACGGTCACAACTTTACCCTGATGTGTCGCGATCTTTCTGAACtctc          <  1:3551473/150‑1 (MQ=255)
                                                                                                                                              ggCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCAGCTGCCGGACGGTCACAACTTTACCCTGATGTGTCGCGATCTTTCTGAACTCTCg         >  1:2422590/1‑150 (MQ=255)
                                                                                                                                                    gTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCAGCTGCCGGACGGTCACAACTTTACCCTGATGTGTCGCGATCTTTCTGAACTCTCGACCTAt   <  1:2172185/150‑1 (MQ=255)
                                                                                                                                                     tGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCAGCTGCCGGACGGTCACAACTTTACCCTGATGTGTCGCGATCTTTCTGAACTCTCGACCTAtt  <  2:3043800/150‑1 (MQ=255)
                                                                                                                                                     tGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCAGCTGCCGGACGGTCACAACTTTACCCTGATGTGTCGCGATCTTTCTGAACTCTCGACCTAtt  <  2:3106537/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGA‑TTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCAGCTGCCGGACGGTCACAACTTTACCCTGATGTGTCGCGATCTTTCTGAACTCTCGACCTATT  >  AP012306/1128263‑1128561

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