Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,915,397 A→C 7.2% N5K (AAT→AAG atpB ← F0 sector of membrane‑bound ATP synthase, subunit a

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,915,3970AC 7.2% 934.8 / 53.1 292N5K (AAT→AAGatpBF0 sector of membrane‑bound ATP synthase, subunit a
Reads supporting (aligned to +/- strand):  ref base A (128/143);  new base C (14/7);  total (142/150)
Fisher's exact test for biased strand distribution p-value = 1.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

Only 100 of 301 total aligned reads displayed.

AACAGACCCAGCACCACCGAGAAGAACATGGAGTCAATATTGATTGTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAACGCCACCACCAGTAACACCAACATC  >  CP009273/3915248‑3915546
                                                                                                                                                     |                                                                                                                                                     
aacaGACCCAGCACCACCGAGAAGAACATGGAGTCAATATTGATTGTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCata                                                                                                                                                       >  1:5101758/1‑150 (MQ=255)
aacaGACCCAGCACCACCGAGAAGAACATGGAGTCAATATTGATTGTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCata                                                                                                                                                       >  1:2511828/1‑150 (MQ=255)
   acccccAGCACCACCGAGAAGAACATGGAGTCAATATTGATTGTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATAttt                                                                                                                                                    <  2:1933204/147‑1 (MQ=255)
   aGACCCAGCACCACCGAGAAGAACATGGAGTCAATATTGATTGTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATAttt                                                                                                                                                    <  1:2650464/150‑1 (MQ=255)
    gACCCAGCACCACCGAGAAGAACATGGAGTCAATATTGATTGTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATAtttt                                                                                                                                                   <  1:3755160/150‑1 (MQ=255)
      cccAGCACCACCGAGAAGAACATGGAGTCAATATTGATTGTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCATGCGGCGTCATATTTTCt                                                                                                                                                 >  2:4052622/1‑150 (MQ=255)
       ccAGCACCACCGAGAAGAACATGGAGTCAATATTGATTGTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTg                                                                                                                                                <  2:4687203/150‑1 (MQ=255)
          gCACCACCGAGAAGAACATGGAGTCAATATTGATTGTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAg                                                                                                                                             <  1:3261888/150‑1 (MQ=255)
            accaccGAGAAGCACATGGAGTCAATATTGATTGTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGcc                                                                                                                                           <  2:479265/150‑1 (MQ=255)
            accaccGAGAAGAACATGGAGTCAATATTGATTGTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGcc                                                                                                                                           <  2:1392313/150‑1 (MQ=255)
            accaccGAGAAGAACATGGAGTCAATATTGATTGTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGcc                                                                                                                                           >  2:4418196/1‑150 (MQ=255)
              caccGAGAAGAACATGGAGTCAATATTGATTGTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCat                                                                                                                                         >  1:4901182/1‑150 (MQ=255)
              caccGAGAAGAACATGGAGTCAATATTGATTGTCCAGAAGGTGGCTGGGGGGGTTTGTGGATCCCCCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTACTCCTGCGGCGGCATATGTTCTGAAGCCat                                                                                                                                         >  2:85956/1‑150 (MQ=255)
                   agaagaACATGGAGTCAATATTGATTGTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGc                                                                                                                                    >  2:3572880/1‑150 (MQ=255)
                   agaagaACATGGAGTCAATATTGATTGTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGc                                                                                                                                    <  1:3060496/150‑1 (MQ=255)
                    gaagaaCATGGAGTCAATATTGATTGTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGcc                                                                                                                                   <  1:2970570/150‑1 (MQ=255)
                     aagaaCATGGAGTCAATATTGATTGTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCt                                                                                                                                  >  2:2630863/1‑150 (MQ=255)
                     aagaaCATGGAGTCAATATTGATTGTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCt                                                                                                                                  <  1:1231746/150‑1 (MQ=255)
                       gaaCATGGAGTCAATATTGATTGTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCttt                                                                                                                                >  1:4497825/1‑150 (MQ=255)
                        aaCATGGAGTCAATATTGATTGTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCtttt                                                                                                                               <  1:1469498/150‑1 (MQ=255)
                         aCATGGAGTCAATATTGATTGTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTa                                                                                                                              >  2:791271/1‑150 (MQ=255)
                          cATGGAGTCAATATTGATTGTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTAc                                                                                                                             >  1:1936908/1‑150 (MQ=255)
                          cATGGAGTCAATATTGATTGTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTAc                                                                                                                             <  2:1316631/150‑1 (MQ=255)
                           aTGGAGTCAATATTGATTGTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTAcc                                                                                                                            <  1:2051721/150‑1 (MQ=255)
                            tGGAGTCAATATTGATTGTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTAccc                                                                                                                           >  2:882398/1‑150 (MQ=255)
                              gAGTCAATATTGATTGTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCtt                                                                                                                         <  1:4961737/150‑1 (MQ=255)
                                gTCAATATTGATTGTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTttg                                                                                                                       >  2:2731925/1‑150 (MQ=255)
                                 tCAATATTGATTGTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCATGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTttgt                                                                                                                      <  2:4629821/150‑1 (MQ=255)
                                  cAATATTGATTGTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGACTTTTACCCTttgtt                                                                                                                     >  2:5079994/1‑150 (MQ=255)
                                   aaTATTGATTGTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGTGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTttgttg                                                                                                                    <  2:291110/150‑1 (MQ=255)
                                    atatTGATTGTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATCTTTTCTGAAGCCATGATGCCTTTTACCCTTtgttgt                                                                                                                   >  1:2170484/1‑150 (MQ=255)
                                      atttattgTCCAGAAGGTGGCTGGGGGTTTTTGTGGATCCACCTGCGAGAATGTACGCATGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTa                                                                                                                 <  1:5199705/150‑1 (MQ=255)
                                      attgattgTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTa                                                                                                                 >  2:1169074/1‑150 (MQ=255)
                                      attgattgTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTa                                                                                                                 >  2:545327/1‑150 (MQ=255)
                                      attgattgTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTa                                                                                                                 <  2:293834/150‑1 (MQ=255)
                                          attgTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTa                                                                                                             >  2:3147176/1‑150 (MQ=255)
                                            tgTCCTGAAGGTGGCTGGGGGGTTTTGTGGATCCTCCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCTGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACa                                                                                                           <  2:2814349/150‑1 (MQ=255)
                                             gTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAg                                                                                                          <  2:687840/150‑1 (MQ=255)
                                              tCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGc                                                                                                         >  2:2334376/1‑150 (MQ=255)
                                               ccAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGcc                                                                                                        <  1:1042530/150‑1 (MQ=255)
                                                cAGAAGGTTGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCg                                                                                                       <  2:5111022/150‑1 (MQ=255)
                                                  gatggTGGCTGGGGGGTTTTGTGGATCCACAAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGt                                                                                                     <  1:1731714/147‑1 (MQ=255)
                                                        ggCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGt                                                                                               <  2:442531/150‑1 (MQ=255)
                                                         gCTGGGGGGTTGTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATCTTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTa                                                                                              <  2:4849872/150‑1 (MQ=255)
                                                          cTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTAt                                                                                             <  1:4217875/150‑1 (MQ=255)
                                                            ggggggTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTCATTACAGCCTGTGCCAGTATCt                                                                                           >  2:134425/1‑150 (MQ=255)
                                                                 gTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAacc                                                                                      <  1:1411031/150‑1 (MQ=255)
                                                                    ttGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATCTTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAaccacc                                                                                   >  1:2883359/1‑150 (MQ=255)
                                                                    ttGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAaccacc                                                                                   <  1:1784933/150‑1 (MQ=255)
                                                                         gATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGAAc                                                                              >  2:3934421/1‑150 (MQ=255)
                                                                          aTCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCAcc                                                                             >  2:422457/1‑150 (MQ=255)
                                                                           tCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGGTAATCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTATGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCa                                                                            >  2:1282855/1‑150 (MQ=255)
                                                                            ccaccaGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCaa                                                                           <  1:114511/150‑1 (MQ=255)
                                                                            ccaccaGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCaa                                                                           <  1:4941271/150‑1 (MQ=255)
                                                                               ccaGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAAc                                                                        <  1:2523379/150‑1 (MQ=255)
                                                                               ccaGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAAc                                                                        <  1:3131282/150‑1 (MQ=255)
                                                                                caGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAAcc                                                                       <  1:1498180/150‑1 (MQ=255)
                                                                                 aGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAAccc                                                                      <  1:1885084/150‑1 (MQ=255)
                                                                                  gCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCa                                                                     <  2:3381261/150‑1 (MQ=255)
                                                                                    gagaATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATCTTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACg                                                                   >  1:1379191/1‑150 (MQ=255)
                                                                                      gaATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTa                                                                 >  1:3710080/1‑150 (MQ=255)
                                                                                        aTGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTacc                                                               <  1:4181888/150‑3 (MQ=255)
                                                                                        aTGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAAc                                                               <  1:770937/150‑1 (MQ=255)
                                                                                          gTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGa                                                             >  2:4521921/1‑150 (MQ=255)
                                                                                            aCGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATc                                                           <  2:3072840/150‑1 (MQ=255)
                                                                                              gCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAg                                                         <  2:833254/150‑1 (MQ=255)
                                                                                               cAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGc                                                        <  1:4418273/150‑1 (MQ=255)
                                                                                                aGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCg                                                       >  2:1662539/1‑150 (MQ=255)
                                                                                                    ccAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCaa                                                   >  1:889691/1‑150 (MQ=255)
                                                                                                      aGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGa                                                 <  2:2489536/150‑1 (MQ=255)
                                                                                                      aGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGa                                                 <  1:542856/150‑1 (MQ=255)
                                                                                                        cTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCACGAAt                                               >  1:3996292/1‑150 (MQ=255)
                                                                                                        cTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAAt                                               <  2:1993909/150‑1 (MQ=255)
                                                                                                         tGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATa                                              >  2:3440255/1‑150 (MQ=255)
                                                                                                          gAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATAc                                             >  2:2226740/1‑150 (MQ=255)
                                                                                                           aaGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATAcc                                            >  2:5080648/1‑150 (MQ=255)
                                                                                                            aGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCg                                           >  1:3099663/1‑150 (MQ=255)
                                                                                                              gTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGcc                                         >  2:4283584/1‑150 (MQ=255)
                                                                                                                     aGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTACCGATCAGCGGCAAGAATACCGCCTTTaatc                                  >  2:4928856/1‑148 (MQ=255)
                                                                                                                     aGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTaaaa                                  >  2:4131704/1‑150 (MQ=255)
                                                                                                                      ggtggtGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAAc                                 <  1:1622276/150‑1 (MQ=255)
                                                                                                                      ggtggtGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAAc                                 <  2:1719799/150‑1 (MQ=255)
                                                                                                                      ggtggtGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAAc                                 >  1:4690589/1‑150 (MQ=255)
                                                                                                                       gtggtgTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAAcc                                <  2:770055/150‑1 (MQ=255)
                                                                                                                           tgtCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCa                            >  2:2704943/1‑150 (MQ=255)
                                                                                                                                tATGTAATCCTGCGGCGTCATCTTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAACGcc                       >  1:273842/1‑150 (MQ=255)
                                                                                                                                  tGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAACGccac                     >  2:2471225/1‑150 (MQ=255)
                                                                                                                                    tAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAACGccacca                   >  1:905898/1‑150 (MQ=255)
                                                                                                                                    tAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAACGccacca                   >  1:3578645/1‑150 (MQ=255)
                                                                                                                                       tacTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAACGccaccacca                <  2:3681223/148‑1 (MQ=255)
                                                                                                                                        ccTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAATCCGCCAACGCCACCACCAg               >  2:3630222/1‑150 (MQ=255)
                                                                                                                                        ccTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAACGCCACCACCAg               >  1:2796678/1‑150 (MQ=255)
                                                                                                                                          tGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAACGCCACCACCAGTa             >  1:3935168/1‑150 (MQ=255)
                                                                                                                                           gcggcgTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAACGCCACCACCAGTaa            <  2:3853052/150‑1 (MQ=255)
                                                                                                                                            cggcgTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAAGCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAACGCCACCACCAGTAac           <  2:458000/150‑1 (MQ=255)
                                                                                                                                             ggcgTCATCTTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTACCGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAACGCCACCACCAGTAaca          >  2:1031132/1‑150 (MQ=255)
                                                                                                                                             ggcgTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAACGCCACCACCAGTAaca          >  1:3620656/1‑150 (MQ=255)
                                                                                                                                             ggcgTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAACGCCACCACCAGTAaca          >  1:1348918/1‑150 (MQ=255)
                                                                                                                                               cgTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAACGCCACCACCAGTAACAcc        >  2:2986947/1‑150 (MQ=255)
                                                                                                                                                   atatTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCACGAATACCGCCTTTAAAACCGCCAACGCCACCACCAGTAACACCAACa    >  2:522258/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AACAGACCCAGCACCACCGAGAAGAACATGGAGTCAATATTGATTGTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAACGCCACCACCAGTAACACCAACATC  >  CP009273/3915248‑3915546

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: