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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,696,6120AT94.9% 1096.3 / 35.6 316E41D (GAA→GAT) addadenosine deaminase
Reads supporting (aligned to +/- strand):  ref base A (8/8);  new base T (177/123);  total (185/131)
Fisher's exact test for biased strand distribution p-value = 6.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

Only 100 of 318 total aligned reads displayed.

TTTTATATCCAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCA  >  CP009273/1696463‑1696758
                                                                                                                                                     |                                                                                                                                                  
ttttATATCCAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGaa                                                                                                                                                    <  2:678113/150‑1 (MQ=255)
  ttATATCCAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATac                                                                                                                                                  >  1:1220633/1‑150 (MQ=255)
   tataTCCAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATaca                                                                                                                                                 >  1:4181359/1‑150 (MQ=255)
    atatCCAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATCTCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATacac                                                                                                                                                <  1:3607438/150‑1 (MQ=255)
     tatCCAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACt                                                                                                                                               <  2:2040673/150‑1 (MQ=255)
        ccAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGAt                                                                                                                                            <  2:4380200/150‑1 (MQ=255)
         cAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGAtt                                                                                                                                           >  2:608408/1‑150 (MQ=255)
          aGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTc                                                                                                                                          >  2:2984188/1‑150 (MQ=255)
          aGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTc                                                                                                                                          <  2:2428614/150‑1 (MQ=255)
           gAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTcc                                                                                                                                         >  2:1630422/1‑150 (MQ=255)
               aaGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCAc                                                                                                                                     >  1:1463644/1‑150 (MQ=255)
                  agTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGtt                                                                                                                                  <  1:1424116/150‑1 (MQ=255)
                  agTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGtt                                                                                                                                  >  2:420946/1‑150 (MQ=255)
                    tGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCa                                                                                                                                >  1:3345670/1‑150 (MQ=255)
                    tGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCa                                                                                                                                <  2:3777105/150‑1 (MQ=255)
                       gACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGt                                                                                                                             >  1:1671979/1‑150 (MQ=255)
                       gACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGt                                                                                                                             >  1:4075352/1‑150 (MQ=255)
                        aCCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTc                                                                                                                            <  1:2636981/150‑1 (MQ=255)
                         ccATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCa                                                                                                                           >  1:4657455/1‑150 (MQ=255)
                          cATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGAGGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCAt                                                                                                                          <  2:3489394/150‑1 (MQ=255)
                           aTGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCAtt                                                                                                                         >  2:3349388/1‑150 (MQ=255)
                            tgattgatACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTg                                                                                                                        >  1:1743614/1‑150 (MQ=255)
                             gattgatACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGc                                                                                                                       >  2:545783/1‑150 (MQ=255)
                                 gatACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAAc                                                                                                                   <  1:1927468/150‑1 (MQ=255)
                                 gatACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAAc                                                                                                                   <  2:2370217/150‑1 (MQ=255)
                                    accaccCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGaa                                                                                                                >  2:4189556/1‑150 (MQ=255)
                                       accCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAAccc                                                                                                             <  2:2328738/150‑1 (MQ=255)
                                         ccTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGa                                                                                                           >  1:3400775/1‑150 (MQ=255)
                                           tGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATc                                                                                                         >  1:2389445/1‑150 (MQ=255)
                                            gCCATTGACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCt                                                                                                        <  2:1876476/150‑1 (MQ=255)
                                            gCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCt                                                                                                        <  1:2492164/150‑1 (MQ=255)
                                              cATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTgg                                                                                                      >  2:1716198/1‑150 (MQ=255)
                                              cATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTgg                                                                                                      >  1:3276596/1‑150 (MQ=255)
                                               aTTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGt                                                                                                     >  1:2274325/1‑150 (MQ=255)
                                               aTTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGt                                                                                                     <  2:383601/150‑1 (MQ=255)
                                                ttAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTg                                                                                                    <  2:88058/150‑1 (MQ=255)
                                                 tAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGa                                                                                                   >  2:4630954/1‑150 (MQ=255)
                                                 tAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGa                                                                                                   <  1:4491556/150‑1 (MQ=255)
                                                 tAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGa                                                                                                   >  2:902930/1‑150 (MQ=255)
                                                  aaCTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAg                                                                                                  >  2:104549/1‑150 (MQ=255)
                                                     tGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCtt                                                                                               >  2:1214180/1‑150 (MQ=255)
                                                     tGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCtt                                                                                               <  1:2070474/150‑1 (MQ=255)
                                                     tGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCtt                                                                                               >  2:4333116/1‑150 (MQ=255)
                                                      gATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCttt                                                                                              <  1:2308958/150‑1 (MQ=255)
                                                       atatCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTc                                                                                             >  2:44364/1‑150 (MQ=255)
                                                         atCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTg                                                                                           >  2:1766591/1‑150 (MQ=255)
                                                           ccATCGCCACCTTGATGGCAACATTCGTCCCCAGACCCTTCTTGAACTTGGCCGCCAGTCTAAGATCTCGCTTCCTGCACAATCCCTGGATCCACTGATTCCCCACGTGCAGGTCATTGCCAACGAACCCGCTCTGGTGAGCTTGCTGAc                                                                                         >  2:4753227/1‑150 (MQ=255)
                                                           ccATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGAc                                                                                         <  2:526486/150‑1 (MQ=255)
                                                           ccATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGAc                                                                                         <  1:607040/150‑1 (MQ=255)
                                                             aTCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCa                                                                                       >  2:2251757/1‑150 (MQ=255)
                                                                 ccACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACt                                                                                   <  2:2728429/150‑1 (MQ=255)
                                                                    ccTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGa                                                                                >  1:1628141/1‑150 (MQ=255)
                                                                     cTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGAc                                                                               <  1:397431/150‑1 (MQ=255)
                                                                      ttGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACt                                                                              >  2:1345691/1‑150 (MQ=255)
                                                                          tGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTgggg                                                                          <  1:385037/150‑1 (MQ=255)
                                                                           ggCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGc                                                                         >  2:3636729/1‑150 (MQ=255)
                                                                           ggCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGc                                                                         >  2:534718/1‑150 (MQ=255)
                                                                             cAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGt                                                                       >  2:2734953/1‑150 (MQ=255)
                                                                             cAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGt                                                                       >  1:4350099/1‑150 (MQ=255)
                                                                                cATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTaa                                                                    <  2:3535806/150‑1 (MQ=255)
                                                                                 aTTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTaaa                                                                   >  1:2887859/1‑150 (MQ=255)
                                                                                     gTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTtc                                                               >  2:1473735/1‑150 (MQ=255)
                                                                                      tCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTtct                                                              <  2:3599299/150‑1 (MQ=255)
                                                                                        cccAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCg                                                            >  2:1747104/1‑150 (MQ=255)
                                                                                              ccATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCctctc                                                      >  1:1668607/1‑150 (MQ=255)
                                                                                                  tCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGa                                                  >  1:796484/1‑150 (MQ=255)
                                                                                                   cTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGAt                                                 <  2:3213120/150‑1 (MQ=255)
                                                                                                    ttGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATg                                                <  1:4454435/150‑1 (MQ=255)
                                                                                                    ttGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATg                                                >  1:623445/1‑150 (MQ=255)
                                                                                                         cTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGt                                           <  1:2961410/150‑1 (MQ=255)
                                                                                                         cTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGt                                           <  2:2148122/150‑1 (MQ=255)
                                                                                                          ttGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTc                                          >  2:1110594/1‑150 (MQ=255)
                                                                                                                  cAGTATAATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTg                                  >  2:2747281/1‑150 (MQ=255)
                                                                                                                     tataATATCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCa                               >  1:2732731/1‑150 (MQ=255)
                                                                                                                      ataataTCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCAt                              >  1:2816758/1‑150 (MQ=255)
                                                                                                                       taataTCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCAtt                             >  2:708300/1‑150 (MQ=255)
                                                                                                                          tatCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGa                          <  1:4380396/150‑1 (MQ=255)
                                                                                                                          tatCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTGAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGa                          >  2:2287475/1‑150 (MQ=255)
                                                                                                                          tatCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGATCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGa                          >  2:2565000/1‑150 (MQ=255)
                                                                                                                           atCTCGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGaa                         <  1:4595037/150‑1 (MQ=255)
                                                                                                                             ctcGCTTCCTGCACAATCCCTGGATACACTGATTCCTCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGAGTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGaaaa                       <  2:482742/150‑1 (MQ=255)
                                                                                                                             ctcGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGaaaa                       >  1:3191187/1‑150 (MQ=255)
                                                                                                                             ctcGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGaaaa                       <  2:2872566/150‑1 (MQ=255)
                                                                                                                              tcGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAAc                      >  2:3388961/1‑150 (MQ=255)
                                                                                                                              tcGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAAc                      >  2:3958168/1‑150 (MQ=255)
                                                                                                                               cGCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACa                     <  2:3043631/150‑1 (MQ=255)
                                                                                                                                gCTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACAt                    >  1:4784927/1‑150 (MQ=255)
                                                                                                                                 cTTCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACAtt                   <  1:1855525/150‑1 (MQ=255)
                                                                                                                                   tCCTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGa                 >  2:3980788/1‑150 (MQ=255)
                                                                                                                                    ccTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGaa                <  1:913346/150‑1 (MQ=255)
                                                                                                                                     cTGCACAATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAg               >  1:3916674/1‑150 (MQ=255)
                                                                                                                                        cacaATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATg            >  1:1575805/1‑150 (MQ=255)
                                                                                                                                        cacaATCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATg            >  1:1148186/1‑150 (MQ=255)
                                                                                                                                           aaTCCCTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAg         <  1:4113979/150‑1 (MQ=255)
                                                                                                                                              cccTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGccc      >  1:3095919/1‑150 (MQ=255)
                                                                                                                                              cccTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGccc      >  2:3917537/1‑150 (MQ=255)
                                                                                                                                              cccTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGccc      >  2:2903490/1‑150 (MQ=255)
                                                                                                                                                cTGGATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGt    >  2:1989162/1‑150 (MQ=255)
                                                                                                                                                  ggATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCCGTTGAAAACATTGAAGATGCAGCCCGTCa  >  2:1190650/1‑150 (MQ=255)
                                                                                                                                                  ggATACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCa  >  1:1060126/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TTTTATATCCAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGAACTTGGCCGCCAGTATAATATCTCGCTTCCTGCACAATCCCTGGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCA  >  CP009273/1696463‑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: ‑