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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,696,6120AT73.2% 460.8 / 203.9 239E41D (GAA→GAT) addadenosine deaminase
Reads supporting (aligned to +/- strand):  ref base A (37/27);  new base T (100/74);  total (137/102)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 240 total aligned reads displayed.

TTATATCCAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACG  >  CP009273/1696465‑1696760
                                                                                                                                                   |                                                                                                                                                    
ttATATCCAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATac                                                                                                                                                    >  1:79484/1‑150 (MQ=255)
   tatCCAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACt                                                                                                                                                 >  1:709345/1‑150 (MQ=255)
         gAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTcc                                                                                                                                           >  1:3764594/1‑150 (MQ=255)
          aaaaaGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTccc                                                                                                                                          >  2:389319/1‑150 (MQ=255)
          aaaaaGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTccc                                                                                                                                          >  1:3660773/1‑150 (MQ=255)
          aaaaaGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTccc                                                                                                                                          >  2:3682385/1‑150 (MQ=255)
            aaaGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCa                                                                                                                                        >  1:2509607/1‑150 (MQ=255)
              agagTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACg                                                                                                                                      >  2:1634674/1‑150 (MQ=255)
              agagTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACg                                                                                                                                      >  2:2977660/1‑150 (MQ=255)
                  tGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCa                                                                                                                                  >  1:4056635/1‑150 (MQ=255)
                   gCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAg                                                                                                                                 >  2:1361054/1‑150 (MQ=255)
                     gACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGt                                                                                                                               <  1:4619738/150‑1 (MQ=255)
                      aCCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTc                                                                                                                              >  1:2028198/1‑150 (MQ=255)
                        cATGATTGATACCAACCTGACATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCAt                                                                                                                            <  2:2246332/150‑1 (MQ=255)
                          tgattgatACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTg                                                                                                                          <  2:1628634/150‑1 (MQ=255)
                          tgattgatACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTg                                                                                                                          <  1:3586992/150‑1 (MQ=255)
                           gattgatACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGc                                                                                                                         <  1:112615/150‑1 (MQ=255)
                           gattgatACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGc                                                                                                                         >  2:3215093/1‑150 (MQ=255)
                             ttgatACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCa                                                                                                                       >  1:2770598/1‑150 (MQ=255)
                             ttgatACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGCCCGCCAGTATAAAATCTCGCTTCCTGCACAATCCCTGGATAAACTGATTCACCACTATCAGCTCAATGACa                                                                                                                       >  1:2345619/1‑150 (MQ=255)
                              tgatACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCaa                                                                                                                      >  2:2188187/1‑150 (MQ=255)
                                 tACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGa                                                                                                                   >  2:496245/1‑150 (MQ=255)
                                  accaccCTGCCATTAACTGATATCCATCTCCACCTTGATGGCAAAATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCTCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGaa                                                                                                                  >  2:488537/1‑150 (MQ=255)
                                     accCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAAccc                                                                                                               >  1:5115364/1‑150 (MQ=255)
                                     accCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAAccc                                                                                                               <  1:2430570/150‑1 (MQ=255)
                                     accCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAAccc                                                                                                               <  1:478911/150‑1 (MQ=255)
                                         tGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATc                                                                                                           >  2:1020663/1‑150 (MQ=255)
                                          gCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCt                                                                                                          <  1:5148761/150‑1 (MQ=255)
                                           ccATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTg                                                                                                         <  2:3202662/150‑1 (MQ=255)
                                           ccATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTg                                                                                                         <  2:4972604/150‑1 (MQ=255)
                                             aTTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGt                                                                                                       >  2:1938765/1‑150 (MQ=255)
                                               tAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGa                                                                                                     <  1:3810099/150‑1 (MQ=255)
                                                 aCTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGc                                                                                                   >  1:4668854/1‑150 (MQ=255)
                                                     atatCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCTTTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTc                                                                                               <  1:3054624/150‑1 (MQ=255)
                                                       atCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTg                                                                                             >  2:347832/1‑150 (MQ=255)
                                                       atCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTg                                                                                             >  2:5010045/1‑150 (MQ=255)
                                                        tCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGa                                                                                            >  1:1854366/1‑150 (MQ=255)
                                                         ccATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGAc                                                                                           <  2:2798137/150‑1 (MQ=255)
                                                         ccATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGAc                                                                                           >  2:2380580/1‑150 (MQ=255)
                                                            tCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCaa                                                                                        >  2:583935/1‑150 (MQ=255)
                                                             cGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCaaa                                                                                       >  1:3410581/1‑150 (MQ=255)
                                                                  ccTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGa                                                                                  >  2:4133813/1‑150 (MQ=255)
                                                                   cTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGAc                                                                                 <  1:4087446/150‑1 (MQ=255)
                                                                     tGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTg                                                                               <  2:1182520/150‑1 (MQ=255)
                                                                      gATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTgg                                                                              <  1:3112216/150‑1 (MQ=255)
                                                                       aTGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTggg                                                                             >  2:1740639/1‑150 (MQ=255)
                                                                       aTGGCAACAGTCGTCCCCAGCCCTTTCGTGAACTTGGCCGCCAGTATAATCTCTCGCTTCCTGCACAATCCCTGGAGCCACTGATTCCCAACGGTCAGGTCATTGCCAACGAACCCGATCGGGTGAGCTTTCTGACCAAACTTGACTggg                                                                             <  1:3185584/150‑1 (MQ=255)
                                                                         ggCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGc                                                                           <  1:1584003/150‑1 (MQ=255)
                                                                               aTTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTaaa                                                                     >  1:11851/1‑150 (MQ=255)
                                                                                 tCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGt                                                                   <  1:5030820/150‑1 (MQ=255)
                                                                                  cGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGtt                                                                  <  2:3442164/150‑1 (MQ=255)
                                                                                     ccccAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTtctc                                                               >  2:1685850/1‑150 (MQ=255)
                                                                                     ccccAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTTAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTtctc                                                               >  2:1806012/1‑150 (MQ=255)
                                                                                        cAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGcc                                                            >  2:1449702/1‑150 (MQ=255)
                                                                                        cAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGcc                                                            >  2:954228/1‑150 (MQ=255)
                                                                                           aCCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCctct                                                         >  1:3458771/1‑150 (MQ=255)
                                                                                               ttCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTg                                                     <  1:2248706/150‑1 (MQ=255)
                                                                                                tCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGa                                                    >  2:4359377/1‑150 (MQ=255)
                                                                                                 cTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCANAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGAt                                                   >  1:4790170/1‑150 (MQ=255)
                                                                                                  ttGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATg                                                  >  2:4254942/1‑150 (MQ=255)
                                                                                                     aaCTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCt                                               >  2:4523520/1‑150 (MQ=255)
                                                                                                     aaCTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCt                                               >  2:4523958/1‑150 (MQ=255)
                                                                                                     aaCTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCt                                               >  2:5122611/1‑150 (MQ=255)
                                                                                                     aaCTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCt                                               >  2:5123955/1‑150 (MQ=255)
                                                                                                     aaCTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCt                                               >  1:981286/1‑150 (MQ=255)
                                                                                                       cTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGt                                             <  1:4017707/150‑1 (MQ=255)
                                                                                                        ttGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTc                                            <  1:1013159/150‑1 (MQ=255)
                                                                                                        ttGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTc                                            >  1:1440320/1‑150 (MQ=255)
                                                                                                        ttGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTc                                            <  2:3301133/150‑1 (MQ=255)
                                                                                                        ttGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTc                                            <  2:695514/150‑1 (MQ=255)
                                                                                                         tGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTcg                                           <  2:1397925/150‑1 (MQ=255)
                                                                                                         tGGCCGCCAGTATAATATCTCGCCTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTcg                                           <  1:205725/150‑1 (MQ=255)
                                                                                                           gccgccAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTcgcc                                         <  1:1701050/150‑1 (MQ=255)
                                                                                                              gccAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCcgcg                                      <  1:4895716/150‑1 (MQ=255)
                                                                                                               ccAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGt                                     >  2:1550361/1‑150 (MQ=255)
                                                                                                                  gTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGc                                  >  1:1722026/1‑150 (MQ=255)
                                                                                                                   tataATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTAGAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCa                                 <  2:4725158/150‑1 (MQ=255)
                                                                                                                    ataataTCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCAt                                >  2:1159171/1‑150 (MQ=255)
                                                                                                                     taataTCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCAtt                               <  1:3624206/150‑1 (MQ=255)
                                                                                                                     taataTCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCAtt                               >  1:3472140/1‑150 (MQ=255)
                                                                                                                       atatCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTg                             >  2:2033409/1‑150 (MQ=255)
                                                                                                                       atatCTCGCTTCCTGCACAATCCCTGGATACACTGAGTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTg                             <  1:4721153/150‑1 (MQ=255)
                                                                                                                       atatCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTg                             >  2:754122/1‑150 (MQ=255)
                                                                                                                          tctcGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGaaa                          >  2:2720512/1‑150 (MQ=255)
                                                                                                                          tctcGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGaaa                          >  2:1564066/1‑150 (MQ=255)
                                                                                                                                  ccTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGaa                  <  2:1485117/150‑1 (MQ=255)
                                                                                                                                  ccTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGaa                  <  1:1316830/150‑1 (MQ=255)
                                                                                                                                  ccTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGaa                  >  2:2549249/1‑150 (MQ=255)
                                                                                                                                   cTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAg                 >  2:3098829/1‑150 (MQ=255)
                                                                                                                                   cTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAg                 <  1:1043201/150‑1 (MQ=255)
                                                                                                                                    tGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGa                >  2:612551/1‑150 (MQ=255)
                                                                                                                                     gCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGAt               >  1:5095230/1‑150 (MQ=255)
                                                                                                                                     gCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGAt               <  1:3585282/150‑1 (MQ=255)
                                                                                                                                     gCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGAt               >  1:2641978/1‑150 (MQ=255)
                                                                                                                                      cacaATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATg              >  1:3997438/1‑150 (MQ=255)
                                                                                                                                         aaTCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAg           >  1:3230271/1‑150 (MQ=255)
                                                                                                                                            cccTGGATACACTGATTCCCCCCGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGccc        <  2:2458618/150‑1 (MQ=255)
                                                                                                                                              cTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTCATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGt      >  1:877424/1‑150 (MQ=255)
                                                                                                                                                 gAAACTCTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGACCGTCAc   <  2:3979722/150‑1 (MQ=255)
                                                                                                                                                  aaaCACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACg  >  2:119453/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTATATCCAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACG  >  CP009273/1696465‑1696760

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