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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,696,6120AT98.3% 1017.1 / 4.7 291E41D (GAA→GAT) addadenosine deaminase
Reads supporting (aligned to +/- strand):  ref base A (2/3);  new base T (163/123);  total (165/126)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 291 total aligned reads displayed.

TTATATCCAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCA  >  CP009273/1696465‑1696758
                                                                                                                                                   |                                                                                                                                                  
ttATATCCAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATac                                                                                                                                                  <  1:1103662/150‑1 (MQ=255)
 tataTCCAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATaca                                                                                                                                                 <  2:4370640/150‑1 (MQ=255)
  atatCCAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATacac                                                                                                                                                >  2:176825/1‑150 (MQ=255)
  atatCCAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATacac                                                                                                                                                >  2:112354/1‑150 (MQ=255)
   tatCCAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACt                                                                                                                                               >  2:171155/1‑150 (MQ=255)
    atCCAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTg                                                                                                                                              <  1:4250006/150‑1 (MQ=255)
        aGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCGGGATACACTGATTc                                                                                                                                          >  2:315458/1‑150 (MQ=255)
            aaaGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCa                                                                                                                                      <  1:839041/150‑1 (MQ=255)
              agagTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACg                                                                                                                                    >  1:535516/1‑150 (MQ=255)
               gagTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGt                                                                                                                                   <  2:1428450/150‑1 (MQ=255)
               gagTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGt                                                                                                                                   >  2:4207533/1‑150 (MQ=255)
                agTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGtt                                                                                                                                  >  2:4007763/1‑150 (MQ=255)
                       ccATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCa                                                                                                                           >  1:4306685/1‑150 (MQ=255)
                            attgatACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGcc                                                                                                                      >  1:2655150/1‑150 (MQ=255)
                            attgatACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGcc                                                                                                                      >  1:1947165/1‑150 (MQ=255)
                             ttgatACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCa                                                                                                                     <  1:2984948/150‑1 (MQ=255)
                              tgatACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCaa                                                                                                                    >  1:4369486/1‑150 (MQ=255)
                               gatACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAAc                                                                                                                   <  2:2856466/150‑1 (MQ=255)
                                atACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACg                                                                                                                  >  2:2937771/1‑150 (MQ=255)
                                atACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACg                                                                                                                  >  2:3564359/1‑150 (MQ=255)
                                 tACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGa                                                                                                                 >  2:4201400/1‑150 (MQ=255)
                                  accaccCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGaa                                                                                                                >  1:621963/1‑150 (MQ=255)
                                    caccCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAAcc                                                                                                              >  1:419483/1‑150 (MQ=255)
                                     accCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAAccc                                                                                                             <  2:4599657/150‑1 (MQ=255)
                                        cTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGAt                                                                                                          <  2:412958/150‑1 (MQ=255)
                                         tGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATc                                                                                                         >  2:1422171/1‑150 (MQ=255)
                                           ccATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTg                                                                                                       <  2:542336/150‑1 (MQ=255)
                                                aaCTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAg                                                                                                  <  1:3691196/150‑1 (MQ=255)
                                                aaCTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAg                                                                                                  <  2:3790354/150‑1 (MQ=255)
                                                aaCTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAg                                                                                                  <  2:192005/150‑1 (MQ=255)
                                                    gATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCttt                                                                                              <  1:2326830/150‑1 (MQ=255)
                                                    gATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCttt                                                                                              >  1:375858/1‑150 (MQ=255)
                                                    gATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCttt                                                                                              >  2:2567912/1‑150 (MQ=255)
                                                      tatCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCt                                                                                            <  2:87449/150‑1 (MQ=255)
                                                        tCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGa                                                                                          >  1:966049/1‑150 (MQ=255)
                                                        tCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGa                                                                                          <  1:2574401/150‑1 (MQ=255)
                                                          cATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGAcc                                                                                        >  2:2740456/1‑150 (MQ=255)
                                                           aTCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCa                                                                                       <  2:812030/150‑1 (MQ=255)
                                                            tCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCaa                                                                                      >  1:3625868/1‑150 (MQ=255)
                                                             cGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCaaa                                                                                     >  2:792673/1‑150 (MQ=255)
                                                               ccACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACt                                                                                   >  2:1771364/1‑150 (MQ=255)
                                                               ccACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACt                                                                                   >  2:4497997/1‑150 (MQ=255)
                                                                cACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACtt                                                                                  <  2:3676349/150‑1 (MQ=255)
                                                                 aCCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTg                                                                                 >  2:717150/1‑150 (MQ=255)
                                                                  ccTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGa                                                                                >  1:1699020/1‑150 (MQ=255)
                                                                   cTTGATGGCAACATTAGTCACAAGACAATTCTTGAAATTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCATGGATACACTGATTCACCACGTTCAGGACATTGCCAACGAACACGATCTGGTGAGCTTTCTGACCAAACTTTAc                                                                               >  2:963265/1‑150 (MQ=255)
                                                                    ttGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACt                                                                              <  1:1295121/150‑1 (MQ=255)
                                                                    ttGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACt                                                                              <  1:793307/150‑1 (MQ=255)
                                                                        tGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTgggg                                                                          >  2:3091750/1‑150 (MQ=255)
                                                                          gCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCg                                                                        <  1:2500018/150‑1 (MQ=255)
                                                                           cAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGt                                                                       <  2:2153238/150‑1 (MQ=255)
                                                                           cAAAATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCTCTTCCTGCACAATCCCTGGATACAATGATTCCCCACGTTAAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGt                                                                       <  1:877811/150‑1 (MQ=255)
                                                                            aaCATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGtt                                                                      >  2:2368180/1‑150 (MQ=255)
                                                                             aCATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTa                                                                     >  1:322701/1‑150 (MQ=255)
                                                                             aCATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCACACGTTCAGGTCATTGCCAACGAACCCGATCTGTTGAGCTTTCTGACCAAACTTGACTGGGGCGTTa                                                                     >  2:1461242/1‑150 (MQ=255)
                                                                               aTTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTaaa                                                                   >  2:1060692/1‑150 (MQ=255)
                                                                               aTTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTaaa                                                                   >  1:1588358/1‑150 (MQ=255)
                                                                                ttCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAg                                                                  >  2:2907296/1‑150 (MQ=255)
                                                                                     ccccAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTtctc                                                             >  2:2832766/1‑150 (MQ=255)
                                                                                     ccccAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTtctc                                                             <  2:1313801/150‑1 (MQ=255)
                                                                                        cAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGcc                                                          <  1:4205413/150‑1 (MQ=255)
                                                                                         aGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCct                                                         >  2:1677275/1‑150 (MQ=255)
                                                                                         aGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCct                                                         <  1:3294570/150‑1 (MQ=255)
                                                                                          gACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCctc                                                        <  2:3909237/150‑1 (MQ=255)
                                                                                           aCCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCctct                                                       <  1:2228461/150‑1 (MQ=255)
                                                                                            ccATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCctctc                                                      >  1:2724/1‑150 (MQ=255)
                                                                                             cATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCctctct                                                     <  1:3510138/150‑1 (MQ=255)
                                                                                                 cTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGAt                                                 <  1:2277016/150‑1 (MQ=255)
                                                                                                 cTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGAt                                                 >  2:425637/1‑150 (MQ=255)
                                                                                                       cTTGGCCGCCAGTACAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGt                                           >  2:1896226/1‑150 (MQ=255)
                                                                                                        ttGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTc                                          <  2:3489212/150‑1 (MQ=255)
                                                                                                          ggCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTcgc                                        >  2:402666/1‑150 (MQ=255)
                                                                                                           gccgccAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTcgcc                                       <  2:522288/150‑1 (MQ=255)
                                                                                                           gccgccAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTGAAAGTTCTCGCCTCTCTTGATGCCTGTcgcc                                       >  1:2039460/1‑150 (MQ=255)
                                                                                                               ccAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGt                                   >  2:357665/1‑150 (MQ=255)
                                                                                                                cAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTg                                  <  1:1032240/150‑1 (MQ=255)
                                                                                                                cAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTg                                  >  1:1005246/1‑150 (MQ=255)
                                                                                                                cAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTg                                  <  2:106646/150‑1 (MQ=255)
                                                                                                                    ataataTCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCAt                              >  1:4006614/1‑150 (MQ=255)
                                                                                                                    ataataTCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCATGGCAt                              >  1:3396523/1‑150 (MQ=255)
                                                                                                                    ataataTCTCGCGTCCTGCACAATCCCTGGATACACTGATTCCCCAGGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCAt                              <  2:844673/150‑1 (MQ=255)
                                                                                                                             cGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACa                     >  1:4604168/1‑150 (MQ=255)
                                                                                                                              gCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACAt                    <  1:3388394/150‑1 (MQ=255)
                                                                                                                                ttCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTg                  >  1:4585332/1‑150 (MQ=255)
                                                                                                                                ttCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTg                  <  1:1880373/150‑1 (MQ=255)
                                                                                                                                  ccTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGaa                >  2:2648699/1‑150 (MQ=255)
                                                                                                                                    tGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGa              >  2:3226475/1‑150 (MQ=255)
                                                                                                                                    tGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGa              >  2:2528051/1‑150 (MQ=255)
                                                                                                                                    tGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGa              >  1:512241/1‑150 (MQ=255)
                                                                                                                                        caATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCa          <  2:4495891/150‑1 (MQ=255)
                                                                                                                                        caATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCa          <  1:3195715/150‑1 (MQ=255)
                                                                                                                                           tCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGcc       <  2:3629786/150‑1 (MQ=255)
                                                                                                                                            cccTGGATTCACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGccc      <  1:2392587/150‑1 (MQ=255)
                                                                                                                                             ccTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCg     >  1:1201974/1‑150 (MQ=255)
                                                                                                                                               tGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTc   <  1:2120017/150‑1 (MQ=255)
                                                                                                                                                ggATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCa  >  1:3961032/1‑150 (MQ=255)
                                                                                                                                                ggATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCa  >  1:124858/1‑150 (MQ=255)
                                                                                                                                                ggATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCa  <  2:1899578/150‑1 (MQ=255)
                                                                                                                                                ggATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCa  >  2:3290230/1‑150 (MQ=255)
                                                                                                                                                ggATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCa  >  2:972453/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTATATCCAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCA  >  CP009273/1696465‑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: ‑