Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,696,612 A→T 93.9% E41D (GAA→GAT add → adenosine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,696,6120AT93.9% 941.4 / 39.7 278E41D (GAA→GATaddadenosine deaminase
Reads supporting (aligned to +/- strand):  ref base A (9/8);  new base T (148/112);  total (157/121)
Fisher's exact test for biased strand distribution p-value = 8.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-01

Only 100 of 278 total aligned reads displayed.

ATATCCAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCA  >  CP009273/1696467‑1696758
                                                                                                                                                 |                                                                                                                                                  
atatCCAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATacac                                                                                                                                                >  1:5001713/1‑150 (MQ=255)
 tatCCAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACt                                                                                                                                               >  2:1767770/1‑150 (MQ=255)
  atCCAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTg                                                                                                                                              >  2:3773617/1‑150 (MQ=255)
  atCCAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTg                                                                                                                                              >  1:1467367/1‑150 (MQ=255)
      aGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTc                                                                                                                                          >  2:958704/1‑150 (MQ=255)
      aGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTc                                                                                                                                          >  2:4576481/1‑150 (MQ=255)
        aaaaaGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTccc                                                                                                                                        >  2:3865062/1‑150 (MQ=255)
         aaaaGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTcccc                                                                                                                                       <  2:2696597/150‑1 (MQ=255)
            agagTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACg                                                                                                                                    >  2:3100136/1‑150 (MQ=255)
              agTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGtt                                                                                                                                  >  1:192066/1‑150 (MQ=255)
              agTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGtt                                                                                                                                  >  1:912792/1‑150 (MQ=255)
                  cGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAgg                                                                                                                              >  2:4320818/1‑150 (MQ=255)
                   gACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGt                                                                                                                             >  1:4367670/1‑150 (MQ=255)
                     ccATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCa                                                                                                                           >  2:3132407/1‑150 (MQ=255)
                     ccATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCa                                                                                                                           <  1:4085022/150‑1 (MQ=255)
                       aTGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCTCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCAtt                                                                                                                         <  1:523975/150‑1 (MQ=255)
                       aTGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCAtt                                                                                                                         <  1:4030193/150‑1 (MQ=255)
                       aTGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCAtt                                                                                                                         >  1:3312783/1‑150 (MQ=255)
                         gattgatACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGc                                                                                                                       >  1:4456080/1‑150 (MQ=255)
                             gatACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAAc                                                                                                                   >  1:2760295/1‑150 (MQ=255)
                             gatACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAAc                                                                                                                   >  2:5010958/1‑150 (MQ=255)
                                 ccaccCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGaaa                                                                                                               >  2:3706821/1‑149 (MQ=255)
                                   accCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAAccc                                                                                                             <  1:2531503/150‑1 (MQ=255)
                                   accCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAAccc                                                                                                             <  1:2531362/150‑1 (MQ=255)
                                     ccTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGa                                                                                                           <  2:785774/150‑1 (MQ=255)
                                        gCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCt                                                                                                        <  2:4797140/150‑1 (MQ=255)
                                        gCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCt                                                                                                        >  1:4331307/1‑150 (MQ=255)
                                        gCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCt                                                                                                        <  2:2730535/150‑1 (MQ=255)
                                         ccATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTg                                                                                                       <  1:3492777/150‑1 (MQ=255)
                                         ccATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTg                                                                                                       <  1:2916555/150‑1 (MQ=255)
                                          cATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTgg                                                                                                      <  2:1697345/150‑1 (MQ=255)
                                           aTTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGt                                                                                                     >  2:2007537/1‑150 (MQ=255)
                                           aTTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGt                                                                                                     <  2:1640555/150‑1 (MQ=255)
                                             tAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGa                                                                                                   >  1:2854409/1‑150 (MQ=255)
                                                  gATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCttt                                                                                              <  1:3726044/150‑1 (MQ=255)
                                                     atCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTg                                                                                           >  1:1678482/1‑150 (MQ=255)
                                                      tCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGa                                                                                          <  2:2282486/150‑1 (MQ=255)
                                                        cATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGAcc                                                                                        >  1:2974826/1‑150 (MQ=255)
                                                          tCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCaa                                                                                      <  2:2570550/150‑1 (MQ=255)
                                                           cGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCaaa                                                                                     >  1:4623434/1‑150 (MQ=255)
                                                               aCCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTg                                                                                 <  1:3231828/150‑1 (MQ=255)
                                                                ccTTGATGGCAACCTTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGa                                                                                <  2:4398982/150‑1 (MQ=255)
                                                                  ttGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACt                                                                              <  1:3611499/150‑1 (MQ=255)
                                                                     aTGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTggg                                                                           <  2:3348364/150‑1 (MQ=255)
                                                                          aaCATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGtt                                                                      <  1:4890505/150‑1 (MQ=255)
                                                                          aaCATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGtt                                                                      <  1:1305047/150‑1 (MQ=255)
                                                                          aaCATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGtt                                                                      >  1:3340513/1‑150 (MQ=255)
                                                                           aCATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTa                                                                     >  1:2875484/1‑150 (MQ=255)
                                                                                 gTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTtc                                                               >  2:1814068/1‑150 (MQ=255)
                                                                                 gTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTtc                                                               >  1:1761594/1‑150 (MQ=255)
                                                                                 gTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTtc                                                               >  2:1217819/1‑150 (MQ=255)
                                                                                    cccAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCg                                                            <  2:2772641/150‑1 (MQ=255)
                                                                                      cAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGcc                                                          >  1:2633039/1‑150 (MQ=255)
                                                                                         aCCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCctct                                                       <  2:2975137/150‑1 (MQ=255)
                                                                                         aCCATTCTTGAACGTGGCCTCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCTGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAAATTGACTGGTGCGTTAATGTTCTAGCctct                                                       <  2:4886280/150‑1 (MQ=255)
                                                                                          ccATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCctctc                                                      >  1:3500577/1‑150 (MQ=255)
                                                                                           cATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCctctct                                                     >  2:1173985/1‑150 (MQ=255)
                                                                                           cATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCctctct                                                     <  2:1043361/150‑1 (MQ=255)
                                                                                           cATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCctctct                                                     <  1:4918407/150‑1 (MQ=255)
                                                                                            aTTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCtt                                                    <  2:4136901/150‑1 (MQ=255)
                                                                                             ttCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTg                                                   >  1:702378/1‑150 (MQ=255)
                                                                                              tCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGa                                                  >  1:4696981/1‑150 (MQ=255)
                                                                                               cTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGAt                                                 <  1:213056/150‑1 (MQ=255)
                                                                                               cTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGAt                                                 <  1:284849/150‑1 (MQ=255)
                                                                                                ttGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATg                                                <  2:531617/150‑1 (MQ=255)
                                                                                                  gAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGcc                                              <  2:4208550/150‑1 (MQ=255)
                                                                                                  gAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGcc                                              >  2:275136/1‑150 (MQ=255)
                                                                                                   aaCTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCt                                             >  2:2103260/1‑150 (MQ=255)
                                                                                                    aCTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTg                                            >  1:116046/1‑150 (MQ=255)
                                                                                                    aCTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTg                                            <  2:4367670/150‑1 (MQ=255)
                                                                                                      ttGGCCGCCAGTATAGTATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTc                                          <  2:3647389/150‑1 (MQ=255)
                                                                                                         gccgccAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTcgcc                                       <  1:4283316/150‑1 (MQ=255)
                                                                                                            gccAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCcgcg                                    >  2:1551022/1‑150 (MQ=255)
                                                                                                              cAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTg                                  >  1:2390284/1‑150 (MQ=255)
                                                                                                              cAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTg                                  <  2:4563723/150‑1 (MQ=255)
                                                                                                               aGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTgg                                 >  1:64161/1‑150 (MQ=255)
                                                                                                               aGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTgg                                 >  1:1440448/1‑150 (MQ=255)
                                                                                                                  ataataTCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCAt                              <  2:2362156/150‑1 (MQ=255)
                                                                                                                   taataTCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCAtt                             <  2:4475728/150‑1 (MQ=255)
                                                                                                                   taataTCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCAtt                             <  2:1426755/150‑1 (MQ=255)
                                                                                                                       atCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCCTTTGaa                         >  1:4505592/1‑150 (MQ=255)
                                                                                                                          tcGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCCTTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCAGTTGAAAAc                      >  2:4986952/1‑150 (MQ=255)
                                                                                                                          tcGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAAc                      >  2:2170955/1‑150 (MQ=255)
                                                                                                                              ttCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTg                  >  2:5134018/1‑150 (MQ=255)
                                                                                                                              ttCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTg                  >  1:224005/1‑150 (MQ=255)
                                                                                                                               tCCTTCACAATCCCTGTATACACTGATTCCCCACGTTCAGGTCATTTCCAACGAACCCGATCTTGTGAGCTTTCTGACCAAACTTGACTGGGGAGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCTCGTGGCATTTGAAAACATTGa                 <  1:3960118/150‑1 (MQ=255)
                                                                                                                                 cTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAg               >  2:4231415/1‑150 (MQ=255)
                                                                                                                                  tGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGa              >  1:1002589/1‑150 (MQ=255)
                                                                                                                                   gCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGAt             >  1:1950712/1‑150 (MQ=255)
                                                                                                                                    cacaATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATg            >  2:4009764/1‑150 (MQ=255)
                                                                                                                                       aaTCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAg         >  1:1285686/1‑150 (MQ=255)
                                                                                                                                       aaTCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATGTGAAAACATTGAAGATGCAg         >  2:3486721/1‑150 (MQ=255)
                                                                                                                                         tCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGcc       <  2:3301399/150‑1 (MQ=255)
                                                                                                                                         tCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGcc       >  2:4936148/1‑150 (MQ=255)
                                                                                                                                         tCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGcc       <  2:1387447/150‑1 (MQ=255)
                                                                                                                                           ccTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCg     <  2:2218469/150‑1 (MQ=255)
                                                                                                                                           ccTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGATCTCGCCTCTCTTGATGACTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCg     >  1:3424281/1‑150 (MQ=255)
                                                                                                                                             tGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTc   >  2:2298883/1‑150 (MQ=255)
                                                                                                                                             tGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTc   >  1:3875626/1‑150 (MQ=255)
                                                                                                                                              ggATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCa  >  1:3099821/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
ATATCCAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCA  >  CP009273/1696467‑1696758

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: