Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,837,532 C→A 68.0% P70Q (CCG→CAG)  ade → cryptic adenine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,837,5320CA68.0% 343.1 / inf 368P70Q (CCG→CAG) adecryptic adenine deaminase
Reads supporting (aligned to +/- strand):  ref base C (61/57);  new base A (109/141);  total (170/198)
Fisher's exact test for biased strand distribution p-value = 1.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 369 total aligned reads displayed.

GTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACC  >  CP009273/3837383‑3837678
                                                                                                                                                     |                                                                                                                                                  
gTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTcc                                                                                                                                                    <  1:2790417/150‑1 (MQ=255)
    gCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCTGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCt                                                                                                                                                <  2:1875595/150‑1 (MQ=255)
    gCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCt                                                                                                                                                <  2:4107406/150‑1 (MQ=255)
     ccGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCtt                                                                                                                                               <  2:1715426/150‑1 (MQ=255)
      cGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCttt                                                                                                                                              <  1:1096702/150‑1 (MQ=255)
       gtgtCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTg                                                                                                                                             <  2:1029566/147‑1 (MQ=255)
       gTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTg                                                                                                                                             >  2:2324917/1‑150 (MQ=255)
       gTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTg                                                                                                                                             >  1:3570446/1‑150 (MQ=255)
       gTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTg                                                                                                                                             <  2:4764409/150‑1 (MQ=255)
           cccGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGc                                                                                                                                         <  1:5223244/150‑1 (MQ=255)
            ccGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATATCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCg                                                                                                                                        >  1:214904/1‑150 (MQ=255)
                   gACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGAt                                                                                                                                 <  2:3597940/150‑1 (MQ=255)
                     cGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGc                                                                                                                               <  2:93708/150‑1 (MQ=255)
                        tGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTcg                                                                                                                            <  2:2405350/150‑1 (MQ=255)
                        tGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTcg                                                                                                                            >  2:4027506/1‑150 (MQ=255)
                          ttGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTcgcg                                                                                                                          <  2:67208/150‑1 (MQ=255)
                           tGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCgcgg                                                                                                                         >  1:2629097/1‑150 (MQ=255)
                              cGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGgcgc                                                                                                                      >  2:2388951/1‑150 (MQ=255)
                               gATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCa                                                                                                                     >  2:3476290/1‑150 (MQ=255)
                                     attattGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAAcggcg                                                                                                               >  1:1386252/1‑150 (MQ=255)
                                     attattGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAAcggcg                                                                                                               >  1:2389059/1‑150 (MQ=255)
                                      ttatttatAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATAGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAAcggcgg                                                                                                              >  2:2795241/1‑150 (MQ=255)
                                          tGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGcc                                                                                                          >  1:404941/1‑150 (MQ=255)
                                           gATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCa                                                                                                         <  2:1326603/150‑1 (MQ=255)
                                               aTGTCTCTATTCTCGACCTGAACAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGt                                                                                                     >  1:3712776/1‑150 (MQ=255)
                                                 gTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGttt                                                                                                   >  1:1597703/1‑150 (MQ=255)
                                                  tctcTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTa                                                                                                  <  1:507810/150‑1 (MQ=255)
                                                   ctctctTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTAt                                                                                                 <  1:3973133/150‑1 (MQ=255)
                                                    tctATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTAtt                                                                                                <  2:364474/150‑1 (MQ=255)
                                                     ctATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTg                                                                                               >  1:667657/1‑150 (MQ=255)
                                                     ctATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTg                                                                                               >  2:452706/1‑150 (MQ=255)
                                                      tATTCTCTACCTGATCATTGGCTGAGTTATTTCCGGCACTTTTGTGTTTATAGGACGTTCCTTTGCCGGTGTTGGCTCAGAATACACTGATTCTCAGGCTTTGCAGCGGATTGTTGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGa                                                                                              <  1:2290221/150‑1 (MQ=255)
                                                      tATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGa                                                                                              >  2:807530/1‑150 (MQ=255)
                                                         tctcGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGc                                                                                           >  2:3769194/1‑150 (MQ=255)
                                                           tcGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTc                                                                                         >  2:4651379/1‑150 (MQ=255)
                                                            cGACCTGATCTATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCa                                                                                        <  2:2244002/150‑1 (MQ=255)
                                                                cTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTg                                                                                    <  1:907197/150‑1 (MQ=255)
                                                                      aaTGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATAtt                                                                              >  2:3263096/1‑150 (MQ=255)
                                                                      aaTGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATAtt                                                                              <  1:4271900/150‑1 (MQ=255)
                                                                       aTGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTg                                                                             <  1:2035022/150‑1 (MQ=255)
                                                                       aTGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTg                                                                             <  1:3430925/150‑1 (MQ=255)
                                                                         ggCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGaa                                                                           <  1:2591250/150‑1 (MQ=255)
                                                                          gCGGAGATATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAAt                                                                          <  2:287256/150‑1 (MQ=255)
                                                                             gagaAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCa                                                                       <  2:5120401/150‑1 (MQ=255)
                                                                                 aaTTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCat                                                                   >  1:4514734/1‑150 (MQ=255)
                                                                                  aTTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCatg                                                                  <  2:1538203/150‑1 (MQ=255)
                                                                                  aTTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCatg                                                                  >  2:2620458/1‑150 (MQ=255)
                                                                                    ttCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCatgat                                                                <  1:4972006/150‑1 (MQ=255)
                                                                                     tCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCatgatg                                                               <  1:3395875/150‑1 (MQ=255)
                                                                                       cGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGAc                                                             <  1:1705696/150‑1 (MQ=255)
                                                                                          cccAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGcc                                                          <  2:2106859/150‑1 (MQ=255)
                                                                                            cAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCgg                                                        >  1:3057663/1‑150 (MQ=255)
                                                                                              aTTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTc                                                      >  2:72561/1‑150 (MQ=255)
                                                                                              aTTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTc                                                      >  1:1776853/1‑150 (MQ=255)
                                                                                              aTTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTc                                                      >  1:1861047/1‑150 (MQ=255)
                                                                                                 gtgATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACt                                                   <  2:1141695/150‑1 (MQ=255)
                                                                                                 gtgATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACt                                                   <  2:4127863/150‑1 (MQ=255)
                                                                                                  tgATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACtt                                                  >  2:4186246/1‑150 (MQ=255)
                                                                                                   gATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACttt                                                 <  2:125277/150‑1 (MQ=255)
                                                                                                   gATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACttt                                                 <  1:4419849/150‑1 (MQ=255)
                                                                                                    aTTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACtttt                                                <  2:1822668/150‑1 (MQ=255)
                                                                                                      tAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGa                                              <  1:4746445/150‑1 (MQ=255)
                                                                                                      tAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGa                                              <  1:787482/150‑1 (MQ=255)
                                                                                                       aaaGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGaa                                             >  2:1402180/1‑150 (MQ=255)
                                                                                                       aaaGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGaa                                             >  2:4369511/1‑150 (MQ=255)
                                                                                                       aaaGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGaa                                             >  1:5202779/1‑150 (MQ=255)
                                                                                                         aGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAAc                                           <  1:4676298/150‑1 (MQ=255)
                                                                                                          ggACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAAcc                                          <  1:4027506/150‑1 (MQ=255)
                                                                                                          ggACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAAcc                                          >  1:4773849/1‑150 (MQ=255)
                                                                                                           gACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCg                                         >  1:840104/1‑150 (MQ=255)
                                                                                                             cGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCt                                       <  1:643605/150‑1 (MQ=255)
                                                                                                             cGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCt                                       <  2:3101334/150‑1 (MQ=255)
                                                                                                             cGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCt                                       <  2:3375003/150‑1 (MQ=255)
                                                                                                                tACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTAcc                                    >  1:4917034/1‑150 (MQ=255)
                                                                                                                tACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTAcc                                    >  2:1930161/1‑150 (MQ=255)
                                                                                                                 aCATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTAccc                                   >  1:1569063/1‑150 (MQ=255)
                                                                                                                 aCATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTAccc                                   >  1:2213720/1‑150 (MQ=255)
                                                                                                                  cATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCt                                  >  2:3697316/1‑150 (MQ=255)
                                                                                                                  cATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCt                                  <  2:2540745/150‑1 (MQ=255)
                                                                                                                  cATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCt                                  >  1:97707/1‑150 (MQ=255)
                                                                                                                   aTTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTg                                 <  2:4855179/150‑1 (MQ=255)
                                                                                                                         ggTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCcgcgc                           >  2:2024840/1‑150 (MQ=255)
                                                                                                                            gtTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGc                        >  2:3002481/1‑150 (MQ=255)
                                                                                                                             ttGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGcc                       <  1:2340754/150‑1 (MQ=255)
                                                                                                                              tGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCt                      <  1:3793514/150‑1 (MQ=255)
                                                                                                                               ggCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTg                     <  2:624591/150‑1 (MQ=255)
                                                                                                                                 cgcAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTgac                   <  1:4482885/150‑1 (MQ=255)
                                                                                                                                 cgcAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTgac                   >  1:3268018/1‑150 (MQ=255)
                                                                                                                                  gcAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTgacg                  <  2:5001873/150‑1 (MQ=255)
                                                                                                                                   cAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTgacga                 >  2:1436455/1‑150 (MQ=255)
                                                                                                                                   cAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCAGATTGAATCCAGCATGATGACGCCGGTCACGTATGAATCCGCTACCCTGCCGCGCGGCCTgacga                 >  2:3534398/1‑150 (MQ=255)
                                                                                                                                    aGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTgacgac                >  1:1055277/1‑150 (MQ=255)
                                                                                                                                     gAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGAcc               >  2:4964050/1‑150 (MQ=255)
                                                                                                                                     gAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGAcc               <  1:2430879/150‑1 (MQ=255)
                                                                                                                                     gAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGAcc               >  2:3165885/1‑150 (MQ=255)
                                                                                                                                         acacTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTa           >  1:2807556/1‑150 (MQ=255)
                                                                                                                                          cacTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTAt          <  1:2906048/150‑1 (MQ=255)
                                                                                                                                           acTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTAtt         >  1:1165954/1‑150 (MQ=255)
                                                                                                                                            cTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTAttt        <  2:5150084/150‑1 (MQ=255)
                                                                                                                                              gATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGc      >  1:757069/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACC  >  CP009273/3837383‑3837678

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