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.8% 840.4 / 0.3 243E41D (GAA→GATaddadenosine deaminase
Reads supporting (aligned to +/- strand):  ref base A (0/3);  new base T (131/109);  total (131/112)
Fisher's exact test for biased strand distribution p-value = 9.65e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 243 total aligned reads displayed.

ATATCCAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCA  >  CP009273/1696467‑1696758
                                                                                                                                                 |                                                                                                                                                  
atatCCAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATacac                                                                                                                                                >  2:3548829/1‑150 (MQ=255)
   tCCAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGa                                                                                                                                             <  1:235805/150‑1 (MQ=255)
    ccAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGAt                                                                                                                                            >  1:2913748/1‑150 (MQ=255)
        aaaaaGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTccc                                                                                                                                        >  1:1765125/1‑150 (MQ=255)
           aaGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCAc                                                                                                                                     >  1:2780353/1‑150 (MQ=255)
            agagTGCGACCATGATTGATACCACCCTGCCATTAACTTATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACg                                                                                                                                    <  2:1053831/150‑1 (MQ=255)
                 gCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAg                                                                                                                               >  2:3486819/1‑150 (MQ=255)
                       aTGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCAtt                                                                                                                         <  1:3138463/150‑1 (MQ=255)
                          attgatACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGcc                                                                                                                      >  2:788471/1‑150 (MQ=255)
                          attgatACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGcc                                                                                                                      >  1:3800214/1‑150 (MQ=255)
                           ttgatACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCa                                                                                                                     >  2:2839073/1‑150 (MQ=255)
                            tgatACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCaa                                                                                                                    >  1:1883134/1‑150 (MQ=255)
                            tgatACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCaa                                                                                                                    <  2:1866489/150‑1 (MQ=255)
                            tgatACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCaa                                                                                                                    <  2:1866568/150‑1 (MQ=255)
                              atACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACg                                                                                                                  >  1:2985485/1‑150 (MQ=255)
                                   accCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAAccc                                                                                                             <  1:254065/150‑1 (MQ=255)
                                    cccTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCg                                                                                                            >  1:1653137/1‑150 (MQ=255)
                                       tGCCATTAACTGATATCCATCGCCACCTTGATGTCAACATTAGTCCCGAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGTTTCCCCACGTTCAGGTCATTGCCAACGAACCCGATc                                                                                                         <  2:200038/150‑1 (MQ=255)
                                       tGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATc                                                                                                         <  2:199676/150‑1 (MQ=255)
                                        gCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCt                                                                                                        <  2:336085/150‑1 (MQ=255)
                                          cATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTgg                                                                                                      <  2:1061181/150‑1 (MQ=255)
                                           aTTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGt                                                                                                     >  2:3960419/1‑150 (MQ=255)
                                           aTTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGGTCAGGTCATTGCCAACGAACCCGATCTGGt                                                                                                     >  2:3960522/1‑150 (MQ=255)
                                              aaCTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAg                                                                                                  >  2:2126933/1‑150 (MQ=255)
                                              aaCTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAg                                                                                                  >  1:1122572/1‑150 (MQ=255)
                                                 tGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCtt                                                                                               >  2:1659368/1‑150 (MQ=255)
                                                  gATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCttt                                                                                              <  1:4217670/150‑1 (MQ=255)
                                                   atatCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTc                                                                                             >  2:38076/1‑150 (MQ=255)
                                                     atCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTg                                                                                           <  1:1247518/150‑1 (MQ=255)
                                                       ccATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGAc                                                                                         <  1:2227736/150‑1 (MQ=255)
                                                           cGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAAGGAACCCGATCTGGTGAGCTTTCTGACCaaa                                                                                     <  1:1954263/150‑1 (MQ=255)
                                                            gCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAAc                                                                                    >  2:2735308/1‑150 (MQ=255)
                                                              cACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACtt                                                                                  <  1:603310/150‑1 (MQ=255)
                                                              cACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACtt                                                                                  <  2:1600575/150‑1 (MQ=255)
                                                               aCCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTg                                                                                 >  2:523031/1‑150 (MQ=255)
                                                                     aTGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTggg                                                                           >  1:3325942/1‑150 (MQ=255)
                                                                     aTGGCAACATTCGTCCCCAGACCATTATTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTggg                                                                           <  2:2957342/150‑1 (MQ=255)
                                                                       ggCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGc                                                                         <  2:1029884/150‑1 (MQ=255)
                                                                          aaCATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGtt                                                                      >  1:4081403/1‑150 (MQ=255)
                                                                             aTTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTaaa                                                                   >  1:1958573/1‑150 (MQ=255)
                                                                               tCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGt                                                                 <  1:244809/150‑1 (MQ=255)
                                                                               tCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGt                                                                 <  1:244759/150‑1 (MQ=255)
                                                                                cGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCACTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGtt                                                                <  2:3292361/150‑1 (MQ=255)
                                                                                  tCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTtct                                                              >  1:4631969/1‑150 (MQ=255)
                                                                                  tCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTtct                                                              >  2:2292866/1‑150 (MQ=255)
                                                                                  tCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTtct                                                              >  1:396712/1‑150 (MQ=255)
                                                                                   ccccAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTtctc                                                             <  1:341277/150‑1 (MQ=255)
                                                                                    cccAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCg                                                            >  1:796403/1‑150 (MQ=255)
                                                                                      cAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCTAACTTGACTGGGGCGTTAAAGTTCTCGcc                                                          <  2:2694647/150‑1 (MQ=255)
                                                                                      cAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGcc                                                          >  1:3702563/1‑150 (MQ=255)
                                                                                      cAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGcc                                                          <  1:3202943/150‑1 (MQ=255)
                                                                                       aGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCct                                                         >  2:2962689/1‑150 (MQ=255)
                                                                                       aGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCct                                                         >  1:706353/1‑150 (MQ=255)
                                                                                       aGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCct                                                         >  1:3025964/1‑150 (MQ=255)
                                                                                           cATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCctctct                                                     >  2:4090098/1‑150 (MQ=255)
                                                                                           cATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCctctct                                                     <  2:693734/150‑1 (MQ=255)
                                                                                             ttCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTg                                                   >  2:1455923/1‑150 (MQ=255)
                                                                                             ttCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTg                                                   >  1:1027485/1‑150 (MQ=255)
                                                                                             ttCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCttt                                                   <  1:184826/150‑2 (MQ=255)
                                                                                                 tGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGc                                               >  2:377607/1‑150 (MQ=255)
                                                                                                 tGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGc                                               >  2:908272/1‑150 (MQ=255)
                                                                                                  gAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGcc                                              <  2:3594933/150‑1 (MQ=255)
                                                                                                  gAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGcc                                              >  1:536804/1‑150 (MQ=255)
                                                                                                  gAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCCGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGcc                                              >  1:3395435/1‑150 (MQ=255)
                                                                                                  gAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTCGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGcc                                              <  1:2673673/150‑1 (MQ=255)
                                                                                                   aaCTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCt                                             <  1:3283161/150‑1 (MQ=255)
                                                                                                    aCTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTg                                            >  2:2468232/1‑150 (MQ=255)
                                                                                                       tGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTcg                                         <  2:230177/150‑1 (MQ=255)
                                                                                                         gccgccAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTcgcc                                       >  2:2879539/1‑150 (MQ=255)
                                                                                                           cgccAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTcgccgc                                     >  2:2053809/1‑150 (MQ=255)
                                                                                                           cgccAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTcgccgc                                     >  2:2053875/1‑150 (MQ=255)
                                                                                                            gccAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTTATGCCTGTCGCcgcg                                    >  1:2887309/1‑150 (MQ=255)
                                                                                                               aGTATAATATCTCCCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGACTCTCTAGATGCCTGTCGCCGCGTgg                                 >  2:4375212/1‑150 (MQ=255)
                                                                                                                  ataataTCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCAt                              >  1:3620261/1‑150 (MQ=255)
                                                                                                                   taataTCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCAtt                             <  1:2094380/150‑1 (MQ=255)
                                                                                                                   taataTCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCAtt                             <  1:4425135/150‑1 (MQ=255)
                                                                                                                     atatCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTg                           >  2:3079706/1‑150 (MQ=255)
                                                                                                                      tatCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGa                          <  1:3486819/150‑1 (MQ=255)
                                                                                                                       atCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGaa                         >  2:2405299/1‑150 (MQ=255)
                                                                                                                       atCTCGCTTCCTGCACAATCCCTGGATACACTGAGTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGaa                         <  2:1104865/150‑1 (MQ=255)
                                                                                                                        tctcGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGaaa                        >  1:885395/1‑150 (MQ=255)
                                                                                                                          tcGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAAc                      >  1:1215182/1‑150 (MQ=255)
                                                                                                                             cTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACAtt                   >  2:207896/1‑150 (MQ=255)
                                                                                                                             cTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACAtt                   >  2:1687699/1‑150 (MQ=255)
                                                                                                                             cTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCAGTTGAAAACAtt                   >  2:1242979/1‑150 (MQ=255)
                                                                                                                                ccTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGaa                <  1:2848712/150‑1 (MQ=255)
                                                                                                                                ccTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGaa                <  1:1590011/150‑1 (MQ=255)
                                                                                                                                ccTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGaa                <  2:3155602/150‑1 (MQ=255)
                                                                                                                                  tGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGa              <  1:4106061/150‑1 (MQ=255)
                                                                                                                                   gCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGAt             >  1:859209/1‑150 (MQ=255)
                                                                                                                                    cacaATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATg            >  2:2912950/1‑150 (MQ=255)
                                                                                                                                    cacaATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATg            >  2:411996/1‑150 (MQ=255)
                                                                                                                                      caATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCa          <  2:1962109/150‑1 (MQ=255)
                                                                                                                                      caATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCCGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCa          <  2:1765125/150‑1 (MQ=255)
                                                                                                                                       aaTCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAg         <  1:326252/150‑1 (MQ=255)
                                                                                                                                           ccTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCg     >  1:1279489/1‑150 (MQ=255)
                                                                                                                                           ccTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCg     <  2:4174901/150‑1 (MQ=255)
                                                                                                                                           ccTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCg     <  2:4306727/150‑1 (MQ=255)
                                                                                                                                             tGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTc   >  2:2687392/1‑150 (MQ=255)
                                                                                                                                              ggATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCa  <  1:2594700/150‑1 (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: