Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,365,607 G→T 69.4% I109I (ATC→ATA)  nanT ← sialic acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,365,6070GT69.4% 279.1 / 227.7 222I109I (ATC→ATA) nanTsialic acid transporter
Reads supporting (aligned to +/- strand):  ref base G (40/28);  new base T (80/74);  total (120/102)
Fisher's exact test for biased strand distribution p-value = 3.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.21e-01

Only 100 of 224 total aligned reads displayed.

GACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCAGAGATCAG  >  CP009273/3365460‑3365756
                                                                                                                                                   |                                                                                                                                                     
    gCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAg                                                                                                                                                 >  2:2265860/1‑150 (MQ=255)
     cccccACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGc                                                                                                                                                <  2:2578041/150‑1 (MQ=255)
            agagAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGgcgc                                                                                                                                         >  2:1940345/1‑150 (MQ=255)
            agagAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACAGGGTGATGTAGCCTGgcgc                                                                                                                                         >  2:1940043/1‑150 (MQ=255)
             gagaAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCa                                                                                                                                        <  1:2261377/150‑1 (MQ=255)
                 aGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGc                                                                                                                                    <  1:2343308/150‑1 (MQ=255)
                 aGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGc                                                                                                                                    <  1:871872/150‑1 (MQ=255)
                 aGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGc                                                                                                                                    >  2:351044/1‑150 (MQ=255)
                 aGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACAGAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGc                                                                                                                                    >  2:350218/1‑150 (MQ=255)
                  gCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGcc                                                                                                                                   >  1:3180647/1‑150 (MQ=255)
                  gCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGcc                                                                                                                                   <  1:4926307/150‑1 (MQ=255)
                   ccTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCg                                                                                                                                  >  2:107664/1‑150 (MQ=255)
                    cTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGc                                                                                                                                 >  2:1807458/1‑150 (MQ=255)
                      gAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAg                                                                                                                               <  2:4025201/150‑1 (MQ=255)
                      gAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAg                                                                                                                               >  1:463869/1‑150 (MQ=255)
                      gAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAg                                                                                                                               <  1:4271728/150‑1 (MQ=255)
                       aaaTCAAAAAACCACTGGCTTTGTTACGCAGATGATTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGAATTCACCCGCCATCCCCAAGCCGAAGACCAGACGAGCGATAAACATGGAGATGTATCCTGGCGCAAAGACGCAgg                                                                                                                              >  1:2562327/1‑150 (MQ=255)
                        aaTCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGc                                                                                                                             >  2:2979443/1‑150 (MQ=255)
                         aTCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGcc                                                                                                                            <  1:905997/150‑1 (MQ=255)
                          tCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAGGCCa                                                                                                                           >  2:1900257/1‑150 (MQ=255)
                           cAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAg                                                                                                                          >  2:62378/1‑150 (MQ=255)
                                aaCCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTcc                                                                                                                     >  2:3647071/1‑150 (MQ=255)
                                 aCCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTccc                                                                                                                    <  2:160106/150‑1 (MQ=255)
                                       ggCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCga                                                                                                              >  2:2066381/1‑150 (MQ=255)
                                         cTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCgaga                                                                                                            >  2:2121935/1‑150 (MQ=255)
                                         cTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCgaga                                                                                                            >  2:4098518/1‑150 (MQ=255)
                                          tttGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCgagaa                                                                                                           >  1:2408945/1‑150 (MQ=255)
                                          tttGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCgagaa                                                                                                           >  1:2193840/1‑150 (MQ=255)
                                                aCGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAAc                                                                                                     <  1:2775924/150‑1 (MQ=255)
                                                  gCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGa                                                                                                   >  1:2981239/1‑150 (MQ=255)
                                                   cAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGAt                                                                                                  <  1:1234937/150‑1 (MQ=255)
                                                      aTGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCt                                                                                               >  2:226230/1‑150 (MQ=255)
                                                        gTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTgg                                                                                             <  2:830025/150‑1 (MQ=255)
                                                         ttttGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGt                                                                                            <  2:1767979/150‑1 (MQ=255)
                                                          tttGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTg                                                                                           <  2:1490388/150‑1 (MQ=255)
                                                            tGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGAc                                                                                         <  1:893096/150‑1 (MQ=255)
                                                             ggCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGAcc                                                                                        <  1:1920455/150‑1 (MQ=255)
                                                               ccAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCAt                                                                                      >  2:5171499/1‑150 (MQ=255)
                                                                 aGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTg                                                                                    >  2:1947479/1‑150 (MQ=255)
                                                                    tttCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCa                                                                                 <  2:4739298/150‑1 (MQ=255)
                                                                     ttCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAg                                                                                >  2:3266174/1‑150 (MQ=255)
                                                                        aaTGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAgacg                                                                             <  2:2895/150‑1 (MQ=255)
                                                                         aTGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAgacga                                                                            <  1:4570374/150‑1 (MQ=255)
                                                                         aTGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAgacga                                                                            <  1:3485117/150‑1 (MQ=255)
                                                                         aTGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAgacga                                                                            >  2:2699501/1‑150 (MQ=255)
                                                                         aTGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAgacga                                                                            >  2:29770/1‑150 (MQ=255)
                                                                         aTGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAgacga                                                                            >  2:4574902/1‑150 (MQ=255)
                                                                         aTGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAgacga                                                                            >  2:104227/1‑150 (MQ=255)
                                                                           gACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGacgacg                                                                          <  1:5171499/150‑1 (MQ=255)
                                                                              aTAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGccc                                                                       >  2:4199657/1‑150 (MQ=255)
                                                                                 ggTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAATGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTa                                                                    <  2:4803829/150‑1 (MQ=255)
                                                                                       gcTGGAACCGTATTCACCCGCCATCCCCATGCCGCTGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTc                                                              <  2:3162564/150‑1 (MQ=255)
                                                                                       gcTGGAACCGTATTCACCCGCCAGCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTc                                                              <  1:5059999/150‑1 (MQ=255)
                                                                                         tGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCAc                                                            <  2:5019221/150‑1 (MQ=255)
                                                                                           gAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCAccc                                                          <  1:83291/150‑1 (MQ=255)
                                                                                             aCCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCAt                                                        <  2:2627008/150‑1 (MQ=255)
                                                                                             aCCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCAt                                                        <  1:3534700/150‑1 (MQ=255)
                                                                                             aCCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCAt                                                        >  1:2933625/1‑150 (MQ=255)
                                                                                               cGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCCTAg                                                      >  2:2425630/1‑150 (MQ=255)
                                                                                                  aTTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAgcgc                                                   >  2:4959207/1‑150 (MQ=255)
                                                                                                   ttCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGcc                                                  <  1:807042/150‑1 (MQ=255)
                                                                                                   ttCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGcc                                                  <  1:4995382/150‑1 (MQ=255)
                                                                                                   ttCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGcc                                                  <  2:754513/150‑1 (MQ=255)
                                                                                                     cACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGa                                                <  1:3696113/150‑1 (MQ=255)
                                                                                                      aCCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAg                                               >  2:4158245/1‑150 (MQ=255)
                                                                                                       cncGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGc                                              >  1:1945155/1‑150 (MQ=255)
                                                                                                             aTCCCCATGCCGATTACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGTTGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAgg                                        >  2:1694095/1‑150 (MQ=255)
                                                                                                             aTCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAgg                                        >  2:4400559/1‑150 (MQ=255)
                                                                                                                cccATGCCGATGACCAGACGATCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGgccg                                     >  1:38568/1‑150 (MQ=255)
                                                                                                                cccATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGgccg                                     >  1:38620/1‑150 (MQ=255)
                                                                                                                  cATGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGgccgcc                                   <  1:408104/150‑1 (MQ=255)
                                                                                                                   aTGCCGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGccgccg                                  <  2:3862848/150‑1 (MQ=255)
                                                                                                                      ccGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAAc                               <  1:3244746/150‑1 (MQ=255)
                                                                                                                      ccGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAAc                               >  2:3789863/1‑150 (MQ=255)
                                                                                                                      ccGATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAAc                               <  2:3157336/150‑1 (MQ=255)
                                                                                                                        gATGACCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCa                             >  1:3311677/1‑150 (MQ=255)
                                                                                                                            aCCAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCga                         >  1:105689/1‑150 (MQ=255)
                                                                                                                             ccAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCgag                        <  1:4438853/150‑1 (MQ=255)
                                                                                                                             ccAGACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATACTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCgag                        <  2:1284703/150‑1 (MQ=255)
                                                                                                                                gACGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCgagaga                     >  1:1785996/1‑150 (MQ=255)
                                                                                                                                gACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATa                                                       <  2:548209/116‑1 (MQ=255)
                                                                                                                                gACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATa                                                       >  1:548209/1‑116 (MQ=255)
                                                                                                                                gACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCgagaga                     <  2:2780169/150‑1 (MQ=255)
                                                                                                                                gACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCgagaga                     <  1:4199657/150‑1 (MQ=255)
                                                                                                                                  cGAGCGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATa                   <  1:2424197/150‑1 (MQ=255)
                                                                                                                                      cGATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAgg               >  1:1407969/1‑150 (MQ=255)
                                                                                                                                       gATAAACATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGc              >  1:1073021/1‑150 (MQ=255)
                                                                                                                                          aaaCATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGc           >  2:635645/1‑150 (MQ=255)
                                                                                                                                           aaCATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCa          <  1:4382904/150‑1 (MQ=255)
                                                                                                                                             cATGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCaga        <  1:1501737/150‑1 (MQ=255)
                                                                                                                                             cATGGTGATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCaga        >  1:4510195/1‑150 (MQ=255)
                                                                                                                                              aTGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCagag       <  2:1207674/150‑1 (MQ=255)
                                                                                                                                              aTGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCagag       <  1:613522/150‑1 (MQ=255)
                                                                                                                                              aTGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCagag       <  1:1374238/150‑1 (MQ=255)
                                                                                                                                               tGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCAgaga      <  1:2948990/150‑1 (MQ=255)
                                                                                                                                               tGGTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCAGTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCAgaga      <  1:4131075/150‑1 (MQ=255)
                                                                                                                                                ggTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCAGAGAt     <  1:226321/150‑1 (MQ=255)
                                                                                                                                                ggTTATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCAGAGAt     >  1:2183499/1‑150 (MQ=255)
                                                                                                                                                 gTGATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCAGAGATc    <  2:1411271/150‑1 (MQ=255)
                                                                                                                                                   gATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCAGAGATCAg  <  2:5004114/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCAGAGATCAG  >  CP009273/3365460‑3365756

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: ‑